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^7  The  Applied  Physics  Laboratory  isi  a division  of  The  Johns  Hopkins 
University,  dedicated  to  serving  the  4iation  by  performing  research  and 
development  in  important  areas  of  National  defense  and  welfare.  APL  is 
located  on  a 35&-acre  site  in  HowarcJ  County,  Maryland.  About  one-half  of 
it  2500-member  staff  are  scientists  ai>d  engineers. 

'4 

The  principal  mission  of  APL  is  "to  further  national  defense  through 
the  application  of  science  and  technology,  primarily  for  the  Department  of 
the  Navy.  With  the  encouragement  of  the  Department  of  Defense,  APL  also 
undertakes  on  a noninterfering  basis  to  apply  its  technical  capabilities  to 
nationally  important  nondefense  areas  to  which  it  can  make  a unique 
contribution.  ~ 

APL  was  established  by  The  Johns  Hopkins  University  in  1942  at  the 
urgent  request  of  the  Office  of  Scientific  Research  and  Development,  to  imple- 
ment  the  concept  that  science  and  technology  would  more  effectively  further 
national  defense  by  addressing  broad  tactical  situations  and  requirements 
rather  than  by  limiting  efforts  to  specific  technical  objectives.  APL  was 
challenged  by  the  assignment  of  applying  new  technology  to  the  defense  of 
the  Fleet  against  air  attack.  The  radio  proximity  fuse,  whicf'  vvas  developed 
and  produced  in  one  year,  made  a major  c'u.  on  to  American  victories 
in  Europe  and  the  Pacific. 

After  World  War  II,  the  University,  at  i ng  of  the  Navy,  agreed  to 
continue  its  sponsorship  of  APL  as  a unique  national  resource  with  a broad 
mission  to  enhance  the  security  of  the  Fleet  through  application  of  new 
science  and  technology.  This  relationship  has  proved  fruitful  in  successful 
development  and  deployment  of  guided  missile  systems.  Fleet  defense  combat 
systems,  sea-based  strategic  deterrent  systems,  aircraft  defense,  a precise 
worldwide  positioning  and  navigation  system,  and  other  areas  in  which  APL 
has  made  important  technical  contributions. 

From  the  beginning,  technical  advances  in  defense  systems  led  to  spinoffs 
in  the  civil  area,  for  example,  high-speed  aircraft  and  electronics.  However, 
it  became  increasingly  apparent  that  these  new  techniques  and  capabilities 
could  be  advantageously  applied  more  directly  to  important  civilian  areas. 
APL  conceived,  developed,  and  initially  operated  the  Navy  Navigation  Satel- 


lite  System  (Transit),  which  supplies  precise  worldwide  position  and  naviga- 
tion data  by  means  of  a network  of  earth-orbiting  satellites.  A simplified 
version  of  this  system  is  now  widely  used  by  commercial  shipping.  The 
Transit  System  has  produced  precise  worldwide  geodetic  measurements  and 
basic  new  knowledge  of  the  shape,  gravitational  field,  and  motion  of  the 
earth.  Under  NASA  sponsorship,  APL  has  built  and  flown  numerous  scientific 
satellites  and  space  experiments  that  have  helped  to  revolutionize  our  knowl- 
edge and  understanding  of  the  earth,  the  solar  system,  and  the  stellar  universe. 

Other  civil  areas  in  which  APL  is  advantageously  applying  technology 
developed  in  defense  programs  include  civil  air  safety  and  traffic  control, 
harbor  and  inland  waterway  traffic  management,  automated  urban  transit  sys- 
tems, the  environmental  impact  of  electric  power  plants  and  other  energy 
facilities,  fire  prevention  and  control,  and  the  development  of  new  energy 
sources  and  systems. 

A highly  successful  collaborative  program  between  APL  and  the  Johns 
Hopkins  Medical  Institutions  is  conducted  in  the  field  of  biomedical  engineer- 
ing. Technology  developed  in  space,  electronics,  and  computer  application 
programs  is  being  successfully  applied  to  medical  research,  techniques,  and 
instruments  for  diagnosis,  therapy,  and  improved  health  care  delivery  systems. 

To  gain  the  insights  required  for  future  technical  advances  and  as  a part 
of  its  function  as  a University  laboratory,  APL  conducts  a vigorous  research 
program  in  the  basic  science;  underlying  its  evolving  technology.  The  quality 
of  this  work  is  attested  by  its  recognition  in  scientific  publications,  honors, 
and  awards. 

A Center  of  the  Evening  College  of  the  University  located  at  APL  serves 
the  Washington-Baltimore  area  by  offering  accredited  graduate  courses  in 
science  and  technology.  Most  of  the  faculty  are  drawn  from  the  APL  Senior 
Staff.  Other  areas  of  scholarly  cooperation  between  APL  and  other  divisions 
of  the  University  include  predoctoral  and  postdoctoral  fellowships,  visiting 
professorships,  joint  appointments,  joint  seminars,  and  collaborative  research. 

^ This  volume  presents  a representative  sample  of  the  civilian-oriented  and 
basic  research  programs  that  currently  comprise  about  20%  of  the  total  APL 
effort.  All  results  of  these  programs  are  fully  available  to  the  public. 
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INTRODUCTION 


The  APL  program  in  space  science  and  technology  began  in  the  immediate 
post  World  War  II  perirxl  when  a group  then  led  by  Dr.  James  Van  Allen 
started  to  study  cosmic  radiation  and  the  upper  atmosphere  using  instrumented 
sounding  rockets.  The  discovery  in  1957  by  APL  scientists  that  the  position 
of  the  Russian  Sputnik  satellite  could  be  tracked  by  observing  the  droppler 
shift  of  its  radio  signal  led  to  the  conception  of  a satellite  navigation  sys- 
tem based  on  the  doppler  tracking  principle.  The  Navy  Navigation  Satellite 
System  (Transit)  developed  by  APL  has  provided  precise,  worldwide,  all- 
weather  navigational  data  to  the  Navy  since  1964  and  to  commercial  users 
since  1967.  It  includes  a constellation  of  satellites;  tracking,  receiving,  and 
computing  stations;  equipment  and  software;  and  user  receiving  and  comput- 
ing equipment. 

The  Transit  program  has  been  continually  upgraded  since  its  inception. 
Recent  improvements  include  a single-axis  active  disturbance  compensation 
system  (DISCOS)  to  provide  improved  orbit  prediction  and  improved  hybrid 
microelectronic  circuitry  packaging.  These  are  described  in  this  section.  Other 
improvements,  not  discussed  in  detail,  yield  increased  satellite  and  ground 
equipment  reliability  and  greater  navigational  accuracy. 

The  Transit  program  required  more  accurate  knowledge  of  the  shape  of 
the  earth  for  navigational  accuracy  and  provided  the  measurement  tool  for 
acquiring  this  knowledge.  The  comprehensive  geodesy  program  conducted  by 
APL  with  the  Transit  System  improved  the  accuracy  of  the  geoid  model  by 
about  two  orders  of  magnitude  and  led  to  APL  participation  in  other  space 
science  programs.  Articles  in  this  section  describing  recent  experiments  on  the 
interactions  of  charged  particles  with  the  earth’s  magnetic  field  are  representa- 
tive of  APL’s  continuing  space  physics  research. 

APL’s  special  capabilities  in  attitude  stabilization  and  control,  ultrastable 
oscillators,  precision  tracking  and  other  technologies,  and  component  and 
system  design  and  development  are  used  in  numerous  spare  programs  for 
NASA  and  DoD  as  well  as  the  Navy.  Particularly  notable  at  this  time  are  the 
latest  members  of  a series  of  small  astronomy  satellites  (SAS)  and  geodetic 
research  satellites  (GEOS)  that  APL  has  designed,  built,  and  launched  for 
NASA.  'The  third  members  of  each  series,  SAS-C  and  GEOS-C,  were  com- 
pleted and  launched  this  year.  Both  have  been  highly  successful  and  are  add- 
ing significant  new  knowledge  in  their  areas  of  research.  Key  design  features, 
launch,  and  early  postlaunch  operation  of  these  scientific  satellites  are  de- 
scribed in  this  section. 

Rain  attenuation  may  be  a problem  in  high-frequency  satellite-to-earth 
conununication.  Techniques  to  predict  rain  attenuation  statistics  are  described. 
Current  programs  to  improve  time  and  frequency  standards  and  the  automated 
receipt,  storage,  and  handling  of  high-speed  data  are  also  representative  of 
the  many  APL  activities  contributing  to  increased  capabilities  in  space  science 
and  technology. 


DISCOS  SYSTEM  FOR  THE  TIP-II  SATELLITE 


A Disturbance  Compensation  System  (DISCOS)  to 
correct  for  drag  effects  has  been  incorporated  in  the 
Transit  Improvement  Program  (TIP-II)  Navy  Navigation 
Satellite.  The  system  wiil  force  the  spacecraft  to 
follow  a predictable  orbit,  aliowing  the  generation  of 
more  accurate  ephemeris  data  for  navigation  fixes 
over  a longer  period  of  time. 


BACKGROUND 

The  present  Transit  satellites  require  ephemeris 
updating  at  12-hour  intervals  because  of  the  disturb- 
ing effects  of  solar  radiation  pressure  and  atmospheric 
drag.  To  produce  a longer  period  of  accurate  ephem- 
eris, DISCOS  has  been  incorporated  in  the  TIP 
spacecraft. 

The  first  DISCOS,  installed  in  Triad  (TIP-I), 
compensated  for  disturbances  in  all  three  axes  (Ref. 
1).  A spherical  "proof  mass”  in  an  enclosed  chamber 
followed  a nearly  pure  gravitational  orbit.  Its  position 
was  sensed,  and  cold-gas  thrusters  positioned  the  satel- 
lite about  the  proof  mass,  forcing  the  satellite  to 
follow  a "drag-free”  trajectory.  Although  the  three- 
axis  DISCOS  produced  an  exceptionally  accurate  orbit 
(Ref.  2),  Triad  was  a difficult  spacecraft  to  fabricate. 
The  gravitational  force  between  the  proof  mass  and 
all  components  in  Triad  was  disturbing  enough  so 
that  bookkeeping  of  the  mass  and  of  the  position  of 
each  component  was  necessary  to  compensate  for  mass 
attraction  (Ref.  3). 

The  single-axis  DISCOS  used  on  TIP-II  was 
devised  to  achieve  a drag-free  satellite  in  a simpler 


fashion.  A cylindrical  proof  mass  was  electromag- 
netically  suspended  on  an  axis  oriented  in  the  direction 
of  the  spacecraft  velocity  vector.  While  free  to  move 
without  friction  along  the  axis,  the  proof  mass  was 
constrained  in  radial  motion  by  the  suspension  force. 
In  this  way,  disturbance  forces  affecting  the  spacecraft 
velocity  could  be  compensated  for  by  firing  a set  of 
thrusters  while  the  internal  mass  attraction  force  was 
nearly  normal  to  the  axis,  a direction  opposed  by  the 
suspension  mechanism.  An  incrementally  command- 
able  biasing  system  provided  trimming  of  any  re- 
sidual force  components.  This  reduced  the  motion 
of  the  proof  mass  to  a gravitational  trajectory  without 
the  need  for  mass  attraction  bookkeeping,  thereby 
simplifying  fabrication  of  the  spacecraft. 


DISCUSSION 

The  Navy  Navigation  Satellite  Program  required 
the  capability  of  maintaining  the  projected  orbital 
ephemeris  with  an  accuracy  of  ±85  m for  a seven-day 
period.  This  required  that  long-term  acceleration  in 
the  along-track  direction  be  compensated  to  within 
1.58  X 10^’*  g (1  g = acceleration  of  gravity  = 
980  cm/s’).  To  provide  this  capability,  the  single- 
axis DISCOS  was  designed  and  fabricated  at  APL. 
The  system  consists  of  a sensor  pjositioned  at  the  space- 
craft center  of  mass  and  connected  via  cabling  to  an 
electronics  package  located  in  the  spacecraft  main 
body.  A set  of  Teflon  solid  propellant  thrusters  devel- 
oped at  Fairchild  Industries  provides  the  thrust  neces- 
sary to  counteract  disturbance  forces. 
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Fig.  1 Cutaway  view  of  siagle-axis  DISCOS  sensor. 


To  meet  a design  criterion  of  along-track  force  com- 
pensation to  within  10~"  g,  the  sensor  was  configured 
as  a cylindrically  symmetric  enclosure,  shielding  the 
proof  mass  from  external  forces  while  producing  little 
mass  attraction  along  the  axis  (Fig.  1).  A thermal 
control  system  stabilizes  the  sensor  temperature  to 
within  ±0.5°  of  its  nominal  70°  ± 4°F  set  point. 
Thermal  gradients  of  less  than  ±2°F  within  the 
sensor  neutralize  radiometric  pressure  differentials  on 
the  ends  of  the  proof  mass. 

The  proof  mass,  formed  into  a cylindrical  shell,  is 
the  sensed  element  that  follows  a purely  gravitational 
trajectory.  It  is  made  of  pure  aluminum  with  low 
magnetic  susceptibility,  and  its  surface  is  vacuum 
deposited  with  2000  nm  of  gold  to  distribute  electro- 
static potential  evenly.  A caging  system  protects  the 
proof  mass  during  launch  and  releases  it  during 
gravity-gradient  boom  deployment. 

The  frictionless  levitation  force  is  provided  by  a 
2.08-kHz  square-wave  current  passing  through  the 
axis,  generating  eddy  current  repulsion  forces  in  the 
proof  mass  (Ref.  4).  Commandable  to  four  levels,  the 
suspension  current  produces  forces  on  the  order  of 
10^'^  g to  balance  cross-track  disturbances.  As  shown 
in  Fig.  2 for  a sine-wave  current,  the  suspension  force 
constant  is  a function  of  current  frequency  and  radial 
offset  of  the  proof  mass  on  its  axis.  By  using  a square 
wave,  the  force  is  increased  by  a factor  of  2.19.  A 
proof  mass  weight  of  12.3  g produces  an  acceleration 
constant  of  6.72  X 10~“  g/A’'.  The  suspension  cur- 
rent is  returned  through  evenly  spaced  radial  conduc- 
tors in  the  end  caps  and  resistive  wires  parallel  to  the 
axis,  preventing  imaging  effects  from  biasing  the  proof 
mass. 


PROOF  MASS  RADIAL  DISPLACEMENT 
(center  offset/cylindef  radius) 


Fig.  2 Theoretical  eddy  current  suspension  Force  for 
sine-wave  excitation  of  aluminum  proof  mass  (1.04-cm 
ID  by  1.70-cm  OD). 


Two  optical  systems  are  used  to  detect  proof-mass 
position.  An  axial  system  determines  along-track  posi- 
tion, and  a transverse  system  detects  radial  offset  about 
the  axis.  The  axial  optical  system  uses  parabolic  mir- 
rors to  collimate  light  from  a 100  000-hour-life  tung- 
sten bulb,  pass  the  light  across  the  ends  of  the  proof 
mass,  and  refocus  it  on  a set  of  photocells.  The  differ- 
ence in  current  produced  between  the  two  photocells 
is  a linear  measure  of  proof-mass  position  along  the 
axis.  Both  coarse- range  (±10  mm)  and  fine-range 
(±2  mm)  proof-mass  positions  are  telemetered  to 
ground  receivers,  providing  position  resolution  of  0.08 
and  0.015  mm,  respectively.  The  data  are  used  in 
aeronomy  studies  and  to  provide  in-orbit  thruster 
calibration. 

The  transverse  system  uses  the  light  from  the  axial 
system  bulbs.  Spherical  mirrors  at  the  axis  ends  col- 


Fig.  3 Block  diagram  showing  TIP-II  single-axis 
DISCOS  function. 
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limate  the  bulb  light,  pass  it  through  the  proof -mass 
center,  and  refocus  it  on  a photocell.  This  signal,  a 
measure  of  radial  position  of  the  proof  mass,  is  differ- 
entiated to  determine  radial  velocity.  The  result  is  used 
to  modulate  the  suspension  current,  providing  active 
damping  of  proof -mass  dynamics. 

The  axial  position  signal  is  the  input  to  the  derived 
rate  modulator  (DRM)  control  electronics,  which  is 
the  plant  controller  of  the  system.  Relative  motion  of 
the  proof  mass  is  analyzed  by  the  DRM  and,  allowing 
for  hysteresis  to  provide  damping  for  system  stability, 
signals  are  generated  to  fire  the  thrusters  at  the  proper 
time. 

The  thrusters  consist  of  a pair  of  opposing  Teflon 
pulse  plasma  engines,  each  of  which  produces  an  im- 
pulse of  44  dyn-s.  Upon  command  of  the  DRM,  a 
high-voltage  discharge  is  initiated  across  the  end  sur- 
face of  a Teflon  stick,  ionizing  particles  that  are  then 
accelerated  by  an  electromagnetic  field  to  produce  the 
impulse.  This  thrust  reacts  against  the  spacecraft  to 
counteract  disturbances  (Fig.  3). 

A biasing  system  is  required  to  balance  any  axial 
components  of  proof-mass  disturbance  forces.  The 
system  consists  of  two  opposing  coils,  each  configured 


to  provide  a constant  force  over  a proof -mass  range  of 
± 10  mm.  Over  the  normal  operating  range  of  ±2.25 
mm,  a variation  of  less  than  0.05%  exists  in  the  bias- 
ing force.  Again,  a 1 .04-kHz  sine-wave  current  pass- 
ing through  the  selected  coil  produces  eddy  current 
repulsion  of  the  proof  mass.  Coil  current  is  command- 
able  to  a maximum  of  10— » g,  in  10— ’*-g  increments, 
providing  in-orbit  tailoring  of  the  proof -mass  velocity. 
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HYBRID  MICROCIRCUITS  FOR  TIP-11 
AND  TIP-III  SATELLITES 


The  Transit  Improvement  Program  Navy  Navigation 
Satellite  designs  (TIP-II  and  TIP-III)  required  ad- 
vanced electronic  packaging  techniques  In  order  to 
meet  the  goafs  ot  size,  weight,  power  consumption, 
and  reliability.  Through  the  effective  use  of  hybrid 
microelectronic  packaging,  APL  has  assisted  the  TIP 
system,  circuit,  mechanical,  and  packaging  designers 
In  meeting  their  goals. 


BACKGROUND 

An  APL  hybrid  microcircuit  for  spacecraft  use 
typically  consists  of  a miniature  ceramic  substrate 
that  supports  deposited  aluminum  conductors,  de- 
posited chromium  resistors,  and  deposited  silicon 
monoxide  insulators.  Such  active  elements  as  trans- 
sistors,  diodes,  and  integrated  circuits  or  such  special 
passive  elements  as  capacitors,  inductors,  and  high- 
16  value  resistors  are  attached  using  eutectic  soldering. 


special  epoxies,  or  solder  reflow  techniques.  Tliese 
elements  are  "wired  in”  using  thermocompression 
wire  bonding  or  ultrasonic  wire  bonding  techniques. 
Then  the  assembly  is  sealed  hermetically  in  an  inert 
atmosphere.  Rigorous  electrical  testing  and  reliability 
screening  ensure  a high-quality  device  for  the  space 
environment. 

DISCUSSION 

APL  has  designed,  assembled,  tested,  and  screened 
more  than  1000  hybrid  microcircuits  used  in  the  com- 
mand system,  telemetry  system,  computer/memory, 
and  attitude  stabilization  system  of  the  TIP-II  and 
TIP-III  satellites.  In  addition,  assemblies  for  the 
command  receiver,  5-MHz  oscillator,  400-MHz 
transmitter,  150-MHz  transmitter,  and  frequency 
multiplier/phase  modulator  system  and  approximately 
100  miniature  circuit  substrates  were  assembled. 


Figures  1 and  2 illustrate  a relatively  simple  ap- 
plication of  the  technique.  The  high-current  core- 
memory driver  circuit  shown  schematically  in  Fig. 
1 was  converted  to  microcircuit  form.  If  the  circuits 
were  built  using  miniature  discrete  components  on 
printed  circuit  boards,  the  X-Y  memory  driver  shown 
might  occupy  more  than  1 sq.  in.  of  board  area  and 
require  0.30  in.  of  board  height.  The  hybrid  shown 
in  Fig.  2 occupies  less  than  0.20  sq.  in.  of  board  area 
and  less  than  0.080  in.  of  board  height.  The  signifi- 
cant features  of  the  technique  are  identified  in  Fig.  2. 

A more  complex  example  of  the  technology,  the 
digit  driver/differential  comparator  used  in  the  TIP 
memory,  is  depicted  in  Figs.  3,  4,  and  5.  The  sche- 
matics of  four  digit  drivers  (Fig.  4)  and  four  differ- 
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Fig.  I Schematic  of  X-Y  memory  driver. 


Fig.  2 X-Y  memory  driver  hybrid  microcircuit. 
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Fig.  6 TIP  frequency  multipliei/phase  modulator. 


ential  comparators  (Fig.  4)  have  been  incorporated 
into  a single  1.23  by  1.23  in.  flat  package  (Fig.  3) 
comprising  16  transistors,  20  diodes,  68  resistors, 
eight  capacitors,  and  four  monolithic  integrated  cir- 
cuits. More  conventional  technology  would  require 
approximately  10  sq.  in.  compared  with  about  2.0 
sq.  in.  to  perform  this  function. 

Other  custom  miniature  circuits  and  assemblies  have 
been  developed  for  TIP  using  hybrid  technology. 
Typically  these  circuits  are  slightly  larger  than  the 
hermetic  hybrids  and  are  built  on  ceramic  substrates. 
Conductor  patterns  are  defined  by  etching  chromium/ 
gold  conductors;  miniature  uncased  components  are 
added  separately  Using  soHer  reflow  techniques.  Fig- 
ure 6 illustrates  the  TIP  frequency  multiplier/phase 
modulator  that  has  been  assembled  using  diese 


techniques.  The  technology  usually  is  applied  to  RF 
circuitry  where  components  such  as  inductors,  trans- 
formers, capacitors,  and  strip-line  or  microstrip  ele- 
ments are  required.  Circuits  with  capabilities  up  to 
several  gigahertz  have  been  fabricated  successfully. 

The  complex  and  sophisticated  TIP  satellite  elec- 
tronics design  and  rigorous  reliability  requirements 
were  successfully  acconunodated  within  the  available 
volume  and  weight  constraints.  This  could  not  have 
been  done  without  the  use  of  the  hybrid  microcircuit 
fabrication  and  padcaging  techniques. 
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CHARACTERISTIC  ENERGY  SPECTRA 
OF  1 TO  500  eV  ELECTRONS  OBSERVED 
IN  THE  AURORAL  ZONES  FROM 
ATMOSPHERE  EXPLORER-C 


Data  obtained  from  the  photoelectron  experiment 
on  board  the  Atmosphere  Explorer-C  (AE-C)  satellite 
have  been  used  in  what  is  believed  to  be  the  first 
general  classification  of  low-energy  (1  to  500  eV)  elec- 
tron fluxes  in  the  high-latitude  ionospheric  regions 
with  high-resolution  energy  spectra  measurements.  In 
addition  to  the  energy  spectra  already  associated  with 
the  normal  day  airglow  (Ref.  1),  four  major  c/asses  of 
low-energy  electrons  have  been  identified.  Three  of 
these  are  lalatively  stable  while  the  fourth  has  large 
temporal  variations. 


BACKGROUND 

The  joint  APL/JHU  Department  of  Qiemistty 
photoelectron  experiment  on  AE-C  was  designed  pri- 
marily to  measure  the  energy  distribution  of  photo- 
electrons in  the  daytime  ionosphere  (Ref.  2),  The 
instrument  was  designed  to  detect  electrons  with  ex- 
tremely low  energies  (down  to  1 eV)  and  to  have 
very-high-energy  resolution  (AE/E  = 2.5%),  with 
on-board  tape  recorders  providing  continuous  data 
collection. 

The  experiment  and  its  mounting  on  the  space- 
craft have  previously  been  described  (Ref.  3).  Brief- 
ly, two  electrostatic  analyzer  heads  are  mounted  on 
opposite  sides  of  one  end  of  the  cylindrical  space- 
craft. The  spacecraft  operates  in  either  a despun  or 
a 4-rpm  spinning  mode.  In  the  despun  mode,  one 
sensor  head  looks  up  away  from  the  earth  while  the 
other  looks  down  toward  the  earth.  There  is  a large 
variety  of  possible  modes  of  operation,  including 
variable  sweep  rate  and  energy  scan  range. 

The  results  presented  here  were  determined  from 
analysis  of  the  data  obtained  in  the  0 to  500  eV  high- 
sweep-rate  mode.  In  this  mode,  one  16-point  spec- 
trum is  produced  every  250  ms.  The  16  ener^  bins 
are  distributed  linearly  along  the  0 to  500  eV  sweep 
range. 


the  spacecraft  in  a circular,  250-km  orbit.  The  con- 
stant altitude  near  the  peak  of  the  F region  of  the 
ionosphere  has  also  been  extremely  helpful  in  sorting 
electron  precipitation  phenomena. 


DISCUSSION 

Although  the  high-latitude  fluxes  of  1 to  500  eV 
electrons  observed  from  AE-C  are  highly  variable, 
it  has  been  possible  to  organize  the  observed  events 
into  four  major  classes : 

1.  Fluxes  of  electrons  having  a peak  in  energy 
below  100  eV,  approximately  isotropic  over 
the  upper  2w  hemisphere,  and  located  at  high 
invariant  latitudes  on  the  dayside,  that  are  ap- 
parently gaining  access  to  the  ionosphere 
through  the  dayside  magnetospheric  cusps; 

2.  Structured  fluxes  of  electronr  showing  one  or 
more  discrete  but  variable  peaks  in  energy, 
reminiscent  of  "inverted-V”  p^tenomena  at 
higher  energies,  observed  at  all  magnetic  local 
times; 

3.  FIuxk  that  show  a monotonically  decreasing 
differential  flux  versus  energy  dependence  pro- 
portional to  where  £ > 30  eV  and  1 < 
o < 2,  characteristic  of  auroral  secondary 
electrons;  and 

4.  Intense  fluxes  of  electrons  with  rapidly  dung- 
ing and  random  energy  spectra,  characteristic 
of  passage  through  auroral  arcs  or  breakup 
regions. 

Fluxes  of  types  3 and  4 depend  on  predpi- 
tation  of  high-energy  auroral  primary  particles  while 
those  of  types  1 and  2 appear  to  be  a conse- 
quence of  fundamental  processes  that  couple  the 
magnetosphere  and  ionosp^re. 
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The  68®  inclination  of  the  AE-C  orbit  takes  the 
spacecraft  to  geomagnetic  latitudes  as  high  as  80°, 
depending  on  longitude.  Although  the  AE-C  orbit 
was  not  chosen  primarily  for  hi^-latitude  research, 
it  has  been  particularly  effective  since  the  spacecraft 
remains  in  the  auroral  region  for  long  periods  of  time. 
These  passes  have  been  extremely  useful  for  the  ob- 
servation of  relatively  stable  electron  precipitation 
events.  Most  of  the  data  analyzed  were  taken  with 


Figures  1 and  2 show  the  spatial  locations  in  the 
north  polar  region  where  the  electron  fluxes  denoted 
as  class  1 and  2,  respectively,  have  been  observed  by 
AE-C.  The  regions  are  indicated  along  each  orbit 
where  class  1 and  2 fluxes  have  been  observed,  cor- 
responding to  geomagnetic  conditions  characterized 
by  Kp  < 3 (quiet  to  moderately  disturbed)  and  Kp 
> 3 (disturbed).  Because  of  its  inclination,  AE-C 
could  not  reach  geomagnetic  latitudes  higher  than 
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Fig.  1 Orbital  tracks  of  AE— C over  the  northern  high- 
latitude  region  with  the  location  of  cusp  (class  1)  elw- 
trons  for  Kp  < 3 and  Kp  ^ 3. 
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Fig.  2 Orbital  tracks  of  AE-C  over  the  northern  high-  I 

latitude  region  with  the  location  of  structured  (class  2) 

electrons  for  Kp  <i  and  Kp  > 3.  ' 


81°;  the  "inaccessible  region”  is  shaded  in  Figs.  1 
and  2. 

Figures  1 and  2 demonstrate  that  the  class  1 elec- 
tron fluxes  occur  most  often  in  the  dayside  high-lati- 
tude region  associated  with  cusp  phenomena,  whereas 
the  class  2 electrons  are  found  at  nearly  all  local 
times  and  during  almost  every  AE-C  orbit  shown. 
The  distributions  support  the  view  that  the  class  1 
electrons  have  gained  direct  access  through  the  mag- 
netospheric  cusp  from  the  magnetosheath/solar  wind 
plasma  and  that  the  class  2 electrons  have  been  ac- 
celerated and  redistributed  in  the  region  between  the 
source  and  the  auroral  ionosphere. 

The  AE-C  measurements  are  an  extension  to  low 
energies,  lower  altitudes,  and  higher-energy  resolution 
compared  with  previous  experiments  in  the  high-lati- 
tude  region.  These  more  detailed  data  will  help  iden- 
tify the  important  magnetospheric/ionospheric  regions 
and  coupling  mechanisms  and  will  help  resolve  cur- 


rent uncertainties  concerning  the  location  of  the  cusp 
region,  the  nature  of  daytime  auroral  phenomena,  and 
the  source  and  charge  carriers  of  field-aligned  cur- 
rents. 
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FIELD-ALIGNED  CURRENTS  IN 
THE  SOUTHERN  AURORAL  ZONE 
MEASURED  BY  TRIAD 


Auroral  phenomena  have  been  studied  in  the  north 
poiar  regions  since  the  iast  century  by  means  of  their 
dramatic  optical  emissions.  Similar  auroral  phenomena 
have  long  been  expected  to  occur  in  the  southern 
polar  regions,  but  because  of  the  remote  location  of 
appropriate  Antarctic  iand  areas  with  respect  to  popu- 
lated areas,  and  the  very  severe  weather  conditions, 
very  few  measurements  have  been  made  there.  With 
the  establishment  of  a satellite  ground  station  at 
McMurdo,  Antarctica,  in  March  1974,  data  obtained 
from  the  magnetometer  experiment  on  board  the 
Navy/APL  Triad  satellite  are  providing  the  first  con- 
tinuous mapping  of  fieid-aiigned  currents  and  asso- 
ciated phenomena  in  the  southern  auroral  zone. 


BACKGROUND 

The  concept  of  field-aligned  currents  in  the  auroral 
region  was  introduced  by  Birkeland  in  1908  but  re- 
mained a subject  of  controversy  until  1966  when  the 
single-axis  magnetometer  measurements  from  the 
Navy/APL  satellite  1963-38C  provided  the  first 
direct  evidence  for  their  existence  in  the  northern 
polar  regions  (Ref.  1 ) . 

The  first  vector  magnetic  field  observations  in  tiic 
northern  polar  regions  were  made  by  the  Navy/APL 
Triad  satellite,  which  was  launched  into  a circular 
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Fig.  1 Image  of  visual  auroral  arcs  in  the  southern  auroral  zone  obtained  by  the  DAPP  satelUte,  with  the 
orbit  of  the  Triad  satellite  superimposed.  The  magnetic  field  recorded  by  one  of  the  magnetometers  on  board 
Triad  is  plotted  on  the  left.  The  south  magnetic  (dipole)  pole  is  shown  at  the  top  of  the  figure. 


polar  orbit  at  800-km  altitude  on  2 September  1972. 
The  obsen’ation.'i  were  obtained  w'ith  an  on-board 
three-axis  flux-j>ate  magnetometer  at  a rate  of  2.25 
samples  per  axis  per  second  with  a resolution  of  12  y 
(where  the  total  geomagnetic  field  intensity  is  approx- 
imately 0.5  gauss,  or  50  000  y).  The  data,  obtained 
from  a satellite  ground  station  in  College,  Alaska, 
have  been  used  in  several  studies  of  auroral  phenom- 
ena in  the  northern  polar  region  (Refs.  2 through  6). 

Data  recording  was  begun  at  McMurdo,  Antarctica, 
in  March  1974  and  continues  to  the  present  time. 
These  data  have  been  used  for  the  first  time  to  study 
the  amplitudes,  directions,  and  spatial  distribution 
of  field-aligned  currents  in  the  southern  auroral  zone 
(Ref.  7).  The  information  accumulated  by  the  Triad 


magnetometer  experiment  in  the  north  and  south  au- 
roral zones  is  expected  to  play  a very  important  role 
in  future  studies  of  the  major  complex  mechanisms 
that  couple  energy  from  the  sun,  through  the  mag- 
netosphere, to  the  atmosphere  and  ionosphere.  The 
importance  of  the  field-aligned  currents  can  be  appre- 
ciated by  the  fact  that  their  total  intensit)-  varies  be- 
tween 1 million  and  10  million  amperes. 


DISCUSSION 

The  intensity  and  direction  of  currents,  J,  in  the 
auroral  region  are  determined  from  the  Triad  mag- 
netic field  measurements,  B,  with  the  Maxwell  equa- 
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Fig  2 Image  of  visual  auroral  features  and  Australian  city  lights  obtained  by  the  DAPP  satellite,  with  the  orbit 
of  the  Triad  satellite  superimposed.  TTie  magnetic  field  recorded  by  one  of  the  Triad  magnetometers  is  plotted 
on  the  left.  The  south  polar  region  is  located  at  the  bottom  of  this  figure. 
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tion  J = — V X B.  Therefore  magnetic  distur- 

flo 

bances  observed  by  Triad  indicate  the  presence  of 
electric  currents  that,  in  some  cases,  may  be  due  to 
streams  of  energetic  particles  bombarding  the  auroral 
atmosphere  and  producing  optical  auroral  forms.  Fig- 
ure 1 shows  an  example  in  which  a pronounced  mag- 
netic disturbance  was  observed  by  Triad  over  a south 
polar  region  where  a visual  auroral  arc  was  concur- 
rently observed  with  the  optical  imaging  system  on 
board  the  U.S.  Air  Force  DAPP  satellite.  In  this 
example,  a field-aligned  current  equal  to  2 jxA/m‘  was 
observed  by  Triad  near  magnetic  local  midnight  (at 
0030  Magnetic  Local  Time  (MLT))  in  the  same 
region  where  the  DAPP  satellite  had  observed  a 
strong  auroral  arc  approximately  18  minutes  earlier. 

Figure  2 shows  a similar  example  of  aurora,  with 
the  lights  of  major  cities  in  Australia  also  appearing. 
In  this  case,  the  magnetic  disturbance  detected  by 
Triad  indicates  that  a pair  of  upward  and  downward 
current  sheets  with  intensities  of  approximately 
3 /tA/m®  occurred  In  a region  where  a large  auroral 
form  was  photographed  by  the  DAPP  satellite  ap- 
proximately 30  minutes  earlier. 

Not  all  current  systems  detected  with  the  Triad 
satellite  can  be  associated  with  visual  auroral  features, 
principally  because  of  the  threshold  energy  of  the 
current  particles  required  to  produce  optical  emis- 


sions. Consequently,  besides  providing  the  only  con- 
tinuous vector  magnetic  field  measurements  in  both 
hemispheres,  the  Triad  satellite  is  the  most  sensitive 
means  of  detecting  the  energy  input  to  the  earth’s 
auroral  zones. 
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COMPLETION  AND  LAUNCHING  OF 
SMALL  ASTRONOMY  SATELLITE-C 


The  Small  Astronomy  SaMlite-C  (SAS-C)  was 
launched  by  a Scout  rocket  on  7 May  1975  from  the 
San  Marco  Equatorial  Range,  Kenya,  East  Africa,  and 
Is  performing  well  in  orbit.  The  third  astronomy 
satellite  developed  by  APL  for  NASA,  SAS-C  carries 
an  X-ray  experiment  designed  by  The  Massachusetts 
Institute  of  Technology.  The  design,  fabrication,  and 
testing  of  SAS-C  (Fig.  1)  extended  over  a 3>/i-year 
period  and  culminated  in  a successful  astronomy 
mission  to  explore  the  X-ray  sources  of  the  celestial 
sphere. 


BACKGROUND 

SAS-C  completes  the  series  of  three  astronomy 
satellites  built  by  APL  for  NASA  Goddard  Space 


Flight  Center  and  launched  from  the  San  Marco 
Range.  After  achieving  the  desired  orbit,  SAS-C  was 
designated  SAS-3  by  NASA.  SAS-A,  also  called 
Uhuru,  was  launched  on  12  December  1970.  It  car- 
ried an  X-ray  detector  experiment  made  by  American 
Science  and  Engineering  Company.  SAS-B,  launched 
on  16  November  1972,  carried  a gamma-ray  telescope 
developed  by  Goddard.  Much  important  new  informa- 
tion about  the  celestial  sphere  was  obtained  from  these 
satellite  experiments. 

'The  areas  of  X-ray  astronomy  being  investigated 
by  SAS-C  include  the  location  of  X-ray  sources  to 
1 3 arc  seconds,  the  existence  and  identification  of  very 
weak  extragalactic  sources,  the  properties  of  transient 
X-ray  phenomena  associated  with  novas  and  super- 
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Fig.  1 SA&-C  being  tested  by  APL. 


novas,  the  absorption  of  the  low-energy  diffuse  X-ray 
background  by  interstellar  matter,  the  long-  and  short- 
term variations  of  Scorpius  X-1,  the  detailed  energy 
spectrum  of  X-ray  sources  and  the  background  from 
0.1  to  50  keV,  and  the  periodic  time  variations  of  X- 
ray  sources  (e.g.,.  X-ray  pulsars)  with  periods  from 
10"**  to  10*  s.  The  instrumentation  for  the  experiment 
was  designed  and  fabricated  by  the  Center  for  Space 
Research  of  MIT. 

The  spacecraft  control  section,  which  was  designed, 
fabricated,  and  tested  by  APL,  includes  all  the  sup- 
port instrumentation  required  to  provide  attitude  con- 
trol, data  storage  and  transmission,  command  control, 
and  power. 


DISCUSSION 

The  specifications  developed  by  NASA  and  MIT 
for  the  scientific  missions  required  a complete  redesign 
from  those  used  in  the  SA^A  and  SAS-B  satellites. 
The  structure  was  required  to  support  up  to  400  lb  of 


experiment  instrumentation,  and  the  larger  power 
system  involved  a new  solar  array  panel  configuration 
and  a new  deployment  method.  Three-segment  solar 
array  panels  with  triple-hinge  assemblies  enabled  the 
panels  to  be  folded  inside  the  heat  shield  during 
launch  and  deployment.  New  design  techniques  en- 
sured that  the  array  deployed  symmetrically,  without 
colliding  with  any  part  of  the  spacecraft.  Computer 
simulation  of  the  hinge-spring  constants  and  panel 
mass  provided  good  insight  to  its  deployment  char- 
acteristics. An  actual  deployment  test  with  an  engi- 
neering model  verified  the  computer  simulation  and 
proved  the  adequacy  of  the  design. 

Attitude  System.  The  attitude  control  system  of 
SAS-C  (Refs.  1 and  2)  provides  considerably  more 
precise  performance  than  did  those  of  SAS-A  and 
SAS-B.  The  heart  of  the  SAS-C  attitude  system  is  an 
active  closed-loop  control  of  spin  rate  with  redundant 
integrating  rate  gyro  sensors  and  ground  command 
call-ups  for  a specific  spin  rate.  The  desired  spin  rate 
is  transmitted  to  the  satellite  as  a digital  word  that  is 
converted  into  an  analog  voltage  to  be  compared  with 
the  analog  output  of  the  rate  gyro.  The  error  voltage 
is  then  used  to  control  the  speed  of  the  reaction 
momentum  wheel,  thus  adjusting  the  spin  rate  of  the 
outer  satellite  body  to  the  required  value.  Many  other 
features  are  described  in  Refs.  1 and  2. 

Telemetry  System.  The  telemetry  system  is  basic- 
ally a standard  pulse  coded  modulation  (pcm)  digital 
readout  phase  modulating  a vhf  transmitter.  The 
telemetry  system  has  two  modes  of  operation;  a fixed 
format  established  and  programmed  into  a prom 
memory  before  launch,  and  a variable  format  using 
core  storage.  The  variable  format  can  be  changed 
from  the  ground  by  loading  specific  programs  into 
the  core  storage.  In  addition,  bit  rates  can  be  changed 
from  125  bits  per  second  in  doubling  steps  to  16  000 
bits  per  second. 

By  ground  command,  operation  can  be  switched 
from  fixed  to  variable  format.  There  are  two  redun- 
dant sections  in  the  fixed  format  and  two  in  the  varia- 
ble. The  variable  format  was  operated  during  launch. 
In  this  way  a launch  format  was  available  in  one  sec- 
tion and  an  experiment  data  format  in  the  other.  After 
launch,  a single  ground  comnund  was  used  to  switch 
from  the  first  section  to  the  second,  enabling  a diange 
from  launch  housekeeping  d?ta  to  experiment  data.  An 
ultrastable  crystal-controlled  oscillator  is  used  to  gen- 
erate the  clodc  for  the  telemetry  system  as  well  as 
other  frequency-dependent  systems.  TTris  oscillator  has 
a basic  stability  of  10—*®  per  hour  within  the  normal 
operating  satellite  temperature.  Thus  all  the  timing 
and  control  functions  of  the  satellite  have  a stable  fre- 
quency source.  Other  important  features  of  the  telem- 
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etry  system  arc  three  64-channel  analog  commuta- 
tors and  two  l6-<hannel  digital  commutators,  an 
analog-to-digital  converter,  error  coding,  and  tape 
recorders  for  the  storage  of  data. 

Command  System.  The  SAS-C  command  system 
uses  the  standard  NASA  64-bit  fsk  word,  AM  modu- 
lated 50%  by  a sine-wave  bit-synchronizing  signal. 
The  system  consists  of  redundant  antennas,  receivers, 
bit  detectors,  decoder  logics,  power  switching  relay 
coil  matrices,  and  delayed  command  systems  con- 
nected so  that  either  redundant  system  can  be  chosen 
by  proper  coding  of  the  command  word  rather  than 
by  switching  system  connections.  It  provides  56  relay 
commands  on  and  off  and  any  number  of  data  com- 
mands utilizing  24  bits  of  the  command  word  for 
coding.  A total  of  15  commands  of  any  combination 
can  be  stored  in  each  redundant  delayed  command 
system,  giving  a total  of  30  delayed  commands  if 
desired.  Delays  of  from  3 s to  2.4  h can  be  achieved 
with  about  2-s  minimum  separation.  The  resolution 
time  of  execution  of  the  first  command  can  be  on  the 
order  of  0.5  ms  if  all  parameters  affecting  transmission 
time,  etc.  are  considered.  A program  of  delayed  com- 
mands can  be  stored  for  continual  recycling  since  the 
commands  are  not  destroyed  when  used  but  are  re- 
cycled into  storage.  It  is  not  necessary  to  use  all  15 
commands  in  either  of  the  two  storage  systems  since 
zero  bits  can  be  inserted  to  fill  the  complete  storage, 
thus  recycling  all  bits  to  their  original  positions. 

The  delayed  command  system  is  of  special  impor- 
tance to  the  operation  of  the  SAS-C  attitude  control 
system.  It  allows  changes  in  altitude  such  as  -h  z-axis 
precession,  spin  rate  changes,  etc.  to  be  done  at  any 
time  in  the  orbit  rather  than  only  when  the  spacecraft 
is  directly  over  the  control  station.  Other  commands 
can  also  be  executed  by  the  delayed  command  system, 
e.g.,  turning  on  heaters  and  other  housekeeping  sys- 
tems at  the  beginning  or  end  of  a sunlit  period. 

Power  System.  The  SAS-C  power  system  includes 
a solar  cell  array,  a nickel-cadmium  battery,  and  a 
charge  control  system.  The  solar  cell  array  provides  a 
minimum  of  65  W (orbit  average),  and  the  battery 
has  a nominal  capacity  of  9 Ah. 

The  charge  control  system  consists  of  a digital 
coulometer  including  a voltage  limiter.  The  coulometer 
meters  the  ampere  minutes  charge  and  discharge.  The 
coulometer  is  adjusted  so  that  when  101,  105,  110,  or 
125%  (as  chosen  by  ground  control)  of  the  discharge 
power  has  been  returned  during  the  charge  cycle,  the 
unwanted  current  is  shunted  by  the  coulometer  and 
voltage  limiter. 

The  charge  control  circuits  are  redundant  and  can 
be  chosen  by  command  switching.  A low-voltage 


sensing  switch  removes  the  battery  and  all  loads  except 
the  command  system  from  the  bus  when  the  battery 
teaches  a predetermined  value  of  13.2  V.  The  battery 
is  still  connected  to  a small  portion  of  the  solar  cell 
array,  allowing  it  to  charge  at  a trickle  charge  rate. 
The  nominal  battery  voltage  is  16.1  V. 

The  solar  cell  array  blades  can  be  rotated  by  ground 
command  to  collect  the  sun’s  energy  efficiently  as  the 
spacecraft  attitude  changes  with  respect  to  the  sun.  In 
this  way,  a relatively  constant  array  power  can  be  made 
available  at  all  attitudes. 

Thermal  Control.  An  active  thermal  control  sys- 
tem is  used  in  SAS-C  together  with  some  passive 
thermal  insulation.  Thermal  louvers  were  installed 
between  the  heat  dissipating  elements  such  as  the  bat- 
tery and  electronics  books  and  the  thermal  radiator. 
The  louvers,  similar  to  Venetian  blinds,  open  and 
close,  depending  on  the  temperature  of  the  thermal 
sensors.  When  the  louvers  are  open,  heat  is  released 
through  a low  radiation  path  to  the  radiator.  When 
the  louvers  are  closed,  the  radiation  path  resistance  is 
increased  and  heat  is  conserved  by  restriction  to  the 
radiator  paneb.  Thus  a relatively  constant  internal 
temperature  is  maintained  with  changing  heat  loads, 
either  internal  or  external. 

PosTLAUNCH  OPERATIONS.  The  Spacecraft  was 
shipped  to  tho  San  Marco  Equatorial  Range  on  5 
April  1975  and  was  launched  as  scheduled  on  7 May 
1975.  The  first  few  days  of  orbit  operation  indicated 
that  all  systems  were  functioning  as  designed.  Scien- 
tific observations  began  about  four  or  five  days  after 
launch.  It  was  then  discovered  that  the  passive  nuta- 
tion damper  had  stuck  to  one  side  and  no  longer  was 
able  to  oscillate  with  nutation  motion.  After  a week  or 
so,  the  nutation  damper  became  free  and  has  been 
operating  satisfactorily  since. 

The  active  spin  rate  control  has  been  performing 
beyond  expiectations.  A dither  mode  has  been  used 
wherein  the  spin  rate  is  made  to  vary  in  an  oscillatory 
fashion  to  permit  scanning  of  a specific  X-ray  source 
until  it  is  obscured  by  the  earth.  The  dither  mode  can 
be  repeated  as  necd^  by  allowing  the  delayed  com- 
mand program  to  be  recycled. 

A zero  spin  rate  command  is  initiated  at  a specific 
time  by  means  of  the  delayed  command  feature  to 
provide  a constant  space  stabilized  azimuth  angle  for 
a period  of  time.  This  gives  a longer  integration  time 
for  the  data  from  a simple  source  than  does  the  scan- 
ning mode. 


A variable  spin  rate  using  a faster  or  slower  spin 
rate  during  the  dark  side  with  respect  to  the  sunlit  side 
can  be  programmed  during  one  orbit  so  that  the  inte- 
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grated  spin  rate  is  one  revolution  per  orbit.  This 
enables  the  X-ray  sensors  to  scan  the  same  portion  of 
the  sky  during  each  orbit  for  a long  integration  time. 


The  orientation  of  the  spin  axis  is  controlled  or 
maneuvered  by  programming  the  strength  of  the  z 
magnetic  dipole  reacting  against  the  earth's  magnetic 
field.  A desired  declination  and  right  ascension  can 
be  accomplished  to  1 or  2°.  The  z— axis  drift  is  less 
than  0.5°  per  day. 


REFERENCES 

1.  G.  A.  Fountain  and  F.  F.  Mobley.  "SAS-C  Attitude  Control 
System,  Atphea  Pbptcs  Laboratory  Developmrnts  in  Science 
and  Technology.  Fiscal  Year  1974,  APL/JHU  DST-2. 

2.  American  Institute  af  Aeronautics  and  Astronautics  (AlAA). 
Paper  No.  74-901.  August  1974. 

“A  Programmable  Sampling  Format  Telemetry  System."  Pro- 
ceedings  of  1973  International  Telemetering  Conference.  October 
1975.  pp.  350-361. 

4.  M.  R.  Peterwn.  "A  Programmable  Sampling  Format  Telemetry 
System  for  the  Small  Astronomy  Satellile-C"  Applied  Physics 

J97|~Ap£/?HU^I?StA"  Fijca!  Y,tr 

Author:  H.  B.  Riblet 

Support:  NASA  Office  of  Space  Science 


GEOS-C  SPACECRAFT  LAUNCH  AND 
EARLY  ORBIT  OPERATIONS 


On  9 April  1975,  the  GEOS-C  spacecraft,  built  by 
APL,  was  successfully  launched  by  a NASA  Delta  Class 
1410  rocket  from  SAMTEC,  Vandenberg  Air  Force 
Base,  CA.  The  satellite  is  demonstrating  the  feasibility 
of  radar  altimetry  from  satellites  as  well  as  the  feasi- 
bility of  satellite-to-satellite  tracking. 


BACKGROUND 


APL  began  developing  the  GEOS-C  spacecraft  for 
NASA  Wallops  Flight  Center  in  July  1972.  GEOS-C 
represents  a transition  from  the  National  Geodetic 
Satellite  Program,  which  included  the  ANNA  and 
GEOS-1  and  -2  spacecraft  developed  by  APL,  to  the 
emerging  Earth  and  Ocean  Dynamics  Applications 
Program.  In  addition  to  the  radar  altimeter,  this 
spacecraft  carries  an  S-band  transpxmder,  two  C-band 
transponders,  a laser  retroreflector,  and  a dual- 
frequency doppler  beacon.  It  is  being  tracked  by 
ground-based  S-  and  C-band  radars,  by  laser  trackers 
from  both  fixed  and  mobile  sites,  by  the  U.S.  Defense 
Mapping  Agency’s  network  of  satellite  tracking  sta- 
tions (TRANET)  including  Station  1 1 1 at  APL,  and 
by  geoceivers.  This  multiple  tradcing  capability  serves 
to  calibrate  the  radar  altimeter  accurately  and  is  pro- 
viding further  calibrations  of  NASA  and  other 
ground-based  C-band  systems  as  well  as  data  for 


continued  geodetic  studies.  The  laser  retroreflector 
also  supports  long-term  (20-year)  geodetic  studies 
as  well  as  the  San  Andreas  Fault  Experiment. 


DISCUSSION 

Following  liftoff  at  2358  UT  on  9 April  1975 
(Fig.  1),  near-nominal  performance  of  the  Delta 
14 10  launch  vehicle  placed  GEOS-C  into  the  follow- 


ing  orbit: 

Epoch 

75Y10M07D,  00  h 00.00  min  UT 

Apogee 

853.05  km 

Perigee 

829.56  km 

Period 

101.74  min 

Inclination 

114.989° 

Eccentricity 

0.001627 

Semimajor 

axis 

7219.47  km 

These  parameters  are  all  well  within  the  tolerances 
specified  before  launch  and  will  facilitate  carrying 
out  the  scientific  objectives  of  the  mission.  The  orbit- 
ing satellite  has  been  identified  by  NASA  as  GEOS-3 
(international  designation  1975  027A,  Spadat  Cata- 
log No.  7734). 

Prior  to  launch,  GEOS-C  was  given  a comprehen- 
sive series  of  electrical  performance  and  environ- 
mental tests  at  APL  and  at  NASA  Goddard  Space 
Flight  Center.  These  tests  extended  from  late  in  the 
spring  of  1974  through  February  1975.  The  flight 
radar  altimeter  was  received  from  General  Electric 
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Fig.  1 Liftoff  of  Delta  rocket.  Serial  No.  109,  on  9 
April  1975,  which  successfully  launched  the  GEOS-C 
spacecraft  into  the  desired  orbit. 


Co.  in  September  1974  and  integrated  with  the 
spacecraft.  The  test  series  verified  the  ability  of  the 
spacecraft  and  its  experiments  to  meet  all  perform- 
ance retjuirements  and  showed  that  performance 
would  not  be  adversely  affected  by  the  launch  and 
orbital  environments. 

Only  two  significant  anomalies  occurred  during  the 
test  program.  Thermal  balance  tests  determined  that 
the  internal  spacecraft  temperatures  were  too  low. 
Localized  changes  in  thermal  coatings  and  adjust- 
ments to  thermal  insulation  corrected  the  problem. 
During  the  cold-temperature  excursion  portion  of  the 
thermal-vacuum  performance  test,  the  telemetry  trans- 
mitter modulation  index  decreased.  The  transmitter 
was  removed  from  the  spacecraft  and  subjected  to  a 
series  of  special  tests,  which  indicated  that  the  anom- 
aly was  probably  a tuning  problem  within  the  trans- 
mitter induced  by  low  temperatures  within  a narrow 
range.  Comparison  with  the  acceptance  test  results 
for  the  transmitter  indicated  that  the  condition  may 
have  existed  at  that  time  and  was  probably  stable. 

Based  upon  these  conclusions  and  calculations 
28  showing  that  some  link  margin  exists  in  the  presence 


Fig.  2 GEOS-C  mounted  on  Delta  launch  vehicle, 
showing  half  of  the  nose  fairing  installed. 


of  the  anomalous  behavior,  and  the  fact  that  even 
if  the  •transmitter  were  to  fail  altogether  the  Pro- 
ject Approval  Document  objectives  could  be  achieved 
by  using  the  S-band  link,  it  was  decided  to  ffy  the 
transmitter  unchanged.  APL  then  completed  the 
spacecraft  acceptance  tests  and,  following  a NASA 
Consent  to  Slip  Review  on  15  and  14  February  1975, 
shipped  the  spacecraft  to  Vandenberg  Air  Force 
Base.  On  1 April,  the  spacecraft  was  mounted  on 
the  launch  vehicle  (Fig.  2),  and  final  tests  were  con- 
ducted during  the  countdown. 

On  the  first  pass  after  injection  into  orbit  telem- 
etry telltales  revealed  that  all  systems  were  still 
in  the  launch  state  and  the  solar  array  was  producing 
electrical  power.  A test  load  was  injected  into  the 
delayed  command  system  from  the  NASA  Winkfield, 
England,  station  during  this  first  orbit,  and  actuation 
of  the  commands  was  observed  at  the  NASA  Orroral, 
Australia,  station. 
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During  the  second  orbit,  the  gravity  gradient  boom 
was  deployed  approximately  0.76  m (2.5  ft),  and  the 
planned  gravity  capture  scheme  was  initiat^.  On  1 1 
April,  the  boom  was  extended  to  about  1.2  m (4  ft) 
to  excite  some  yaw  motion.  By  12  April,  the  gravity 
capture  sequence  was  completed,  with  the  boom 
fully  extended  and  then  trimmed  to  about  6.5  m (21 
ft  2.45  in.).  On  14  April,  the  pitch  angles  were  less 
than  3.7°,  and  roll  and  yaw  angles  were  less  than  1°. 
Specified  attitude  control  (less  than  1°  for  all  axes) 
was  achieved  during  the  week  of  20  April. 

During  the  gravity  capture  phase,  the  C-band 
transducers  were  activated  and  the  doppler  beacons 
were  turned  on  and  successfully  tracked  by  TRANET. 
Nominal  ope»-ation  of  all  spacecraft  systems  was  veri- 
fied, and  satellite  temperatures  were  found  to  be  with- 
in the  expected  range. 

Successful  operation  of  the  radar  altimeter  in  global 
mode  was  achieved  during  orbits  66,  69,  and  75  (l4 
and  15  April).  The  S-band  antenna  switching  sys- 
tem checked  out  correctly  on  15  April.  During  orbit 
89  (16  April),  the  S-band  satellite-to-satellite  link 
was  verified.  For  this  demonstration  the  satellite  was 
in  an  all-up  condition,  with  the  altimeter  functioning 
and  its  data  being  relayed  via  the  Applications  Tech- 
nology Satellite  6 (ATS-6)  to  the  NASA  Rosman, 
NC,  station.  Altimeter  data  were  also  telemetered  via 
the  VHF  link  directly  to  NASA  stations  at  Quito, 
Ecuador,  and  Merritt  Island,  FL. 

Engineering  evaluation  of  the  satellite  was  com- 
plete on  20  April;  experiment  calibration  commenced 
the  following  day  and  ontinued  until  20  May.  Dur- 
ing this  period.  Wallops  ( inducted  underflights  with 
C-54  aircraft  to  provide  ground  truth. 


The  satellite  continued  to  perform  satisfactorily  un- 
til an  evening  pass  on  12  May  1975,  when  the  com- 
mand system  started  throwing  spurious  relays  when- 
ever OFF  relay  commands  were  sent  through  logic  1 
matrix  1.  The  ON  relay  commands  through  matrix  1 
and  ON  or  OFF  commands  through  matrix  2 are  oper- 
ating properly.  Analysis  of  the  problem  indicated  the 
most  probable  cause  to  be  a wire  shorted  across  a 
diode.  Recommendations  on  how  to  reduce  the  opera- 
tional consequences  of  the  failure  to  a negligible  level 
were  given  to  NASA  and  have  been  implemented. 

Since  20  May  1975,  the  satellite  has  been  used  for 
global  activities  (Phase  C of  the  mission  operations 
plan)  wherein  the  ATS-6  satellite  is  used  to  obtain 
four  arcs  of  Satellite-to-Satellite  Experiment  (SSE) 
data  per  day.  During  each  of  the  SSE  atcs,  radar 
altimeter  data  will  be  scheduled  for  most  portions  of 
the  arcs  that  are  over  ocean  areas.  Also,  C-band, 
S-band,  laser,  and  doppler  data  are  being  obtained 
to  support  orbit  determination  for  the  SSE  and  radar 
altimeter  data  arcs. 

During  Fiscal  Year  1976,  the  GEOS-3  satellite 
will  enter  Phase  D of  the  mission  operations  plan, 
which  consists  of  unique  experiments  and  localized 
grid  activities.  In  addition  to  these  major  activities, 
intermittent  calibration  and  evaluation  activities  will 
be  conducted  at  the  nominal  rate  of  once  per  month. 


Author:  L.  D.  Eckard,  Jr. 

Support:  NASA  Wallops  Flight  Center 
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HYDRCX3EN  MASER  FREQUENCY 
STANDARD  IN  OPERATION  AT  APL  TIME 
AND  FREQUENCY  STANDARDS 
LABORATORY 


The  APL  Time  and  Frequency  Standards  Laboratory 
has  placed  in  operation  a NASA  Model  NP  hydrogen 
maser  frequency  standard.  The  unit,  the  primary  ref- 
erence standard  for  a NASA  hydrogen  maser  develop- 
ment program  at  APL,  permits  the  Laboratory  to  make 
more  accurate  stability  measurements  than  were  pre- 
viously possible.  The  hydrogen  maser  is  also  being 
used  in  the  evaluation  of  precision  time  and  frequency 
instrumentation  developed  for  other  APL  programs. 


BACKGROUND 

Precision  time  and  frequency  measurements  depend 
strongly  on  the  quality  and  stability  of  the  reference 
standards  available..  When  the  reference  standards 
are  of  a quality  comparable  to  the  devices  being 
measured,  tedious  statistical  means  must  be  used  to 
separate  the  errors  associated  with  the  device  from 
those  of  the  reference.  However,  when  the  reference 
standard  is  highly  superior,  the  measurements  can  be 
made  directly  and  less  data-processing  time  is  needed. 

The  atomic  hydrogen  maser  is  the  most  stable  type 
of  frequency  source  known.  Its  short-term  stability, 
or  spectral  purity,  is  outstanding,  making  it  a supe- 
rior reference  standard. 


DISCUSSION 

APL  received  the  NASA  hydrogen  maser  Model 
NP-'  (Fig.  1)  from  Goddard  Space  Flight  Center 
in  early  1975  for  retention  as  a reference  working 
model.  APL,  with  Goddard  assistance,  has  restored 
the  maser  to  working  condition.  The  APL  Time  and 
Frequency  Standards  Laboratory  has  gained  experi- 
ence in  its  operation  and  monitors  it  continuously. 

The  hydrogen  maser  is  an  active  oscillator  powered 
by  energy  from  electron-proton  magnetic  interaction 
of  excited  hydrogen  atoms.  The  working  components 
of  the  hydrogen  maser  are  shown  in  Fig.  2.  Mole- 
cular hydrogen  gas  is  introduced  into  a dissociation 
bulb.  TTiere,  rf  excitation  converts  the  molecular  hy- 
drogen, to  excited  atomic  hydrogen,  which  fjasses 
through  a collimator  (a  number  of  small  parallel 
tubes)  into  a state  selection  magnet.  Atoms  of  the 
desired  atomic  energy  are  focused  into  a beam  that 
passes  through  a hole  in  the  vacuum  enclosure  and 
into  a Teflon-lined  quartz  storage  bulb  that  is  centered 


Fig.  1 NASA  Model  NP-1  hydrogen  maser. 


in  a resonant  microwave  cavity.  The  atoms  moving 
with  thermal  kinetic  energy  collide  with  the  wall  of 
the  storage  bulb  and  are  reflected  back  and  forth  until 
they  finally  leave  by  the  hole  they  entered  and  are 
absorbed  by  the  VACION  vacuum  pump.  While  the 
excited  atoms  are  in  the  storage  bulb  (about  one 
second),  they  give  up  their  excited  state  energy  to 
increase  the  energy  in  the  rf  field  o'  the  resonant 
microwave  cavity.  In  accordance  with  the  uncertainty 
principle,  this  unusually  long  energy  decay  time  (about 
one  second)  provides  the  characteristic  for  a very 
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monochromatic  radiation  field  to  exist  in  the  cavity. 
TTiis  ener^  at  1 420  405  751  Hz  is  coupled  out 
loosely  and  used  to  phase  lock  standards  at  more  con- 
venient working  fretjuencies  such  as  5 MHz. 

The  hydrogen  maser  has  been  used  to  evaluate 
the  performance  of  another  unit  of  the  same  design, 
NP-l/NP—4,  and  to  measure  the  performance  of  the 
cesium  fretjuency  standard  of  the  APL  Time  and  Fre- 
quency Standards  Laboratory.  It  has  also  been  used 
to  measure  the  characteristics  of  two  crystal  oscillators 
designed  for  spacecraft  applications.  The  results  are 
presented  as  Allan  variances  in  Fig.  3.  The  Allan 
variance  is  a statistical  determination  of  the  uncer- 
tainty of  a frequency  measurement  made  over  a speci- 
fied interval  oT  time  using  the  device  under  test  as 
a time  standard. 

As  a part  of  the  maser  development  program,  APL 
is  training  staff  members  to  make  major  repairs  on 
the  NP  type  masers  to  sustain  service  from  these 
designs  in  experimental  field  operations.  The  NP— 4 
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Fig.  3 Allan  variance  measurements  of  precision  fre- 
quency standards  relative  to  a hydrogen  maser. 


maser,  the  second  unit  to  be  worked  on,  was  dis- 
assembled, repaired,  revised,  and  reassembled.  After 
a training  and  education  period  of  five  months,  APL 
feels  qualified  to  build  for  NASA  additional  hydro- 
gen masers  of  an  improved  design  for  field  service. 


Author:  L.  J.  Rueger 

Support:  NASA  Goddard  Space  Flight  Center 


A “STORE  AND  FORWARD"  TELEMETRY 
COMPUTER  PROGRAM  DESIGNED 
USING  COMPUTER  SIMULATION 


The  design  of  e set  of  lntarrupt-drlv»n  software 
algorithms  has  been  accomplished  by  APL  for  NASA 
Goddard  Space  Flight  Center  fGSFC).  The  hey 
algorithm.  Digital  Storage  System  (OSS),  eutometes 
collection,  storage,  and  forwarding  of  hIgh-rate  digital 


telemetry  Information  using  e network  of  minicom- 
puters. Extensive  simulation  testing  by  computer 
permitted  assessment  of  the  Impact  on  system 
throughput  of  verious  herdwere/softwere  tradeoffs 
before  system  Implementation. 
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BACKGROUND 

APL  was  asked  by  GSFC  in  1972  to  investigate  a 
configuration  of  computer  equipment  proposed  to  aug- 
ment the  worldwide  Satellite  Tracking  and  Data  Net- 
work (STADAN).  Increased  remote-site  capabilities 
ate  required  because  the  number  of  transmitting  satel- 
lites and  their  individual  data  rates  have  increased  far 
beyond  the  information  capacity  of  the  present  system. 
Such  disparities  as  a predicted  aggregate  data  input 
rate  of  1 000  000  bits  per  second  to  remote  sites  now 
provided  with  a maximum  data  relaying  rate  of  28  500 
bits  per  second  require  a considerable  on-site  storage 
capacity  for  non-real-time  data,  while  real-time  data 
are  immediately  forwarded. 

The  computer  equipment  complex  under  consid- 
eration for  the  remote  site,  known  as  the  Digital  Data 
Processing  System  (DDPS),  is  depicted  in  Fig.  1. 
The  existing  Univac  642-B  communications  computer 
is  to  be  augmented  by  the  addition  of  three  Digital 
Equipment  Corporation  minicomputers  (PDP-11/ 
20’s),  a large-capatity  rotating  magnetic  drum,  and 
two  high-density  tape  drives.  The  computer  labeled 
"file  processor,"  which  provides  the  data  storage,  was 
the  subject  of  APL’s  in-depth  analytical  investigation. 

APL  was  requested  to  propose  a new  computer  in- 
terface for  the  existing  rotating  magnetic  drum  and  to 
define  software  strategies  for  the  file  processor  that 


would  make  feasible  the  data  storage  concept  using 
multiple  minicomputers. 

DISCUSSION 

The  technique  of  simulation  was  used  to  allow  the 
many  hardware/software  tradeoffs  to  be  explored 
quickly  and  to  provide  a prediction  of  attainable 
performance  prior  to  implementation. 

The  specific  objectives  of  the  simulation  effort  were 

(a)  to  prove  the  correctness  of  software  algorithms, 

(b)  to  explore  the  effect  on  subsystem  performance 
of  various  equipment  alternatives,  and  (c)  to  deter- 
mine performance  parameters  of  the  optimired  sys- 
tem. IBM's  Computer  System  Simulation  Language 
(Ref.  1)  was  chosen  as  the  iimulation  vehicle  since  it 
had  sufficient  generality  to  permit  detailed  modeling 
of  the  projected  minicomputer  control  program,  but 
was  of  a high  enough  level  to  require  only  the  para- 
meterized specification  ol  standard  digital  computer 
peripherals. 

The  design  philosophy  applied  to  the  DSS  control 
program  was  one  of  interrupt-responsiveness,  moti- 
vated by  the  relentless  arrival  of  telemetry  informa- 
tion in  real  time  from  sources  external  to  the  sub- 
system. The  identification  of  key  subsystem  functions 
permitted  a fixed  priority  ranking  among  software 
components  and  suggested  an  optimum  priority  place- 


1.  REAL  TIME  STORE  AND  FORWARD  BACKUP 

2.  INPUT  1 000  000  BITS  PER  SECOND 

3.  OUTPUT  28  500  BITS  PER  SECOND 

4.  13  INPUT  DATA  CHANNELS  MAXIMUM 
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Fig.  1 Digital  Data  Procetsing  Syttem. 
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Fig.  2 Digital  Storage  System  hardware/software  scheduler. 


ment  of  computer  peripherals  on  the  file  processor 
UNIBUS.  This  placement,  combined  with  both  the 
fixed  priority  ranking  of  software  functions  and  the 
PDP-11/20  computer’s  interrupt  architecture,  per- 
mitted responsive,  highly  adaptive  software  inter- 
ruptibility. 

The  general  input/output  scheduling  strategies 
chosen  for  DSS  involve  immediate  initiation  of  an 
operation  to  an  idle  peripheral,  placing  a command 
into  a list  threaded  by  priority  if  the  peripheral  is 
still  engaged  in  a previous  operation,  and  always 
choosing  the  top  element  of  the  preordered  list  upon 
receipt  of  an  operation-complete  interrupt  (Fig.  2). 
A significant  exception  is  applied  to  the  file  processor 
drum  since,  if  its  rotation  is  unaccounted  for,  im- 
mediate initiation  of  an  operation  to  a random  data 
sector  will  cause  unacceptably  long  waiting  times  for 
a correct  sector  to  arrive  under  the  read/write  heads. 
The  strategy  developed  involves  building  many 
threaded  lists,  three  for  each  drum  sector,  representing 
the  storage  functions,  telemetry  data  storage,  network 
data  forwarding,  and  on-line  temporary  tape  storage. 


The  drum  interface  schedules  the  drum  operation- 
initiation  software  via  a position  interrupt,  which 
acknowledges  the  arrival  of  sector  S,  to  which  the 
software  responds  by  initiating  a previously  selected 
operation  for  sector  5-1-1  and  then  preselects  an  op- 
eration for  sector  5 -f  2 from  among  the  three  pos- 
sible functions  for  that  sector.  Furthermore,  the  choice 
of  preselected  operations  is  determined  by  the  amount 
of  core-resident  telemetry  data  waiting  to  be  stored  on 
the  drum. 

The  temporary  clearing  of  telemetry  data  from  the 
drum  to  the  tape  drives  is  necessary  because,  at  the 
predicted  aggregate  data  input  rates  and  the  limited 
forwarding  rates,  the  drum  capacity  is  sufficient  to 
store  only  two  minutes  of  telemetry  information  with- 
out overflowing.  The  drum  space-clearing  function, 
initiated  shortly  before  overflow,  applies  a set  of 
criteria  including  forwarding  status,  message  length, 
and  intermessage  priority  to  select  the  optimum 
message  to  be  transferred  for  drum  space  recovery 
(Fig.  3).  If  the  message  selected  has  bem  completely 
forwarded,  space  may  be  regained  quickly  by  merely 


Fi(.  3 Algorithm  for  adaptive  drum  operation  selec- 
tioo. 


Fig.  4 DDPS  performance  curve. 


cated  that  only  a 3%  increase  in  space-clearing  rate 
was  achieved  by  a 25%  increase  in  tape  speed,  and 
led  to  the  incorporation  of  a multisector  read  capa- 
bility in  the  drum  hardware  interface  design  to  reduce 
the  number  of  clearing-related  interrupts  that  were 
revealed  to  be  the  limiting  factor.  Inspection  of  the 
file  processor  software  utilization  suggests  that 
UNIBUS  loading  may  limit  subsystem  performance 
even  if  a faster  processor  is  used. 

GSFC  has  requested  that  APL  implement  a proto- 
type file  processor  subsystem  to  validate  the  simula- 
tion results  and  to  provide  a test  site  for  the  APL 
drum  interface  design. 


detecting  dnim  data  index  entries.  For  still-resident 
data,  clearing  is  accomplished  at  almost  maximum 
tape-writing  speeds  by  overlapped  read  command 
generation  utilizing  two  tape  buffers.  Finally,  since 
the  computer  operator  must  assist  in  replacing  full 
tapes,  two  tape  drives  are  devoted  to  space  clearing 
to  lessen  operator  fatigue. 

The  results  of  many  simulation  runs  are  sum- 
marized in  Fig.  4.  They  demonstrate  that  the  DSS 
software  design  is  capable  of  meeting  a maximum 
steady-state  input  rate  of  925  000  bits  per  second, 
f^omporison  of  these  results  with  earlier  runs  indi- 
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PREDICTION  OF  RAIN 
ATTENUATION  STATISTICS 
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Attenuation  by  reflection  and  scattering  of  electro- 
magnetic signals  by  rain  Increases  as  the  frequency 
Increases.  Because  of  the  present  emphasis  on  the 
use  of  frequencies  above  10  GHz  for  satellite-earth 
comrmmicetion  systems,  knowledge  of  the  statistics 
of  attemation  caused  by  rain  becomes  Increasingly 
Important.  Techniques  to  predict  attenuation  statistics 
using  radar  reflection  data  have  been  developed  for 
frequencies  up  to  100  GHz  and  for  arbitrary  path 
elevation  angles. 


BACKGROUND 

Rain  may  seriously  attenuate  satellite-earth  com- 
munication signals.  Useful  rain  attenuation  statistics 
should  not  only  provide  an  estimate  of  the  amount  of 
time  per  year  that  the  attenuation  along  an  earth- 
satellite  path  exceeds  a given  .level  but  should  also 
establish  the  optimum  spacing  between  earth  termi- 
nals so  that  a satellite  communicating  with  one  earth 
station  with  high  path  attenuation  could  automatically 
be  switched  to  a second  station  having  significantly 
less  path  attenuation  (space  diversity).  Such  statistics 
can  be  arrived  at  directly  by  means  of  experimental 
satellites  communicating  with  earth  terminals  sep- 
arated by  known  distances  or  indirectly  by  means  of 
radiometry  cr  radar  measurements.  At  present,  space 
diversity  and  attenuation  statistics  are  still  meager,  and 
those  giving  the  dependence  on  path  elevation  angle, 
frequency,  and  climatology  are  practically  nonexis- 
tent. 

Radar  can  measure  rain  reflectivity  profiles  over  a 
large  coverage  area.  By  means  of  stored  radar  data 
in  three  dimensions,  attenuation  statistics  may  be  ob- 
tained over  various  path  elevation  angles  and  fre- 
quencies. 

During  the  summer  of  1973,  radar  reflectivity  pro- 
files of  rain  regimes  were  recorded  on  digital  tape 
at  the  Radar  Atmospheric  Research  Facility  at  Wal- 
lops Island,  VA,  by  the  Radar  Atmospheric  Physics 
Group  of  APL.  These  data  were  obtained  for  various 
rain  types  out  to  139  km  from  the  radar  facility 
(SPANDAR)  and  were  stored  on  digital  tape  in 
bins  for  which  the  radar  reflectivity  is  defined  as  a 
function  of  azimuth  and  range.  By  injection  of  this 
data  base  into  modeling  procedures,  attenuation  and 
space  diversity  statistics  have  been  obtained  for  se- 
lected frequencies  in  the  13-  to  100-GHz  interval 
and  at  path  elevation  angles  from  20  to  90°  (Refs. 
1 and  2). 


Subsequent  analysis  has  also  demonstrated  the  abil- 
ity to  predict  (a)  attenuation  statistics  at  other  path 
elevation  angles  when  given  the  statistics  at  a parti- 
cular path  angle,  and  (b)  attenuation  statistics  at  a 
third  frequency  when  given  the  statistics  at  two  fre- 
quencies. 


DISCUSSION 

A useful  descriptor  in  characterizing  rain  observed 
by  radar  is  the  reflectivity  factor  Z (in  mm*/m’),  a 
quantity  dependent  upon  the  raindrop  size  distribu- 
tion. The  radar  rain  reflectivity,  which  represents  a 
measure  of  the  received  jx)wer  scattered  from  the  rain 
at  a particular  wavelength,  may  be  shown  to  be  propor- 
tional to  Z.  A computer  program  was  develo}^  that 
can  find  profiles  of  the  measured  rain  reflectivity  fac- 
tors Z(/)  along  constructed  representative  earth- 
satellite  paths  passing  through  rain  regions,  where  I 
represents  the  distance  (in  kilometers)  along  a path 
referenced  to  the  ground  terminal  The  corresponding 
profiles  for  the  attenuation  coefficient  i(l)  (in  de- 
cibels per  kilometer)  were  modeled  from  the  relation 

*(/)  = 4Z(0]*  CD 

where  the  values  of  a and  b at  the  individual  fre- 
quencies were  calculated  using  the  raindrop  distri- 
bution measured  by  Joss  et  al.  in  Locarno,  Switzer- 
land, for  thunderstorm  activity  (Ref.  3). 

The  total  attenuation.  A,  expressed  in  decibels  (also 
referred  to  as  the  "fade  depth”)  associated  with  pro- 
pagation along  a representative  path  is  given  by 

A ~ f ' k dl  (2) 

0 

where  dl  is  the  differential  length  along  the  path  and 
/„  is  the  length  along  the  path  beyond  which  the 
measured  reflectivity  is  negligible. 

Probability  estimates  have  been  obtained  by  sam- 
pling the  calculated  attenuations  along  13  080  con- 
stmcted  paths  contained  in  232  azimuths  (great  circle 
planes)  covering  a period  of  14  rain  days.  Only  those 
azimuths  were  selected  for  which  an  attenuation  of  2 
dB  or  greater  at  / 30  GHz  and  B ~ 45®  (equiva- 

lent to  approximately  8 dB  at  100  GHz)  exist^  along 
at  least  one  path  per  azimuth.  Therefore  the  proba- 
bility estimates  obtained  are  "conditional"  in  the  sense  35 
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Fig.  1 Joint  conditional  probability  that  the  anenua- 
tion  at  two  terminals  separated  by  a distance  d exceeds 
the  absdssa  at  a path  elevation  angle  t of  45°  and  at 
/ of  100  GHz. 


tfuit  an  "attenuation  event"  has  occurred.  Figure  1 
shows,  as  an  example,  probability  estimates  at 
/ - 100  GHz  and  - 45".  The  curve  labeled 

— 0 km  represents  the  probability  (single  terminal) 
that  the  attenuation  along  an  earth-satellite  path  ex- 
ceeds the  abscissa.  The  curves  labeled  </  ■=  2, 4, . . .20 
km  represent  the  joint  probabilities  that  two  stations 
separate  a distance,  d,  experience  simultaneous  fades 
that  jointly  exceed  the  abscissa.  The  curve  labeled 
optimum  spacing  represents  the  joint  probability  for 
the  case  that  two  stations  in  "independent”  rain  en- 
vironments experience  fades  that  jointly  exceed  the 
abscissa.  This  joint  probability  is  given  by  the  square 
of  the  probability  levels  of  the  single  terminal.  The 
indicated  points  in  the  curves  of  Fig.  1 represent  the 
probability  levels  calculated  using  the  experimentally 
acquired  rain  reflectivity  data  base.  All  the  curves 
are  relatively  smooth  and  approach  the  optimum 
spacing  curve  as  the  separation  increases.  Similar  prob- 
ability curves  have  bew  obtained  at  13,  18,  25,  and 
30  GHz  at  path  angles  of  30,  45,  and  60®. 

Given  the  fade  depth  Ag{e„)  for  a path  angle 
the  predicted  fade  depth  at  another  angle  9 is 
given  by 

O) 

where  we  associate  with  both  fades  equal  probabilities 
of  being  exceeded.  In  Fig.  2 we  examine  the  validity 
of  Eq.  (3)  over  a path  angle  interval  ranging  from 
20  to  90®  at  / - 100  GHz.  The  solid  curves  reprc- 
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Fig.  2 Comparison  of  fade  statistics  calculated  at  a 
series  of  path  elevation  angles  using  the  esc  $ variation 
at  / - 100  GHz. 


sent  the  single-terminal  probabilities  for  the  various 
path  angles  calculated  independently  using  Eq.  (2). 
The  circled  points  represent  the  predicted  probabilities 
using  Eq.  (3)  given  only  the  values  at  6o  = 45®.  From 
the  table  in  Fig.  2,  the  rms  errors  over  the  probability 
range  considered  for  all  the  cases  are  less  than  10%. 
Comparable  accuracies  have  been  demonstrated  for 
13  to  30  GHz  over  a similar  probability  interval. 

Close  agreement  between  the  fade  statistics  de- 
termined using  Eq.  (3)  and  those  independently  de- 
termined using  Eq.  (2)  implies  that,  over  the  range 
of  probabilities  considered,  the  environment  may  be 
assumed  to  possess  an  equivalent  uniform  reflectivity 
with  an  associated  probability  of  being  exceeded.  Un- 
der such  conditions  of  uniformity,  the  fade  depth  is 
proportional  to  path  length,  which  in  turn  is  propor- 
tional to  CSC  0 as  given  by  Eq.  (3). 

The  techniques  used  to  predict  statistics  at  a third 
frequency  when  given  similar  statistics  at  two  fre- 
quencies are  as  follows.  We  use  the  previous  con- 
clusion that  an  equivalent  uniform  reflectivity  may 
be  defined  in  a region  with  an  associated  probability 
of  being  exceeded.  We  assign  to  the  level  of  the 
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Ai  = Z„*‘  /o  (4) 

At  = <»j  Zo  * /«  C5) 


where  /«  is  the  path  length  through  rain.  We  may  state 
er|uivalently  that  the  attenuations  at  /,  and  /,  are  /f, 


and  and  each  has  associated  with  it  the  same 
probability  P,  that  these  respective  attenuations  are 
exceeded.  Given  the  fade  statistics  at  frequencies 
/,  and  /„  Eqs.  (4)  and  (5)  represent  two  equations 
whose  unknowns  are  Z„  and  /o  at  each  P„  level. 
Hence,  solving  for  these  parameters  as  a func- 
tion of  Po,  the  fade  statistics  at  any  third  frequency 
f,  may  be  obtained  using  a formulation  similar  to 
Eqs.  (4)  or  (5).  Examples  in  Fig.  3 demonstrate  the 
validity  of  the  technique.  The  curves  represent  the 
single-terminal  fade  statistics  calculated  individually 
using  Eq.  (2)  at  the  indicated  frequencies.  The  cir- 
cled points  associated  with  the  18,  23,  and  100  GHz 
curves  represent  the  predicted  levels,  assuming  only  a 
knowledge  of  the  1 3 and  30  GHz  statistics.  Note  that 
the  predicted  levek  have  rms  errors  of  less  than  5%. 
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INTRODUCTION 


The  Johns  Hopkins  University  School  of  Medicine,  the  School  of  Hygiene 
and  Public  Health,  and  the  Applied  Physics  Laboratory  have  collaborated  in 
a program  of  biomedical  research  and  development  since  1965.  An  important 
objective  of  the  program  is  to  apply  the  "know-how”  of  physical  sciences  and 
engineering  and  the  systems  engineering  approach  acquired  by  APL  in  its 
defense  and  space  research  to  the  practi:al  solution  of  problems  in  medical 
research  and  health  care  delivery.  'ITie  success  of  this  program  has  fully  justi- 
fied the  confidence  of  its  founders  that  fruitful  collaboration  of  engineering 
and  medical  specialists  would  be  attained  through  working  side  by  side  in 
solving  real  and  important  problems. 

The  program  has  grown  to  include  nearly  all  the  clinical  departments  and 
several  of  the  basic  science  departments  of  the  medical  divisions.  More 
recently,  the  Johns  Hopkins  Hospital  has  joined  the  collaboration.  Active 
program  areas  include  Nuclear  Medicine,  Ophthalmology,  Neurosensory 
Instrumentation,  Radiography,  Cardiovascular  Systems,  Prosthetic  Systems, 
Biomedical  Engineering,  and  Clinical  Engineering.  The  application  of  state- 
of-the-art  technology  has  contributed  to  basic  biomedical  research,  clinical 
diagnosis,  and  patient  therapy  through  improved  knowledge,  instruments, 
and  techniques. 

The  results  of  biomedical  engineering  research  and  development  are  dis- 
seminated by  about  10  to  20  publications  per  year  in  the  medical  and  bio- 
logical literature,  by  presentations  at  scientific  and  medical  meetings,  and 
through  the  teaching  activities  of  the  Johns  Hopkins  Medical  Institutions. 
Representative  of  the  latter  is  the  Masters  program  in  Clinical  Engineering 
developed  in  collaboration  with  the  Department  of  Biomedical  Engineering 
of  the  Medical  School,  which  is  completing  its  second  year. 

The  accomplishments  reported  here  are  representative  of  the  wide  scope 
of  activities  in  the  program,  from  basic  research  to  improved  health  cate 
delivery  systems. 


JHH  MINI  RECORD  SYSTEM 


T/ie  Mini  Record  Demonstration  was  developed  by 
APL  in  conjunction  with  the  Johns  Hopkins  Hospital 
staff  as  a pilot  effort.  The  primary  objective  of  the 
Mini  Record  Demonstration  is  to  implement  and  evalu- 
ate the  use  of  a low-cost,  medical-summary/problem- 
list  retrieval  system  for  General  Medical  Clinic  out- 
patients. At  the  conclusion  of  the  project,  a fully 
operational  system  was  implemented.  This  system  is 
currently  in  use  at  the  General  Medical  Clinic;  its 
impact  on  health  care  delivery  is  the  subject  of  an 
extended  JHH  evaluation. 


from  inappropriate  new  prescriptions  during 
an  unsch^uled  encounter. 

4.  Assure  the  presence  of  a succinct  and  easily 
maintained  patient  medical-summary /problem- 
list  in  the  Medical  Records  of  all  General 
Medical  Clinic  patients. 

5.  Evaluate  the  impact  of  computer-supported 
medical  abstracts  oji  {Jatient  care  delivery. 
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BACKGROUND 

The  Johns  Hopkins  Hospital  performs  the  services 
of  family  physician  for  a large  portion  of  the  local 
population  as  well  as  for  a fairly  extensive  outpatient 
referral  group  from  suburban  hospitals  through  ap- 
proximately 100  outpatient  clinics  servicing  almost 
100  000  patients  each  year.  The  General  Medical 
Clinic  provides  ongoing  care  to  outpatients  with 
chronic  problems  that  do  not  fall  within  the  purview 
of  specialty  clinics  such  as  Diabetes  Management. 
However,  General  Medical  Clinic  patients  do  visit 
specialty  clinics  and,  in  off  hours,  the  Emergency 
Room. 

The  General  Medical  Clinic  treats  approximately 
7000  patients  a year;  they  constitute  as  varied  a pop- 
ulation as  any  found  at  JHH.  It  was  selected  as  the 
focus  for  this  pilot  effort  because  the  availability  of 
a computer-managed  on-line  Mini  Record  could  be 
most  effectively  tested  on  a group  of  this  size  and 
makeup. 

Specific  objectives  of  the  system  were  to: 

1.  Provide  continuity  in  medical  care  for  those 
patients  who  rely  upon  the  General  Medical 
Clinic  as  their  primary  source  of  health  care 
or  are  currently  being  treated  by  more  than 
one  clinic.  This  should  result  in  better  follow- 
up for  known  problems. 

2.  Facilitate  the  implementation  and  use  of  a tele- 
phone inquiry  service  for  General  Medical 
Clinic  patients.  This  service  would  provide 
emergency  counseling,  avoid  unnecessary  off- 
hour  visits  to  the  Emergency  Room,  and  re- 
assure patients  experiencing  anticipated  com- 
plaints or  reactions. 

3.  Provide  the  minimal  information  required  so 
that  a patient  is  not  placed  at  undue  risk  either 
from  unnecessary  diagnostic  procedures  or 


DISCUSSION 

The  Mini  Record  Demonstration  was  performed 
over  a period  of  five  months  in  three  phases : require- 
ments definition,  program  and  data  base  development, 
and  operational  demonstration.  All  systems  analysis, 
computer  programming,  and  user  training  were  per- 
formed by  APL.  JHH  provided  medical  consultation 
with  regard  to  the  use  and  contents  of  the  Mini  Rec- 
ord and  was  responsible  for  operational  support  and 
system  evaluation.  At  the  end  of  the  APL-supported 
effort,  a fully  operational  system  and  data  base  were 
transferred  to  JHH  where  it  is  now  operational  and 
undergoing  extended  evaluation. 

The  system  maintains  a Mini  Record  (Fig.  1)  for 
each  General  Medical  Clinic  patient.  The  Mini  Rec- 
ord contains  demographic  and  administrative  informa- 
tion, brief  text  descriptions  of  major  medical  prob- 
lems and  history,  identification  of  currently  prescribed 
medications,  and  an  indication  of  visits  to  other  out- 
patient clinics  since  the  last  General  Medical 
Clinic  encounter.  The  Mini  Record  is  available  for 
display  from  interactive  terminals;  a hard  copy  of  the 
most  recent  Mini  Record  is  also  maintained  in  eadi 
patient’s  medical  folder.  Access  to  the  on-line  Mini 
Record  is  by  patient  history  number  or  name.  Termi- 
nals are  located  in  the  General  Medical  Clinic,  the 
Emergency  Room,  and  other  selected  areas. 

Information  for  the  Mini  Record  is  collected  from: 
(a)  an  Encounter  Form,  which  is  filled  out  during 
each  General  Medical  Clinic  visit;  (b)  the  billing 
file,  which  records  all  visits  to  other  clinics;  (c)  dis- 
charge summaries,  which  provide  problem  lists  and 
medications  for  newly  dis^arged  patients  scheduled 
for  follow-up  at  the  General  M^ical  Clinic;  and 
(d)  corrections  and  additions  to  existing  Mini  Rec- 
ords submitted  by  physicians.  Procedures  have  been 
established  to  assure  timely  and  accurate  records. 

Figure  2 illustrates  the  flow  for  processing  and 
using  the  Mini  Record.  General  Medical  Clinic  pa- 
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Fig.  1 Mini  Record  of  hypochetical  patient. 


tients  are  grouped  so  that  each  patient  is  normally  as- 
signed to  a single  “firm’'  for  medical  care.  This  "firm” 
provides  both  inpatient  and  outpatient  care  and  also 
supports  the  telephone  service  for  its  patients.  For  all 
unscheduled  encounters,  the  Mini  Record  is  queried. 
In  many  cases,  the  patient  and  the  Mini  Record  to- 
gether should  provide  sufficient  information  to  assist  in 
treatment.  Where  additional  infornution  is  required, 
the  medical  record  is  requested. 

Evaluation  of  the  Mini  Record  system  is  being 
performed  in  two  phases.  Based  upon  the  APL  dem- 
onstration, the  initial  evaluation  verified  (a)  that  the 
system  was  an  effective  aid  to  the  operation  of  the 
General  Medical  Clinic,  (b)  that  the  system  was  a 
necessary  prerequisite  to  establishing  the  desired  tele- 
phone inquiry  service  for  General  Medical  Clinic 
patients,  and  (c)  that  an  extended  evaluation  period 
would  be  required  in  order  to  assess  the  impact  of  the 
Mini  Record  and  the  telephone  inquiiy  service  upon 
health  care  delivery. 

The  second  phase  of  evaluation,  currently  in  pro- 
gress, is  being  conducted  by  the  hospital  staff  and  is 
considering  the  treatment  of  patients,  the  effectiveness 
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of  the  telephone  inquiry  service,  and  the  net  costs  and 
benefits  of  the  Mini  Record  system.  To  assist  in  this 
evaluation,  a control  group  has  been  identified  for 
detailed  analysis.  Results  of  the  extended  evaluation 
should  be  available  in  about  12  months. 
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EVALUATION  OF  POWERED  MEDICAL  MANIPULATOR 


A powered  medical  martipulator  has  been  fabricated 
to  study  the  man-machine  control  interface  for  highly 
handicapped  spinal-  cord  injury  patients.  Results  of 
initial  clinical  testing  indicate  that  the  device  may  be 
of  value  in  assisting  such  patients  In  some  activities 
of  daily  living  and  work-related  tasks. 


BACKGROUND 

One  handicapped  person  for  whom  few  or  no 
assistive  devices  are  available  is  the  quadriplegic.  He 
has  little  or  no  use  of  his  arms  and  hands  and  must 
depend  on  other  individuals  for  all  of  his  needs.  One 
possible  solution,  that  of  fitting  a powered  orthosis  to 
one  of  the  nonfunctioning  upper  limbs,  has  been 
investigated  by  many  researchers  but  was  found  to  be 
of  little  practical  value  because  of  the  complex  inter- 
action with  the  limb.  The  alternative,  a remote  pow- 
ered manipulator,  appears  to  offer  promise  as  a flexi- 
ble machine  with  which  severely  handicapped  persons 
can  carry  out  personal  and  vocational  tasks  with  some 
independence. 

The  scope  of  the  five-year  APL/JHMI  collaborative 
effort  to  develop  upper-limb  powered  prosthetic  and 
orthotic  systems  was  extended  in  1974  to  include  aids 
for  the  spinal  cord  injury  patient.  A manipulator  and 
an  integrated  work  table  environment  was  designed  to 
provide  the  following  required  functional  capabilities: 

1.  Handling  and  reading  magazines  and  news- 
papers; 

2.  Self-feeding  with  a spoon'or  fork;  and 


3.  Operating  push-button  devices  such  as  the 
Touch-Tone  telephone,  electric  typewriter,  desk 
calculator,  and  tape  recorder  without  equip- 
ment modification. 

The  initial  experimental  mode!  configuration  is 
shown  in  Fig.  1.  It  utilizes  a multidegree-of -freedom 
mechanical  "arm”  integrated  into  a work  table  to 
provide  a wide  range  of  grasping  capability.  The 
manipulator  provides  the  following  motion  capabili- 
ties to  the  patient  in  its  stand-alone  "arm”  arrange- 
ment: 

1.  Shoulder  flexion-extension, 

2.  Shoulder  turntable. 


Fig.  I Experimeotal  model  of  powered  medical  manip- 
nlator/work  table. 


3.  Elbow  flexion-extension, 

4.  Wrist  pronation-supination,  and 

5.  Hook  opening  (can  be  locked  in  any  open  po- 
sition). 

The  manipulator  is  controlled  by  first  generating  a 
pulse  train  of  from  one  to  four  pulses  to  mode-enable 
the  desired  axis  of  motion.  This  is  followed  by  pro- 
portional control  of  an  input  signal  transducer  that 
first  unlocks  the  desired  axis  of  motion  and  then  posi- 
tions the  manipulator  in  its  desired  position.  A single 
pulse  stops  the  motion  and  locks  the  mechanism.  The 
proportional  control  may  be  provided  by  eyebrow/ 
eyeglass  control,  thumb-wheel  control,  or  depression 
of  a control  paddle  with  the  chin  or  head  stick.  Selec- 
tion of  the  appropriate  control  mode  is  based  on  the 
capabilities  of  the  individual  patient.  Initial  clinical 
evaluation  has  included  testing  in  all  three  control 
modes. 


DISCUSSION 

A powered  medical  manipulator,  utilizing  concepts 
and  component  designs  from  previously  developed 
powered  prosthetic  systems,  has  been  designed  and 
fabricated.  Some  of  the  more  important  design  con- 
straints for  the  manipulator  are  that  (a)  the  effort 
required  to  control  the  device  must  be  within  accept- 
able limits,  (b)  output  motions  and  forces  must  be 
limited,  considering  the  safety  interface  between  the 
machine  and  the  patient,  and  (c)  the  electromechani- 
cal design  must  be  reliable  and  acceptably  quiet. 

The  various  personal  and  vocational  devices  are 
located  on  the  work  table  to  minimize  manipulator 
motions  as  well  as  to  allow  maximum  flexibility  to  the 
user.  The  typewriter^'calculator/Touch-Tone  tele- 
phone pad  are  mounted  in  a thin  flat  cart  so  that  the 
manipulator  can  move  them  in  or  out  of  range  of  the 
head  stick  or  mouth  stick.  The  typewriter  may  be 
loaded  by  the  manipulator  by  means  of  special  vacuum 
fingen  permanently  attached  to  the  terminal  device. 
For  patients  with  good  head-stick  or  mouth-stick  con- 
trol, a mechanical  paddle  control  input  may  be  pre- 
ferable to  other  facial  control  techniques,  such  as  the 
eyebrow/eyeglass  control  mode.  If  the  typewriter  tray 
is  moved  a short  distance  to  the  rear  with  the  manipu- 
lator, the  mouth  stick  can  also  reach  the  Touch-Tone 
telej^one  buttons.  In  this  instance,  the  manipulator 
would  be  used  to  pick  up  the  telephone  headset  and 
hold  it  near  the  patient’s  head  for  normal  phone  con- 
versation, and  the  mouth  stick  would  be  placed  in  a 
holder  during  the  actual  conversation. 


Other  work-related  tasks  include  operating  an 
electronic  calculator,  manipulating  paper  work,  and 
handling  magazines.  A magazine  may  be  removed 
from  the  rack  and  placed  in  position  for  reading;  the 
pages  are  turned  and  finally  the  magazine  is  returned 
to-  its  storage  position  in  the  radc.  Figure  2 shows  the 
magazine  pages  being  turned  by  the  manipulator. 


Fig.  2 Manipulator  turning  magazine  pages. 


In  addition  to  the  described  work-environment  ca- 
pabilities, the  system  can  provide  some  self-feeding 
capability  to  the  highly  disabled  patient.  The  type- 
writer is  moved  bade  to  its  storage  position  by  the 
manipulator  to  clear  the  area  directly  in  front  of  the 
patient  for  the  feeding  task.  A plate  is  then  placed  in 
front  of  the  patient  with  food  in  bite-sized  portions  or 
with  other  foods  that  can  be  picked  up  by  a fork  or 
spoon. 

The  medical  manipulator  has  now  been  subjected 
to  limited  clinical  testing  by  two  patients  at  the  Mary- 
land Rehabilitation  Center.  The  first  patient,  a 19- 
year-old  male  who  used  the  device  for  about  eight 
hours,  is  afflicted  with  rheumatoid  arthritis  to  the  extent 
that  his  musculoskeletal  function  is  limited  to  a few 
degrees  of  flexion-extension  of  his  elbows  and  fingers. 
He  quickly  learned  to  master  basic  control  of  the 
manipulator  and  could  perform  all  of  the  designated 
tasks  of  the  system.  However,  since  he  was  studying 
to  be  an  accountant  and  the  machine  could  not  handle 
large  ledgers  and  other  reference  materials,  he  dis- 
continued his  participation  in  the  testing. 

'The  second  p>atient  is  a middle  aged  male  who  has 
used  the  manipulator  for  about  seven  months  and  is 
continuing  to  develop  skills  in  its  use.  This  patient 
suffered  a spinal  cord  injury  about  eight  years  ago 
and  had  no  residual  hand,  wrist,  elbow,  or  shoulder 
motion.  He  initially  used  eyebrow  motion  for  control. 
In  later  clinical  testing  he  greatly  improved  his  per- 
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formance  by  using  a head  stick  to  control  the  mani- 
pulator. Since  he  had  already  developed  skill  in  typing 
using  the  head  stick,  this  control  mode  was  easy  for 
him  to  learn.  He  is  continuing  to  test  the  device  and 
is  becoming  proficient  in  some  tasks  of  interest  to 
him. 

Results  of  the  initial  clinical  testing  have  been 
promising.  The  sponsor  has  requested  APL  to  fabri- 
cate three  additional  models.  Tliese  units  will  incor- 


porate additional  sensors  and  other  features,  to  im- 
prove the  interface  between  the  patient  and  the 
machine. 
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ELEMENTAL  PHENOMENON  OF  VISION 


The  basic  nature  of  a novel  kind  of  visual  phenome- 
non— suprafusion  transients — has  been  uncovered 
through  an  abstract  theory  o1  vision  that  mathemati- 
cally predicts  physiological  suprafusion  responses  of 
the  eye-brain  system.  Retinal  and  cortical  experiments 
confirm  the  theory  and,  through  it,  are  correlated  with 
other  neural  and  behavioral  studies.  Analysis  also 
educes  potential  clinical  utility  for  the  retinal  (electro- 
retinogram,  ERG)  suprafusion  transient. 


BACKGROUND 

A light  that  fluctuates  sufficiently  rapidly  about  a 
constant  mean  luminance  appears  to  be  of  a constant 
brightness  equal  to  that  of  the  mean  luminance.  Long 
ago  this  Talbot-Plateau  law  was  recognized  to  require 
not  only  that  the  fluctuation  be  faster  than  a flicker 
fusion  frequency  but  also  that  it  be  steady.  Neverthe- 
less, in  this  suprafusion  regime  only  passing  investiga- 
tion has  been  made  beyond  that  needed  to  establish 
the  steady-state  precision  of  the  classic  Talbot- 
Plateau  law.  Recently  a novel  technique  for  investigat- 
ing unsteady  modulation  above  fusion  was  developed, 
largely  at  APL,  combining  psychophysical  observa- 
tion, physiological  experimen*^,  and  mathematical 
theory.  A basic  phenomenon  proves  to  be  the  tran- 
sient brightness  change  that  is  perceived  when  the 
light  fluctuation  waveform  is  abruptly  changed,  as  by 
the  sudden  jump  in  waveform  period  depicted  at  the 
top  of  Fig.  1.  The  perception  of  such  "suprafusion 
transients”  was  studied  psychophysical ly  earlier,  and 
its  sensory  basis  was  formulated  into  an  abstract 
theory  (Ref.  1).  The  theory  predicted  the  mathe- 
matical form  of  associated  physiological  responses  of 
the  visual  system  that  have  now  been  detected  in  the 
eye  and  brain  of  monkeys  (Refs.  2 and  3).  The 
analysis  rnd  interpretation  of  these  retinal  (erg)  and 
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Fig.  1 Period-juim>  stimulus  (/i  -*  ti)  and  resulting 
Green’s  functions  for  retina  (ERG)  and  cortex  (VER). 


cortical  (visually  evoked  responses,  ver)  suprafusion 
transients  are  discussed  briefly  here  (cf.  Ref.  4). 


DISCUSSION 

The  elemental  nature  of  suprafusion  transients  is 
manifest  in  the  theoretical  deduction  that  they  repre- 
sent fundamental  Green's  functions  of  the  visual  sys- 
tem. For  example,  the  waveform  period-jump  at  the 
top  of  Fig.  1 is  predicted  to  elicit  suprafusion  re- 
sponses of  the  form 


PC  AMPLITUDES 


Rp(l)  O2-  t^)•Gp(t,0)  (1) 

in  any  linearizable  part,  P,  of  the  system,  where 
<2  — is  the  period- jump  and  Gp(/,0)  is  the  Green's 
function  defined  as  the  response  of  P at  time  / to  an 
impulse  at  time  0.  Equation  (1)  was  derived,  dis- 
cussed, and  accurately  verified  by  human  behavioral 
experiments  (Ref.  1). 

In  the  physiological  studies  (Refs.  2 and  3)  the 
period-jump  stimulus  was  applied  to  monkeys  with 
electrodes  implanted  on  the  visual  striate  cortex  and 
contact  electrodes  on  the  cornea.  A principal-compo- 
nents analysis  has  extracted  from  the  responses  the 
hypothetical  retinal  (erg)  and  cortical  (ver)  Green’s 
functions  shown  in  the  lower  part  of  Fig.  1.  Aside 
from  the  amplitude  difference  and  the  expected  time 
delay  in  transmission  from  retina  to  cortex,  the  two 
waveforms  are  remarkably  similar. 

Figure  2 illustrates  the  good  accord  between  the 
experimental  data  and  the  theoretical  prediction  for 
the  cortical  responses.  Since  for  the  ver  one  may 
distinguish  two  main  parts  of  the  system — the  tradi- 
tional "on-”  and  "off-”  brightness  subsystems — we 
have  P in  Eq.  (1)  equal  to  plus  or  minus  (say), 
depending  on  whether  — /,  is  greater  or  less  than 
0.  Then  Eq.  ( 1 ) predicts  straight  lines  radiating  from 
the  origin.  This  prediction  is  satisfactorily  confirmed 
by  the  data  points  (Fig.  2).  Similar  confirmation  is 
obtained  for  the  retinal  responses. 

With  the  theory  validated,  we  take  the  erg  and 
ver  response  components  in  Fig.  1 as  Green’s  func- 
tions for  the  retina  and  cortex.  Then  their  Fourier 
transformations  must  be  frequency  response  curves. 
'The  calculated  erg  and  VER  transforms  in  Fig.  3 dis- 
play reasonable  correlation  with  other  behavioral  and 
neurophysiological  response  curves. 

Detailed  analysis  now  in  process  of  the  erg  tran- 
sients finds  a substantial  component  that  we  can  only 
associate  with  amacrine  cell  activity  in  the  retina.  If 
so,  the  ERG  suprafusion  transient  affords  increased 
anatomical  resolution  that  could  be  clinically  useful 
for  the  diagnosis-  of  macular  dysfunction. 


Fig.  2 Comparison  of  experiment  (symbols)  with 
theory  (lines)  for  cortical  responses. 
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Fig.  3 Frequeacy  response  calculated  from  Green’s 
functions  (Fig.  1,  VER  and  ERG)  and  measured  by 
flicker  perception  and  cellular  reqxMises. 
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A MONITOR  FOR  INTRACRANIAL 
PRESSURE 


A cause  of  death  in  neurosurgical  patients  and 
others  with  severe  head  trauma  is  an  undetected  rise 
in  intracranial  pressure.  For  the  protection  of  such 
patients,  an  instrument  has  been  constructed  for  clini- 
cal use  that  makes  possible  the  continuous  monitoring 
of  intracranial  pressure  with  the  least  discomfort  and 
risk  to  the  patient.  The  instrument  will  be  evaluated  in 
neurosurgical  patients  in  the  intensive  care  units  at 
several  hospitals. 


BACKGROUND 

There  are  a number  of  conditions  in  which  exten- 
sive brain  damage  or  death  may  be  caused  by  high 
levels  of  intracranial  pressure.  TTie  early  detection  of 
rising  intracranial  pressure  is  important  in  severe 
head  injuries,  in  acute  stroke  (especially  with  bleed- 
ing), in  hydrocephalic  patients,  and  in  the  postopera- 
tive care  of  neurosurgical  patients.  An  in«a..ranial 
pressure  monitor  will  give  early  warning  to  the  at- 
tending physicians  so  that  pressure-relieving  tech- 
niques are  started  before  brain  damage  can  occur. 


DISCUSSION 

TTie  intracranial  pressure  monitoring  device  de- 
signed for  this  purpose  is  an  implantable  sensor  that 
fits  into  a burr  hole  made  in  the  skull.  The  sensor  is 
nontoxic  so  that  if  disintegration  occurs  after  a long 
period  of  time,  the  degradation  products  will  not 
cause  a reaction  in  the  surrounding  tissue.  For  this 
reason  also,  the  pressure  sensor  does  not  contain 
bat  cries,  copper  wires,  or  transistor  devices. 

The  technique  makes  use  of  an  implantable  capsule 
called  a transensor,  which  contains  a radio  frequency 
resonant  circuit  (inductance  and  capacitance)  whose 
natural  frequency  is  affected  by  deformation  of  the 
capsule  caused  by  pressure  changes  in  its  environment. 
Located  outside  the  body  (Fig.  1)  but  in  close  proxi- 
mity to  the  transensor  is  an  external  detector,  or  in- 
terrogator. The  interrogator  radiates  a frequency- 
modulated  VHF  radio  signal  that  loses  energy  to  the 
transensor  when  the  signal  frequency  has  a value 
equal  to  the  transensor’s  resonant  frequency.  A bed- 
side monitor  connected  to  the  sensing  head  records 
the  frequency  at  which  the  energy  loss  occurs  and 
converts  this  measurement  to  an  intracranial  pressure 
reading. 


Fig.  1 Intracranial  pressure-monitoring  system. 


Figure  2 shows  the  steps  in  assembling  the  tran- 
sensor. In  the  second  step  (b),  the  gold-plated  nijcel 
bellows  are  glued  into  the  ceramic  half -sections  with 
epoxy.  In  the  third  step  (c),  the  ceramic  halves  arc 
glued  together,  bringing  the  closed  ends  of  the  bel- 
lows into  opposition  but  separated  by  0.1  nun.  Thus 
an  electrical  capacitance  is  produced.  As  shown  in  the 
third  step  (c),  eight  turns  of  wire  are  wrapped 
around  the  outer  diameter  of  the  ceramic  case  to  form 
the  inductance  coil.  Each  end  of  this  coil  of  wire  is 
soldered  to  one  of  the  bellows,  the  whole  now  making 
up  an  electrically  resonant  circuit.  At  this  stage,  a 
small  vent  hole  is  sealed,  trapping  a volume  of  ait 
under  the  specified  conditions  of  temperature  and 
pressure.  Finally  in  the  fourth  step  (d),  the  ceramic 
capsule  is  enclosed  in  a plastic  case  that  is  filled  with 
medical-grade  silicone  oil  and  sealed. 

Figure  3 shows  a cross  section  of  the  transensor. 
Pressure  from  a body  cavity  deflects  the  thin 
(0.13-mm)  plastic  membrane  inward.  The  pressure 
inside  the  plastic  case  is  transmitted  by  the  silicone  oil 
to  the  flexible  bellows.  The  bellows  are  pushed  closer 
together,  increasing  the  electrical  capacitance.  This 
lowers  the  resonant  frequency,  which  is  readily  de- 
tected outside  the  body  by  an  interrogator  coil  placed 
near  the  skin  over  the  implantation  site.  The  pres- 
sure on  the  transensor  is  determined  from  the  reso- 
nant frequency,  a relationship  known  from  previous 
calibrations. 

Figure  4 shows  the  transensor  and  the  mounting 
flange,  separate  and  assembled.  To  implant  the  tran- 
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Fig.  2 Assembly  of  tnnsensor. 
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Fig.  3 CroM  sectioo  of  passive  pressure  transensor. 


sensor  (Fig.  3),  the  Tiange  is  mounted  in  a burr  hole 
in  the  skull  and  secured  in  place  with  skull  screws. 
By  means  of  the  Krewdriver  shown  in  Fig.  4,  the 
transensor  is  advanced  in  the  threaded  flange  until 
its  sensing  diaphragm  barely  touches  the  dura  mater. 
The  critical  pressure  that  must  be  measured  is  the 
pressure  of  the  cerebrospinal  fluid  in  the  subarachnoid 
space. 


Fig.  4 Fltage  design  for  the  transensor.  The  transensor 
and  threaded  flange  are  shown  both  separately  and  as- 
sembled. The  screwdriver  for  advancing  the  transensor 
in  the  flange  is  shown  in  the  background. 


Currently,  intracranial  pressure  is  measured  by  plac- 
ing the  interrogator  over  the  site  of  the  transensor 
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those  taken  by  the  transensor.  In  both  patients,  the 
pressure  transducer  performed  well  and  playd  a sig- 
nificant part  in  the  management  of  the  patient's  ill- 
ness by  notifying  the  attending  physicians  of  excess- 
ive intracranial  pressure. 
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Fig.  S Method  of  transensor  implantation. 


implant  to  determine  the  transensor  rf  frequency. 
The  frequency  is  read  on  an  oscilloscope  display.  The 
frequency  value  of  the  transensor  is  converted  to  a 
pressure  reading  by  means  of  a calibration  chart.  The 
reading  must  then  be  corrected  with  a second  chart 
for  temperature  and  barometric  pressure.  This  pro- 
cedure, now  done  manually,  will  be  fully  automatic 
next  year. 

Clinical  trials  have  been  initiated  following  animal 
experimentation  and  trials.  Recently,  the  first  two 
patients  have  undergone  transensor  implants  to  moni- 
tor their  intracranial  pressures.  Occasional  spinal  taps 
to  obtain  spinal  fluid  samples  have  also  allowed  direct 
monometric  measurement  of  spinal  fluid  pressure. 
These  pressure  readings  agreed  to  within  5%  with 


FUTURE  PLANS 

The  transensors  are  now  being  made  in  quantity  for 
evaluation  at  several  hospitals.  In  addition,  a new  bed- 
side monitor  will  be  designed  for  24-hour  recording 
of  intracranial  pressure.  The  monitor  will  have  a 
built-in  barometer  to  compensate  for  barometric  pres- 
sure changes  and  will  give  an  alarm  if  the  intracranial 
pressure  reaches  a dangerous  level. 
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CALCULATIONS  AND  MEASUREMENTS 
OF  SIMULATED  ARTERIAL  FLOWS 


Calculations  of  pulsatila  flow  In  a symmetric  branch 
modeling  the  sonic  bifurcation  were  carried  out  to 
elucidate  the  role  of  hemodynamic  factors  In  athero- 
genesis.  A water  tunnel  was  constructed  to  model  the 
computational  geometry;  measured  velocity  fields 
agreed  well  with  the  calculations.  This  validation  thus 
permits  each  technique  to  be  used  to  explore  those 
aspects  of  blood  flow  to  which  it  is  best  suited,  with 
emphasis  on  fluid  mechanical  and  related  processes 
that  can  be  Involved  In  the  development  of  arterial 
disease.  These  Include  the  role  of  wall  distenslbllity, 
particle  and  solute  transport,  and  geometric  effects 
such  as  arterial  narrowing  on  fluid  mechanical  variables. 


BACKGROUND 

Much  current  theory  on  the  etiology  of  athero- 
sclerotic disease  centers  around  certain  hemodynamic 


variables.  Among  the  hemodynamic  phenomena  that 
have  been  suggested  as  being  responsible  for  athero- 
genesis  are  low  shear,  high  shear,  turbulence,  and  flow 
separation.  These  properties  of  the  flow  are  difficult  or 
impossible  to  measure  in  vivo.  Hence  our  approach  is 
via  computational  and  experimental  simulation  models 
that  approximate  the  anatomy  of  interest. 


DISCUSSION 

Interest  first  centered  on  the  aortic  bifurcation, 
which  was  simulated  by  a symmetric  Y-shaped 
geometry  (Fig.  1).  The  governing  equations  and 
method  of  solution  are  described  in  Refs.  1 and  2. 
Nonlinear  flow  solutions  for  peak  Reynolds  numbers 
(based  on  the  full  width  of  the  inlet  channel)  up  to 
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Fig.  1 The  experimental  region. 


400  have  been  computed  for  the  symmetric  90°-Y 
bifurcation  containing  an  "ideal”  (sinusoid  plus  Dc) 
pulsatile  flow.  The  solution  yields  all  of  the  variables 
of  the  flow,  including  instantaneous  streamlines,  ve- 
locity vectors,  vorticity,  vorticity  gradients,  and  wall 
shear.  It  was  shown  that  details  of  the  shear  distribu- 
tion within  the  branch  cannot  be  estimated  from 
considerations  of  mean  shear  levels  proximal  and 
distal  to  the  segment. 

Also,  transient  separation  lasting  about  20%  of  the 
eyrie  was  computed  for  a short  distance  along  the 
outer  wall  just  distal  to  the  corner.  This  suggests  that 
the  replenishment  of  blood  at  the  intimal  surface  will 
proceed  slowly  in  this  region,  to  the  disadvantage  of 
the  lining  endothelium.  If  low  wall  shear  is  a part  of 
the  etiology  of  atherosclerosis,  pulsatility  exacerbates 
the  situation  at  a branch.  Details  of  these  results  can 
be  found  in  Ref.  3. 

A more  realistic  periodic  flux  waveform,  recorded 
by  Dr.  L.  J.  Krovetz  in  a boy's  aortic  bifurcation,  was 
used  next.  Typically,  the  ratio  of  the  peak  in  vivo 
flux  to  the  mean,  R = 5,  is  greater  than  in  the  ideal 
pulsatile  wave  for  which  R = i.  For  the  larger  value 
of  R,  the  wall  shearing  stress  reaciu.:.  higher  peaks, 
and  there  is  an  indication  of  more  intense  recircula- 
tion in  the  transient  separated  region  along  the  outer 
wall  of  the  branch. 

Turning  our  attention  to  approximately  orthogonal 
arterial  junctions  (e.g.,  renals  and  intercOstals),  we 
have  also  computed  flow  within  two-dimensional 
orthogonal  branches  for  an  ideal  pulsatile  inflow.  In 
contrast  to  our  calculations  in  the  Y-branch,  the  flow 
pattern  in  such  T-branches  is  asymmetric  and  depends 
upon  the  flow  partition  However  preliminary  com- 
parisons indicate  certain  qualitative  similarities  within 
the  two  types  of  junction. 

These  numerical  studies  havt  demonstrated  that 
pulsatility  in  general  and  even  the  detailed  shape  of 


the  flux  wave  have  important  effects  on  the  stresses 
experienced  by  the  inner  (intimal)  surface  of  the 
arterial  wall.  'The  results  show  further  that  there  is 
a correlation  between  the  sites  of  hemodynamic  ex- 
trema and  of  atheromatous  lesions  in  branches. 

A water  tunnel  capable  of  producing  well-defined 
flowfields  that  can  validate  the  results  of  numerical 
computations  has  been  designed  and  constructed.  The 
tunnel  can  be  used  to  generate  steady  or  pulsatile 
flows  at  physiological  Reynolds  numbers.  A 90°-Y 
model  bifurcation,  matching  in  plan  geometry  that 
used  in  the  hemodynamic  calculations,  was  designed 
so  that  the  flow  channels  could  be  viewed  vertically 
through  glass  windows  by  a laser  doppler  velocimeter 
(LDV)  and  so  that  surveys  completely  across  the 
channels  could  be  made  without  interference.  The 
flow  cross  section  was  rectangular  with  a nominal 
inlet  aspect  ratio  of  10  (see  Fig.  2).  By  careful  de- 
sign of  the  test  section,  the  flow  velocity  can  be 
measured  at  least  as  close  to  the  wall  as  the  first  com- 
putational mesh  point  (about  400  fjjn).  The  optical 
system  of  the  LDV  has  been  improved  to  the  point 
that  the  spatial  resolution  in  the  lateral  direction  (the 
"important”  coordinate,  measured  normal  to  the  ver- 
tical walls  of  the  channel)  is  much  less  than  100  ^m. 
'This  resolution  capability  has  been  maintained  in 
operation  by  mounting  the  entire  LDV  rigidly  to  a 
translation  table  beneath  the  water  tunnel,  so  that  the 
position  of  the  beam  intersection  within  the  flowing 
fluid  is  known  accurately. 

To  test  the  validity  of  the  computational  results,  a 
pulsatile  input  was  applied  to  the  water  tunnel.  The 


Fig.  2 The  test  bifurcation. 
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EXPERIMENTAL  PROFILE 


Fig.  3 Timenlependent  expenmentAl  and  computed 
velocity  pro&lc*  at  ■ point  in  the  branch,  normalized  by 
average  upatrcam  centerline  velocity.  Note  that  the  ex- 
perimental area  ratio  is  subtly  less  than  the  value  usM 
in  the  numerical  calculadons  and  is  sufficient  to  explain 
the  minor  amplitude  difference. 


time-dependent  velocity  upstream  of  the  branch  at  the 
line  of  symmetry  was  measured.  By  Fourier  decompo- 
sition techniques,  using  solutions  for  plane  parallel 
flow,  this  velocity  was  used  to  calculate  the  inlet  flux 
wave  for  the  numerical  computation. 

Measurements  were  made  in  the  bifurcation  of  the 
velocity  component  parallel  to  the  walls.  Particular 
emphasis  has  been  placed  on  the  velocity  very  near 
the  wall  in  both  parent  and  daughter  in  the  region 
of  the  outside  corner.  Agreement  between  experi- 
mental and  computational  values  is  good  (Fig.  3), 
indicating  that  we  may  have  confidence  in  the  nu- 
merical analysis  as  well  as  the  experimental  tech- 


niques. Since  velocity  profiles  can  be  obtained  quite 
close  to  the  vessel  wall,  these  techniques  can  also  be 
used  to  estimate  wall  shear. 

Since  the  water  tunnel  is  three-dimensional  and  the 
computation  is  two-dimensional,  the  effect  of  the  inlet 
aspect  ratio  on  the  relationship  between  centerline 
velocity  and  flux  was  examined.  Fully  developed  un- 
steady parallel  viscous  flow  (viscosity  coefficient,  v) 
in  straight  rectangular  ducts  admits  an  exact  flow 
solution  in  terms  of  separable  Fourier  components. 
Each  harmonic  component  e*"*  introduces  a charac- 
teristic "viscous  length,"  V v/a>  , that — along  with 
the  dimensions  of  the  duct — determines  the  detailed 
periodic  flowfield.  When  this  analysis  was  applied  to 
the  experimental  situation,  the  aspect  ratio  in  the  tun- 
nel was  large  enough  to  permit  the  estimate  of  flux 
from  a plane  parallel  solution. 
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With  the  encouragement  of  the  Navy,  APL  conducts  applied  research  in 
civil  program  areas  where  techniques  and  expertise  developed  in  its  defense 
activities  can  make  an  especially  favorable  contribution.  In  addition  to  the 
biomedical  engineering  and  space  science  programs  discussed  in  other  sec- 
tions, these  may  be  categorized  as  transportation  systems,  environmental 
engineering,  and,  most  recently,  energy  systems. 

APL  expertise  in  automation  of  radar  signal  data  processing,  interpreta- 
tion, and  presentation  was  applied  to  a prototype  Harbor  Vessel  Traffic  Sys- 
tem developed  and  demonstrated  for  the  U.S.  Coast  Guard,  and  is  being 
utilized  in  continuing  air  traffic  safety  and  control  research  for  the  Federal 
Aviation  Administration.  Navigation  techniques  and  equipment  developed 
for  the  military  have  also  been  applied  to  civil  navigation  problems. 

APL  has  used  its  expertise  in  systems  engineering  and  controls  technology 
to  assist  the  U.S.  Department  of  Transportation  in  establishing  requirements 
and  constraints  of  urban  mass  transit  systems,  evaluating  competitive  system 
concepts  and  prototype  tests,  and  developing  control  systems  and  techniques 
for  the  operation  of  large  fleets  of  automated  vehicles.  These  studies  have 
been  extended  to  other  areas  of  transportation  planning  and  implementation 
including  personal  vehicles  and  highway  safety. 

APL  environmental  engineering  activities  fall  in  two  general  areas;  fire 
research  and  the  evaluation  of  the  environmental  impact  of  power  plant 
siting.  The  APL  fire  research  program,  initiated  by  the  National  Science 
Foundation  and  now  supported  by  the  National  Fire  Prevention  and  Control  i 

Administration,  Department  of  Commerce,  stemmed  from  expertise  acquired 
in  the  research  and  development  of  ramjet  and  rocket  propulsion  systems  for 
missiles.  It  includes  basic  and  applied  combustion  research,  fire  prevention, 
fire  fighting,  health  and  safety,  and  public  education  aspects. 

Since  1972,  APL  has  provided  to  the  State  of  Maryland  evaluation  of  the 
environmental  impact  of  specific  sites  proposed  for  nuclear  or  fossil  fuel 
electric  power  plants.  This  work,  carried  out  jointly  by  APL,  The  Johns 
Hopkins  University,  and  the  Chesapeake  Bay  Institute,  includes  evaluation 
of  air  and  water  pollution  problems,  hydrology,  biological  impacts,  waste 
product  disposal,  radioactivity  problems,  noise,  and  other  community  impacts. 

Also  included  is  assessment  of  new  technologies  for  control  and  reduction  of 
environmental  impacts.  Findings  are  made  available  to  all  interested  parties 
(the  state,  the  utility,  and  the  public)  and  have  been  pronounced  helpful  by 
all  concerned.  Studies  have  now  been  completed  on  oil,  coal,  and  nuclear 
fueled  plants.  As  a result  of  this  work,  Maryland  is  believed  to  have  one  of 
the  most  advanced  power  plant  siting  programs  in  the  country. 

Most  recently  APL  has  begun  study  of  alternative  long-range,  nondeplet- 
ing energy  sources  and  more  efficient  systems  to  meet  future  national  require- 
ments. Exploitations  of  the  thermal  gradient  between  the  surface  and  the 
deep  waters  of  the  ocean  (a  form  of  solar  energy  storage)  and  of  geothermal 
energy  sources  have  been  identified  as  important  areas  to  which  APL  can 
contribute.  In  energy  use  conservation,  more  efficient  building  heating  and 
cooling  and  personal  transportation  are  considered  areas  of  opportunity,  and 
programs  have  recently  been  initiated. 

The  articles  in  this  section  provide  a sample  of  APL’s  broad  and  increas- 
ing interests  in  civil  programs. 
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MOVING  WIRE  TECHNIQUE  STUDIES  OF 
ABLATION,  IGNITION,  AND  EXTINCTION 
OF  POLYMER  FLAMES 
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An  apparatus  and  a group  of  techniques  for  study- 
ing polymer  ffemmebility  have  been  deveioped  at  APL. 
The  method  yields  scientific  information  and  may  also 
be  of  practical  utility  as  a flammability  test 


BACKGROUND 

The  combustion  behavior  of  polymers  during  igni- 
tion and  extinction  is  an  important  practical  problem 
since  many  Ares  are  initiated  and/or  fueled  by  natural 
or  synthetic  polymers.  The  complexity  of  solid-fuel/ 
gaseous-oxidizer  systems  and  the  transient  nature  of 
ignition  and  extinction  have  made  quantitative  studies 
difficult.  Systems  are  commonly  studied  transiently,  or 
critical  time-exposure  parameters  are  determined. 

This  work  describes  a new  method  called  the  mov- 
ing wire  technique  for  studying  these  transitions.  The 
basic  concept  is  the  exchange  of  transient  time  de- 
pendence for  steady-state  position  dependence.  This  is 
accomplished  by  moving  the  material  to  be  ignited 
relative  to  the  laboratory  coordinates.  A polymer- 
covered  wire  sample  is  drawn  across  the  igniting 
flame  (see  Fig.  1 insert).  Positions  across  the  flame 
in  the  direction  of  wire  travel  correspond  to  increas- 
ing residence  times.  The  time  at  a particular  point  in 
the  flame  can  be  controlled  by  changing  the  velocity 
of  the  wire  to  provide  flexible  control  of  ignition  and 
extinction  phenomena.  Precise  measurements  are  pos- 
sible because  the  time  available  for  experiment  is 
limited  only  by  convenience.  The  use  of  moving  sys- 
tems in  fixed  laboratory  coordinates  appears  to  be  new 
in  the  study  of  ignition  and  extinction  phenomena. 


DISCUSSION 

Flame-piloted  ignition  was  investigated  initially. 
Flames  provide  both  high  temperatures  and  non- 
thermal  excesses  of  reactive  species  such  as  atoms  and 
radicals  (Ref.  1).  The  apparatus  could  also  be 
adapted  to  alternative  ignition  sources  such  as  hot  gas 
jets,  radiant  fluxes,  or  electrical  arcs.  Since  polymers 
are  solids,  ignition  requires  (a)  volatilization  of  the 
polymer,  (b)  the  availability  of  a reactant,  normally 
an  oxidizer  such  as  O,,  from  the  air,  (c)  gas-phase 
mixing  of  volatilized  fuel  and  oxidizer  on  the  molecu- 
lar level,  and  (d)  the  local  availability  of  a sufficiently 
high  temperature  and  radical  concentration  to  initiate 
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Fig.  1 Ignition,  extinction,  and  ablation  of  some 
polymers  using  the  moving  wire  technique;  critical  resi- 
dence time  b^avior  of  two  commercial  polymer-cov- 
ered wires,  #30  gauge  copper  covered  with  Teflon  and 
with  polyvinyl  chloride.  Residence  time  at  the  ignition, 
extinction,  or  ablation  transition  is  shown  as  a function 
of  composition  of  the  ignition  gases.  Exposure  time  is 
defined  as  exposed  wire  length  divided  by  wire  velocity. 


the  chain  reaction  processes  of  combustion.  Flames 
can  supply  these  requirements.  Ignition  can  occur  in 
the  wake  of  a flame  if  sufficient  oxygen  is  available 
there,  or  it  can  occur  at  the  boundary  between  a fuel- 
rich  flame  and  an  oxidizing  atmosphere. 

These  situations  are  simulated  by  two  modifications 
of  the  APL  apparatus  (Fig.  1).  The  first  is  an  iso- 
lated flame,  while  the  second  has  an  adjacent  atmos- 
phere. The  compositions  of  both  the  flame  and  the 
atmosphere  are  controlled.  The  substrate  to  be  ignited 
is  in  the  form  of  a polymer-coated  wire  that  is  pulled 
across  the  flame.  Wire  support  is  desirable  and  neces- 
sary if  the  polymer  softens  during  ignition.  The  drive 
for  pulling  the  wire  across  the  flame  is  provided  by 
a mechanism  similar  to  that  used  in  tape  decks.  A 
variable-speed  motor  provides  velocity  control.  The 


speed  is  measured  using  an  optical  counter  mounted 
on  the  drive  shaft.  Fluctuations  measured  with  a dif- 
ferential rate  meter  are  less  than  0.2%.  With  this 
drive  and  control  of  composition  of  the  flame  and 
atmosphere,  exposed  length  of  polymer,  etc.,  the  criti- 
cal transition  times  can  be  reproduced  to  within  a few 
parts  in  a thousand. 

Isolated  Burner  Experiments.  Two  types  of 
behavior  are  observed,  ablation-extinction  transitions 
and  ignition-extinction  transitions.  Ablation-extinction 
occurs  in  oxygen -poor  flames  and  is  reversible.  When 
the  wire  speed  is  decreased,  a characteristic  velocity  is 
found  at  which  visible  attack  begins.  The  involvement 
length  increav'  monotonically  as  velocity  decreases. 
Extinctio^  occuii  at  the  same  velocity  at  which 
abiaDon  began. 

Ignition-extinction  behavior  occurs  in  oxygen-rich 
flame  atmospheres  and  shows  hysteresis.  If  wire  veloc- 
ity is  decreased,  a sharp  ignition  point  is  observed. 
Movies  show  that  ignition  begins  at  the  far  edge  of 
the  flame  and  propagates  back  along  the  wire.  The 
size  of  the  initial  polymer  flame  depends  on  the  con- 
centration of  oxygen,  and  increases  with  increasing 
concentration.  R^ucing  the  velocity  increases  the  in- 
volvement. Conversely,  increasing  wire  velocity  de- 
creases involvement  until  extinction  occurs  at  a critical 
velocity  and/or  flame  size.  The  extinguishment  veloc- 
ity is  higher  than  the  ignition  velocity  so  that  a 
bistable  velocity  region  exists  (Fig.  1).  Temperature 
measurements  in  the  gas  phase  and  on  the  polymer 
surface  have  shown  that  this  behavior  is  reproducible 
and  that  the  sharpness  and  reproducibili^  increase 
with  oxygen  concentration. 

Burner-Atmosphere  Experiments.  If  a flame  is 
bounded  by  an  atmosphere  containing  a higher  free- 
oxygen  concentration,  similar  behavior  is  observed.  In 
the  flame,  free  oxygen  is  defined  as  the  fraction  of 
oxygen  remaining  after  stoichiometric  reaction;  in  the 
atmosphere,  free  oxygen  is  the  fraction  of  oxygen  in 
the  incoming  flow.  Ignition  begins  at  the  flame/ 
atmosphere  interface  but  propagates  along  the  wire 
into  the  atmosphere.  If  the  atmosphere  has  a low 
oxygen  concentration,  the  flame  will  ablatively  attack 


the  wire  prior  to  ignition.  The  behavior  is  similar  to 
that  observed  for  ablation  in  the  isolated  burner  ex- 
periments. Below  a minimum  oxygen  concentration, 
ignition  is  not  observed.  At  higher  oxygen  concentra- 
tions, the  ablative  attack  prior  to  ignition  is  reduced 
until,  above  a certain  level,  no  flame  attack  is  observed 
prior  to  atmospheric  ignition.  If  the  flame  has  a 
higher  free-oxygen  concentration  than  the  atmosphere, 
ignition  will  occur  in  the  flame  rather  than  in  the 
atmosphere,  and  the  behavior  will  be  the  same  as  that 
of  the  isolated  burner.  Typical  ablation,  ignition,  and 
extinction  behaviors  of  polyvinyl  chloride  (PVC)  are 
shown  in  Fig.  1. 

The  moving  wire  apparatus  lends  itself  to  a variety 
of  measurements  of  the  ablation-ignition-extinction 
behavior  of  polymers.  The  processes  can  be  character- 
ized by  measurements  of  temperature  and  composition 
of  the  gas  and  polymer  as  a function  of  position.  Such 
measurements  should  be  interpretable  in  terms  of 
rates  of  elementary  processes  of  heat  transfer  and 
reaction.  This  interpretation  will  require  a quantita- 
tive model,  which  APL  is  attempting  to  develop.  Pre- 
liminary studies  suggest  that  ignition  is  governed 
principally  by  surface  temperature. 

TTie  complete  understanding  of  moving  wire  ex- 
periments will  require  extended  experimental  and 
theoretical  studies,  but  the  resulting  increase  in  under- 
standing of  ignition  and  extinction  should  prove 
worthwhile.  Such  work  should  help  provide  a better, 
more  scientific  understanding  of  flammability  tests 
and  of  realistic  polymer  fire  hazard  situations. 
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EVALUATION  OF  RADAR  PROCESSING 
SUBSYSTEM 


Th«  Radar  Procassing  Subsystem  (RPS)  of  the  Auto- 
mated Radar  Terminal  System  III  (ARTS  III)  has  been 
analyzed  to  determine  the  feasibility  of  the  proposed 
design,  optimum  operating  parameters,  and  design 
simplifications  that  would  reduce  complexity  without 
Impairing  performance.  A system  configuration  is 
recommended  that  significantly  simplifies  the  original 
design  without  sacrificing  performance. 


BACKGROUND 

ARTS  III  is  the  term  used  to  indicate  the  current 
evolutionary  configuration  of  the  Automated  Radar 
Terminal  System.  The  function  of  this  system  is  to 
detect  and  track  aircraft  within  60  nmi  of  an  air 
terminal  while  minimizing  false  contacts  caused  by 
clutter.  The  system  is  automated;  all  detection  and 
tracking  of  aircraft  is  done  by  system  equipment  and 
not  by  an  operator. 

The  radar  display  provides  an  air  traffic  controller 
with  the  positions  and  velocities  of  aircraft  within  the 
coverage  region  while  suppressing  clutter  returns. 
This  allows  the  controller  to  devote  his  full  attention 
to  ensuring  adequate  aircraft  separation  and  mainte- 
nance of  flight  rules.  Since  the  radar  environment 
varies  with  weather  conditions,  ARTS  III  is  adaptive. 
The  ability  of  the  system  to  detect  aircraft  while 
rejecting  clutter  is  a function  of  the  intensity  of  the 
radar  return  from  natural  phenomena. 

In  order  to  assist  in  aircraft  tracking,  major  aircraft 
are  equipped  with  radar  beacons.  When  the  aircraft 
is  illuminated  by  a coded  radar  pulse,  the  beacon 
transmits  a coded  message  back  to  the  radar.  Because 
of  variations  in  aircraft  attitude  relative  to  the  radar, 
it  is  possible  to  receive  a radar  return  from  an  aircraft 
with  no  beacon  return,  and  vice  versa.  In  the  ideal 
case,  both  a radar  return  from  the  body  of  the  aircraft 
and  a beacon  return  are  received  and  used  to  generate 
a trade  for  that  aircraft. 

The  RPS  is  comprised  of  the  Radar  Data  Acquisi- 
tion Subsystem  and  the  associated  radar  software  of 
the  Data  Processing  Subsystem  of  ARTS  III  (Fig.  1). 

The  RPS  receives  linear  moving  target  indicator 
(MTI)  and  logarithmic  video  signals  from  the  Air- 
port ^rveillance  Radar  (ASR)  and  processes  this 
information  to  develop  target  position  reports.  These 
radar  position  reports  are  correlated  with  target  posi- 
tion reports  provided  by  the  Beacon  Data  Acquisition 
System  (BDAS)  to  yield  a composite  target  track  file 
for  each  contact  in  the  ARTS  III  coverage  area. 


DISCUSSION 

The  baseline  RPS  configuration  is  shown  in  Fig.  2. 
There  are  two  inputs,  logarithmic  video  and  linear 
MTI  video  from  the  ASR.  The  output  of  the  RPS  is 
an  estimate  of  the  position  (range  and  azimuth)  of 
each  radar  contact. 

Comparators  Qo  and  Q,  convert  log  video  into  a 
binary  signal  by  comparing  the  input  voltage  from 
the  radar  receiver  to  a variable  threshold  based  on  the 
observed  level  of  receiver  noise.  Whenever  the 
threshold  is  exceeded,  the  output  goes  from  0 to  1 . A 
"1”  signifies  a first-threshold  crossing  and  is  called  a 
hit.  The  threshold  of  Q<,  is  set  to  a higher  value  than 
that  of  Q,  by  the  addition  of  an  offset  to  the  variable 
component.  The  variable  threshold  is  set  by  feeding 
the  output  of  Qi  into  a noise  meter,  which  measures 
the  noise  level  over  a large  portion  of  the  radar  cov- 
erage area. 

The  binary  outputs  of  the  comparators  are  fed  into 
hit  processing  blocks  that  check  for  minimum  and 
maximum  pulse  length  to  ensure  that  the  length  of 
the  first  threshold  crossing  is  characteristic  of  a return 
from  a real  contact. 

The  clutter  monitor  and  clutter  mapping  blocks 
perform  spatial  and  temporal  integration  of  the  com- 
parator outputs  on  a zone  basis.  Each  zone  is  2.8°  in 
azimuth  by  2 nmi  in  range,  encompassing  about  1000 
radar  resolution  cells.  The  clutter  mapping  system 
measures  the  intensity  of  clutter  in  each  of  the  4096 
zones  within  the  radar  coverage  area  (360°  by  60 
nmi). 

Rank  quantizers  Q,  and  Qj  convert  log  normal 
video  and  linear  MTI  video,  respectively,  into  binary 
signals  by  performing  a rank  order  test  on  the  input 
video  signal.  The  rank  order  test  compares  the  relative 
amplitudes  of  the  radar  video  in  25  adjacent  radar 
resoli’.tion  cells.  The  voltage  in  die  center  cell  is  com- 
pared to  the  voltage  in  each  of  the  24  other  cells. 
Esch  time  the  voltage  in  the  center  cell  is  greater  than 
the  voltage  in  one  of  the  other  cells,  the  rank  of  the 
center  cell  is  increased  by  1.  If  the  voltage  in  the 
center  cell  is  greater  than  the  voltage  in  each  of  the 
other  cells,  the  rank  of  the  center  cell  is  25.  If  the 
voltage  in  the  target  cell  is  greater  than  the  voltage 
in  23  of  the  24  other  cells,  the  rank  of  the  center  cell 
is  24.  The  rank  computed  for  the  center  cell  is  com- 
pared to  a rank  threshold,  generally  set  to  23  or  24. 
If  the  threshold  is  exceeded,  a threshold  crossing  is 
declared  for  the  center  cell.  This  test  is  performed 
sequentially  so  that  all  cells  from  1 to  60  nmi  are 
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compared  to  the  24  cells  nearest  in  range  once  per 
scan,  resulting  in  a constant  false  alarm  rate  at  the 
output  of  the  rank  quantizers. 

The  hit  processing  blocks  for  the  rank  quantizers 
function  in  the  same  manner  as  those  for  the  com- 
parators Qo  and  Q,.  Rank  quantizers  and  Q,  drive 
two  weather  maps  via  scan-to-scan  integration  of  the 
output  of  a sequential  observer  (up-down  counter) 
monitoring  the  outputs  of  the  quantizers.  The  weather 
maps  are  used  to  control  the  output  false  alarm  rate 
of  the  quantizers  by  biasing  the  rank  order  test  pro- 
cedure as  a function  of  an  estimate  of  the  degree  of 
azimuthal  video  correlation  present  between  adjacent 
radar  sweeps. 

The  strobe  monitor  determines  the  presence  of 
transient  electromagnetic  interference,  while  the  MTI 
limit-sensing  module  determines  if  the  MTI  video  is 
limiting  (i.e.,  video  output  is  at  the  saturation  level). 

The  output  of  only  one  of  the  three  quantizers 
(Qj,  Q2,  or  Qa)  is  fed  to  the  final  detection  block 
at  any  one  time.  The  choice  is  determined  by  the  hit 
selection  block,  on  the  basis  of  the  combined  outputs 
of  the  strobe  module,  the  MTI  limit  module,  and  the 
two  quantizer  selection  maps.  In  the  presence  of  re- 
ceiver noise,  comparator  Q,  is  selected;  in  light 
clutter,  normal  rank  quantizer  Qj  is  selected;  in 
medium  clutter,  MTI  rank  quantizer  Q,  is  selected; 
and  if  MTI  video  is  limiting,  rank  Quantizer  Q,  is 
selected. 

The  video  selected  by  the  hit  selection  block  is 
integrated  in  azimuth  by  a two-level  second  threshold. 
The  first  level,  called  a predetector,  uses  an  up-down 
counter  similar  to  the  one  used  for  weather  mapping 
to  check  for  azimuthally  correlated  groups  of  radar 
cells  containing  first-threshold  crossings.  When  a pre- 
detection occurs,  the  radar  software  examines  a row 
of  17  azimuthally  adjacent  radar  cells  centered  on  the 


location  of  the  predetection  and  counts  the  number 
of  cells  that  contain  a first-threshold  crossing.  If  8 or 
more  of  the  1 7 cells  contain  a first-threshold  crossing, 
a target  detection  is  declared  (this  is  called  an  8/17 
second  threshold).  The  radar  software  then  proceeds 
to  determine  the  range  and  azimuth  of  the  centroid  of 
the  target.  The  centroid  location  is  passed  to  another 
portion  of  the  ARTS  III  Data  Processing  Subsystem 
for  correlation  with  target  reports  received  from  the 
BDAS. 

The  baseline  RPS  configuration  (Fig.  2)  was  simu- 
lated in  NELIAC  programming  language  for  use  on  a 
Univac  1230  computer.  The  computet  simulation  was 
designed  to  use  standard  Burroughs  variable  quantizer 
recorder  tapes  of  digitized  ASR  video  as  the  input.  A 
concurrent  analytic  simulation  effort  was  undertaken 
to  select  suitable  values  for  the  various  unspecified 
parameters  contained  in  the  system  configuration  and 
to  measure  system  performance  as  influenced  by  these 
parameters.  The  results  of  this  parallel  effort  were 
used  to  evaluate  various  features  of  the  system,  and 
several  configuration  modifications  were  developed. 

The  first  modification  is  essentially  a simplification 
of  the  baseline  RPS  design.  The  second  uses  two 
separate  target  detection  channels.  The  third  presents 
the  basic  characteristics  of  an  all-digital  RPS  as  a 
guide  to  future  development  work. 

Simplified  RPS.  Figure  3 shows  the  simplified 
RPS  configuration.  Either  of  two  video  formats  may 
be  selected  for  target  processing;  log  normal  or  linear 
MTI.  Both  are  processed  by  a rank  quantizer. 

The  hit  processing  blocks  function  in  the  same 
manner  as  in  the  original  system  design.  The  final 
detection  block  consists  of  an  azimuth  sliding  window 
that  looks  at  17  azimuthally  adjacent  radar  cells  and 
counts  the  number  of  cells  that  contain  a hit.  This 
number  is  compared  to  a threshold  on  a real-time 
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Pig.  3 Recommended  simplified  RPS  lonfifunition. 
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Fig.  4 Proposed  modified  RPS  configuration. 


basis.  This  threshold  m is  variable  and  is  set  by  the 
m/\l  estimation  block,  which  estimates  the  azimuthal 
correlation  present  in  each  zone  (2.8°  by  2 nmi) 
within  the  radar  coverage  area.  (If  m = 9,  then  at 
least  9 of  the  17  adjacent  cells  must  contain  a hit  in 
order  to  pass  the  second  threshold.) 

The  second  threshold  for  each  zone  is  stored  in  an 
m/l  7 map,  which  has  a storage  location  for  a second- 
threshold  value  for  each  of  the  4096  zones  that  con- 
stitute the  coverage  area.  This  map  is  also  used  to 
determine  which  quantizer  will  be  selected.  A ground 
map  indicates  areas  where  MTI  video  will  always  be 
selected  in  order  to  eliminate  radar  returns  from  large 
stationary  objects  (buildings,  mountains,  etc.). 

This  system  configuration  represents  a considerable 
savings  in  hardware,  software,  and  complexity  over 
the  baseline  configuration,  without  sacrificing  system 
performance. 

Modified  RPS.  Figure  4 shows  the  proposed  modi- 
fied RPS  configuration.  Again  the  two  basic  video 
formats  are  available  for  processing,  but  instead  of 
selecting  one  of  the  video  formats  as  the  input  to  the 
final  detection  block,  both  formats  are  processed  in 
parallel  through  two  independent  final  detection 
blocks.  Three  maps  are  required  (two  w/17  maps 
and  one  ground  map  not  shown  in  the  figure),  since 
the  second-threshold  m estimates  used  in  the  final 
detection  are  independent  for  the  log  normal  and 
linear  MTI  video  input  channels.  For  certain  values 
of  azimuthal  clutter  correlation,  no  value  of  second 
threshold  will  keep  the  number  of  false  centroids 
reported  per  scan  under  the  saturation  level  of  the 
centroid  processing  algorithm.  Hence,  a "switch”  func- 
tion is  inserted  in  the  normal  video  path  after  the 


position  of  the  centroid  has  been  determined.  The 
switch  prevents  further  normal  video  centroids  from 
being  passed  to  the  radar  software  when  the  value  of 
m selected  for  normal  video  exceeds  16/17  (16  hits 
in  17  azimuthally  adjacent  cells).  This  limiting  value 
is  an  engineering  estimate;  it  may  prove  advantageous 
to  lower  this  limit  after  further  study  or  operational 
usage. 

The  radar  software  (centroid  correlation  block) 
receives  position  estimates  from  either  or  both  posi- 
tion estimation  blocks  and  combines  them  into  a 
single  report  for  track  processing  by  the  ARTS  III 
radar  software. 

The  primary  advantage  of  this  configuration  is  that 
the  problem  of  which  quantizer  output  should  be 
selected  is  eliminated.  The  determination  of  the 
proper  second  threshold  to  be  used  is  put  on  a more 
rational  basis,  and  problems  associated  with  quantizer 
switching  are  eliminated. 

Future  RPS  Development.  During  the  current 
effort  to  improve  and  simpify  the  existing  RPS  design, 
it  is  appropriate  to  consider  the  type  of  system 
that  will  be  desirable  in  the  next  generation  design. 
Such  consideration  will  permit  a broader  perspective 
of  alternative  designs  and  will  lessen  the  probability 
of  future  development  becoming  prematurely  locked 
into  a particular  design  philosophy.  TTie  most  obvious 
trend  in  radar  processors  and  sensors  is  toward  all- 
digitized  systems  whose  basic  features  are: 

1.  Use  of  digital  MTI  as  well  as  digital  log  video' 

2.  Simultaneous  processing  of  both  video  clian- 
nels; 
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3.  Estimation  of  multilevel,  mean-level  azimuthal 
correlation; 

4.  Use  of  microprocessor  technology  to  reduce 
cost  and  complexity; 

5.  Greatly  reduced  input/output  and  processing 
requirements  on  the  general-purpose  computer 
(less  program  maintenance) ; 

6.  Effective  operation  in  correlated  clutter  envi- 
ronments; and 

7.  Improvement  in  target  resolution  and  accuracy. 

Such  a system  presents  a new  design  philosophy 

and  is  currently  under  investigation  for  use  in  various 
Department  of  Defense  programs.  It  represents  a first 


attempt  at  configuring  such  processing  for  FAA  needs 
and  should  be  viewed  as  an  indication  of  possible 
future  trends. 
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DAT — A MOBILE  DATA  ACQUISITION  SYSTEM 


A mobile  Data  Acquisition  Trailer  (DAT)  was  de- 
signed and  constructed  to  provide  a highly  sophisti- 
cated and  versatile  test  control  and  radar  data 
acquisition  facility  that  can  be  operated  at  a host 
radar  site.  Thus,  information  concerning  the  perform- 
ance of  the  radar/beacon  sensor  and  the  character- 
istics of  the  targets  it  sees,  whether  real  or  spurious, 
can  be  obtained.  The  DAT  enables  this  task  to  be 
accomplished  irrespective  of  the  degree  of  automated 
radar  data  processing  currently  operating  at  the  host 
site,  and  has  minimal  impact  on  the  normal  operations 
of  the  host  facility. 


BACKGROUND 

APL  is  working  to  improve  the  quality  of  data 
derivable  from  radar/beacon  sensors  and  is  investigat- 
ing optimal  methods  of  processing  data  from  one  or 
several  sensors.  The  performance  of  the  task  requires 
more  detailed  knowl^ge  of  the  information  provided 
by  the  sensors,  including: 

1 The  number  of  aircraft  detected  by  the  radar 
system; 

2.  The  number  of  aircraft  equippied  with  beacon 
transponders  and  detected  by  the  beacon  inter- 
rogation system; 


The  number  of  false  detections  declared  by 
the  sensors  as  a result  of  precipitation  or  land 
interference,  sensor  noise,  radio  interference, 
etc.;  and 


4.  Detailed  information  on  sensor  output  iii  the 
immediate  vicinity  of  targets  (aircraft  or  false 
detections)  of  interest. 

These  data  can  be  analyzed  to  determine  the  aircraft 
detection  capability  of  a radar/beacon  site,  to  deter- 
mine the  characteristics  of  spurious  data,  and  to  ex- 
amine the  difference  between  true  and  spurious 
detections  in  order  to  devise  methods  of  minimizing 
false  detections.  In  addition,  the  capability  of  simul- 
taneous data  collection  from  several  noncolocated 
sensor  sites  is  desirable. 

APL’s  experience  in  the  research  and  development 
phases  of  U.S.  Navy  search  and  fire  control  radar 
system  programs,  and  its  development  of  the  U.S. 
Navy  AN/SYS-1  Target  Information  Processing  Sys- 
tem, were  applicable  in  the  design  and  construction 
of  the  DAT. 


DISCUSSION 

The  DAT  equipment  is  shown  in  Fig.  1 . The  DAT 
was  designed  to  operate  with  two-dimensional  radar/ 
beacon  search  sensors  or  a tracking  radar.  Sensor 
inputs  are  recorded  directly  on  a wideband  recorder 
and  simultaneously  processed  in  one  of  several  avail- 
able modes  for  storage  on  digital  magnetic  tape. 
Simultaneous  data  collection  at  several  sites  is  accom- 
plished by  placing  wideband  video  recorders  at  each 
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Fig.  1 DAT  data  coUecdoo  equipment. 


site,  coordinating  recording  via  the  DAT  control  and 
communications  system,  and  then  playing  the  record- 
ings through  the  DAT  system  for  digital  processing. 

The  input  select  unit  chooses  the  desired  DAT 
input  from  the  available  data  sources  and  provides 
conditioning  for  all  signals  to  be  processed.  Signals 
selected  from  either  live  inputs  or  prerecorded  analog 
tape  playback  are  fed  to  the  appropriate  interfaces 
where  they  are  amplified,  conditioned,  coded,  and 
formatted  prior  to  distribution  to  the  video  processor, 
computer,  displays,  and  recording  equipment. 

The  Adaptive  Threshold  Video  Processor  (AVP) 
receives  radar  video  or  beacon  video  and  provides 
target  dt':ections  to  the  computer  for  target  tracking. 
Radar  data  require  dual  processing  before  a target  is 
declared.  The  first  process  derives  a threshold  from 
an  analog  sum,  in  range,  of  video  samples  prior  to 
and  following  the  target.  When  the  target  amplitude 
exceeds  the  first  threshold,  a first-threshold  detection 
is  declared.  The  second  process  is  an  azimuth  correla- 
tion (M/N)  performed  on  first-threshold-declared 
targets.  TTiis  requires  that  M first-threshold  crossings 
must  occur  out  of  N radar  dwells  at  the  same  range 
before  a second- threshold  crossing  is  declared.  When 
this  occurs,  the  range  of  the  detected  target  is  sent  to 
the  computer.  The  sensitivity  of  both  processes  is 
selectable  at  the  Display  and  Control  Console. 


A beacon  target  is  declared  when  beacon  video 
passes  a first  threshold  of  pulse  width  and  spacing 
criteria  and  then  the  M/N  azimuth  threshold.  An 
additional  feature  of  the  AVP  provides  video  ampli- 
tude quantizing  when  a data  collect  matrix  is  gener- 
ated by  the  computer. 

The  computer  (Honeywell  DDP-516R)  has  been 
programmed  to  provide  four  basic  modes  of  digital 
data  collection; 

1.  Tracking — The  computer  accepts  targets  de- 
clared by  the  AVP,  performs  scan-to-scan  cor- 
relation, and  classifies  target  tracks  as  fixed  or 
moving,  with  a capacity  of  256  of  each.  Track- 
position  and  velocity  data  are  stored  on  digital 
magnetic  tape.  A typical  moving  track  scenario 
generated  by  this  mode  is  shown  in  Fig.  2. 

2.  Hooked  Track — The  computer  performs  the 
same  functions  as  in  the  tracking  mode.  In 
addition,  a target  of  interest  may  be  selected 
(hooked)  from  the  Display  and  Control  Con- 
sole, causing  a computer-generated  data  matrix 
(40  samples  in  range  by  63  samples  in  azi- 
muth) to  be  centered  around  the  target. 
Quantized  video  amplitude  data  generated  by 
the  AVP  in  this  matrix  are  collected  by  the 
computer  and  stored  with  the  normal  track 
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Fig.  2 Scenario  of  selected  tracks  seen  by  radar  during 
a J-min  interval  (rings  signify  10-nmi  intervals). 


data.  A portion  of  typical  matrix  data  is  shown 
in  Fig.  3. 

3.  Video  Quantizing — ^This  mode  provides  for 
the  collection  of  quantized  video  over  a large 
area  of  interest.  A matrix  generated  by  the 
computer  is  manually  positioned  at  the  Display 
and  G>ntrol  Console.  The  matrix  contains 
either  80  or  160  samples  in  range  and  up  to 
125  dwells  in  azimuth  and  is  stored  on  digital 
tape. 

4.  Track  Radar — This  mode  accepts  track  radar 
measured  range,  azimuth,  and  video  amplitude 
data  and  uses  the  track  radar  coordinates  to 
center  a search  radar  data  matrix.  Data  from 
both  radars  are  stored  on  digital  tape. 

Computer-controlled  modem  communications  with 
APL  can  be  used  to  transmit  real-time  track  data  for 
immediate  analysis  and  display.  Another  computer 
function  provides  on-site  verification  of  DAT-gener- 
ated digital  magnetic  tape  data  via  TV  display. 

Data  from  the  tracking  mode  can  be  analyzed  to 
determine  the  number  of  aircraft  in  the  sensor  space. 
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Fig.  3 A portioa  of  the  data  collect  matrix  for  search  radar  video  from  bird  targets. 


regions  of  the  space  where  the  sensor  is  blind  (by 
flying  a controlled  test  aircraft  through  the  sensor 
space  and  determining  the  regions  where  it  is  not 
detected),  the  probability  of  detecting  aircraft  within 
the  sensor  space,  and  various  data  regarding  an  air- 
craft or  group  of  aircraft  within  the  sensor  space 
(pjosition,  spaced,  and  heading). 

Data  from  the  video  quantizing  mode  can  be  used 
for  detailed  analysis  of  the  video  in  the  immediate 
vicinity  of  targets  of  interest.  Useful  parameters  that 
can  be  determined  include  the  range  and  azimuth 
correlation  of  the  video,  the  mean  amplitude  of  target 
signals,  and  other  statistical  data.  TTiis  information 
can  be  used  to  develop  new  methods  of  detecting  air- 
craft in  the  video  and  reducing  false  detections. 

The  hooked  track  mode  provides  all  of  the  tracking 
mode  data  plus  a reduced  video  quantizing  capability. 
It  has  the  advantage  of  accurately  centering  a com- 
puter positioned  data  matrix.about  a target  as  it  moves 
through  the  sensor  space  and  thus  provides  a detailed 
track  life  history  of  the  motion  and  video  signature 
of  a target. 

The  track  radar  mode  can  be  used  in  accuracy  stud- 
ies by  comparing  data  from  the  track  and  search 
sensors.  Also,  since  the  tracking  radar  provides  con- 
tinuous data  on  a target,  the  characteristics  of  target 
video  as  a function  of  time  can  be  described  more 
accurately. 

The  analog  recorder  used  in  the  DAT  is  the 
AN/SPH-1  Radar  Video  Recorder/Reproducer  Set 
(RAVIR),  which  was  designed  and  developed  by 
APL  for  the  U.S.  Navy.  When  recording  live  data,  it 
provides  a fail-soft  mode  of  data  collection  in  the 
event  of  a failure  in  one  of  the  digital  processing 
components.  Furthermore,  if  data  must  be  collected 
at  several  noncolocated  sites,  only  a recorder  need  be 
installed  at  the  remote  sites  since  the  recorded  analog 
tapjcs  can  be  digitally  processed  in  the  DAT. 


SMALL-CAR  STUDY 


A study  has  b««n  complatad.  In  conjunction  with 
the  Center  for  Metropolitan  Planning  and  Research 
of  The  Johns  Hopkins  University,  that  addressed  the 
problems  associated  with  the  transition  to  a smaller 
private  automobile.  The  results  should  be  helpful  In 
guiding  policy  for  energy  conservation  In  private 
automobile  usage. 


The  Display  and  Control  Console  serves  as  the  pri- 
mary means  of  monitoring  and  controlling  the  jjer- 
formance  of  the  processing  equipment.  By  means  of 
plan  position  indicator  display  and  switching  features, 
a wide  variety  of  data  can  be  displayed.  It  has  the 
capability  of  identifying  aircraft  equipped  with  bea- 
con transponders  and  selectively  displaying  their 
positions.  Tlie  console,  combined  with  the  communi- 
cations system  (not  shown  in  Fig.  1)  enables  the 
DAT  to  be  a data  collection  command  post  while 
controlling  test  aircraft  and  mobile  land  units.  Be- 
cause the  DAT  is  a self-contained  unit  (and  self- 
supporting,  if  required),  and  because  of  the  minimal 
equipment  installation  requited  at  remote  sites  in 
multisite  tests,  the  DAT  can  coordinate  large-scale 
data  collection  efforts  with  minimal  impact  on  the 
operations  and  procedures  of  the  sensor  user. 

The  DAT  has  been  used  in  several  data  analysis 
efforts  by  the  U.S.  Navy  and  Federal  Aviation  Ad- 
ministration. One  such  effort  involved  the  analysis  of 
the  radar  video  of  birds  (Fig.  3)  in  order  that  they 
not  be  mistaken  for  aircraft  or  obscure  the  detection 
of  aircraft  during  pjeriods  of  heavy  bird  migration. 
Every  mode  of  DAT  data  collection  was  used  in  this 
task,  and  subsequent  analysis  resulted  in  several 
recommended  methods  of  rejecting  bird-generated 
false  detections.  Another  effort  involved  the  collection 
of  data  at  several  sensor  sites  that  had  overlapping 
coverage.  The  DAT  controlled  simultaneous  data  col- 
lection at  each  site  and,  by  analysis  of  tracking  mode 
information,  APL  was  able  to  quantify  the  benefits 
achievable  by  processing  the  combined  data  from  sev- 
eral noncolocated  sensors. 
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BACKGROUND 

In  examining  ways  to  alleviate  the  world  energy 
shortage,  the  private  automobile  has  been  singled  out 
as  one  of  the  most  likely  targets  for  energy  conserva- 
tion, mainly  because  of  its  relatively  large  and  ineffi- 


cicnt  use  of  energy.  Various  methods  are  available  to 
improve  the  fuel  efficiency  of  the  private  automobile. 
This  study  focuses  on  a very  promising  method  to 
decrease  curb  weight  and,  implicitly,  size. 

Three  questions  arose  at  the  outset  of  the  investiga- 
tion; How  small  is  small?  How  do  you  sell  such  a 
car?  What  effects  will  result  from  a major  shift  to 
small  cars?  Although  we  have  no  firm  answers  to  any 
of  these  questions,  the  study  has  raised  a number  of 
more  restrictive  questions  that  can  lead  to  a better 
understanding  of  the  problems  associated  with  an 
orderly  transition  ioward  a smaller  private  automo- 
bile. 


DISCUSSION 

How  Small  Is  Small?  Since  fuel  economy  im- 
proves as  weight  is  reduced,  the  question  really  is: 
How  small  a car  will  satisfy  safety,  emission,  comfort, 
aesthetic,  and  performance  criteria  acceptable  to  the 
consumer  and  to  current/proposed  government  regula- 
tions? Three  small  cats  weighing  2000,  1500,  and 
1000  lb,  delivering  25,  35,  and  45  mi/gal,  respec- 
tively, are  addressed.  The  findings  suggest  that  cars 
weighing  1500  lb  or  less  would  be  acceptable  with 
respect  to  performance  and  comfort.  A major  shift 
to  such  cars  by  1977  could  save  approximately  1 
million  barrels  of  oil  per  day  by  1980. 

Since  several  foreign  cars  weighing  less  than  2000 
lb  are  already  in  growing  use  for  both  urban  and 
highway  travel,  it  seems  clear  that  a small,  limited- 
performance  vehicle  weighing  less  than  1 500  lb  would 
be  acceptable  for  urban  use.  The  most  important  prob- 
lem is  safety.  A major  shift  to  small  cars  could  lead 
to  a 15%  increase  in  fatalities  or  serious  injuries 
unless  it  is  accompanied  by  the  institution  of  safety 
measures  such  as  improvement  in  and  more  extensive 
use  of  passenger  restraint  systems,  improved  structural 
design,  and  strict  enforcement  of  sp^  limits. 

Future  emission  goals  also  represent  an  unknown, 
in  that  enforcement  of  the  presently  scheduled  reduc- 
tion of  the  NO,  emission  standard  in  1978,  the  real 
need  for  which  is  being  questioned,  could  nullify 
efforts  to  conserve  fuel  and  hold  car  prices  down.  In 
this  area,  it  is  suggested  that  the  current  California 
standards  be  adopted  as  the  Federal  standards  and 
that  further  reduction  not  be  considered  for  a period 
of  five  years. 

There  do  not  appear  to  be  any  technological  reasons 
why  a small,  low-emission,  economical,  and  safe  pri- 
vate automobile  with  a curb  weight  of  less  than  1 500 
lb  cannot  be  mass-produced.  To  accomplish  this 


objective  will  require  substantial  design  effort  and  the 
introduction  of  specific  legislative  incentives  and  con- 
straints. 

How  Do  You  Sell  It?  Current  automotive  in- 
dustry marketing  is  directed  toward  making  the  con- 
sumer aware  of  fuel  economy  as  evidenced  in  recent 
advertisements.  The  methods  considered  in  this  report 
are  specifically  directed  toward  those  that  can  be  for- 
mulated into  legislation  by  the  Federal  Government. 
An  attempt  to  compile  from  current  literature  what  is 
known  regarding  price  elasticities  (i.e.,  a measure  of 
the  effect  of  a price  change  on  demand  for  a product) 
has  shown  that  the  pertinent  long-run  elasticities  vary 
over  a considerable  range  so  that  the  long-run  market 
response  to  increased  costs  due  to  taxes  and/or  rebates 
cannot  be  predicted  with  great  certainty.  It  is  felt 
that  a shift  to  small  cars,  through  Government  policy, 
might  best  be  effected  in  the  short  run  by  tax/rebate 
policies  designed  to  raise  exponentially  the  original 
investment  cost  of  automobiles  as  their  fuel  economy 
declines,  and  that  large  increases  in  the  cost  of  gaso- 
line must  occur  to  effect  a significant  reduction  in 
gasoline  consumption. 

What  Effect  Will  Result  from  a Major 
Shift  to  Small  Cars?  A major  conversion  to  small 
cars  would  improve  traffic  conditions  and  alleviate 
parking  problems.  If  systems  of  parking  segregated 
by  size  are  instituted,  a potential  for  twofold  increases 
in  both  on-street  and  off-street  parking  capacities  is 
foreseen.  Other  effects  on  the  national  economy  are 
difficult  to  assess,  and  a meaningful  analysis  would 
require  better  data  than  are  currently  available.  A 
"first  approximation”  analysis  indicates  that  effects  on 
the  suppliers  of  raw  materials  for  new  cars  could  be 
great — e.g.,  an  8%  reduction  in  total  steel  sales.  A 
previous  study  led  to  the  conclusion  "that  1 percent 
change  in  new  car  weight  results  in  a greater  than  1 
percent  change  in  auto  manufacturing  employment" 
(Ref.  1).  For  these  reasons,  we  recommend  that  a 
detailed  study  be  initiated  using  such  economics 
methods  as  input-output  analysis  in  order  to  account 
for  the  interrelationships  of  the  various  economic 
sectors. 

Safety.  Three  fundamental  factors  are  involved  in 
determining  the  crashworthiness  of  a vehicle: 

1.  A definable  upper  limit  on  the  short-term 
deceleration  forces  that  can  be  sustained  by  a 
properly  restrained  human  body, 

2.  A required  stopping/crushing  distance  propor- 
tional to  the  square  of  the  initial  velocity,  and 

3.  A difference  in  the  weights  of  two  colliding 
vehicles  that  will  amplify  the  deceleration 
forces  imposed  on  the  lighter  vehicle. 


The  last  factor  is  often  used  to  explain  away  the 
safety  hazards  to  small-car  occupants;  i.e.,  there 
wouldn’t  be  a problem  if  everyone  drove  small  cars. 
However,  the  first  two  factors  refute  this  claim  because 
(a)  the  human  body  can  stand  only  a certain  level  of 
deceleration,  which,  for  a given  vehicle  speed,  dictates 
the  te<)uired  crush  distance,  and  (b)  present  small 
cars  have  inadequate  crush  distance  available  within 
their  frames  for  the  speeds  at  which  they  travel. 

Accident  statistics  show  that  the  risk  of  driver 
injury  when  two  subcompact  cars  collide  is  more  than 
twice  as  great  as  when  two  full-size  cars  collide. 
Nevertheless,  the  risk  of  injury  for  the  small-car 
driver  is  only  one- fourth  to  one-half  as  great  when 
he  collides  with  another  small  car  as  with  a large 
car. 

A significant  improvement  in  small-car  safety  can 
be  achieved  by  addressing  all  three  of  the  fundamental 
factors  listed  above  and  by  taking  action  in  the  follow- 
ing areas: 

1.  Passenger  Restraint  Systems  - Improvements  in 
use  and  effectiveness  of  restraint  systems  are 
even  more  important  for  small  cars  than  for 
large  cars.  For  a 2000-lb  car,  injury  rates  can 
be  halved  by  effective  passive  restraint  systems 
such  as  automatically  activated  belts  or  air  bags. 
The  technology  is  available  to  produce  effective 
and  acceptable  passive  restraint  systems.  It  is 
recommended  that  legislation  be  passed  to  re- 
quire manufacturers  to  provide  them. 

2.  Speed  Limits -The  "safe”  speed  for  a 2000-lb 
car  is  10  to  15  mi/h  less  than  that  of  a 4000-Ib 
car,  assuming  that  passenger  restraint  s)^tems 
are  used.  The  national  speed  limit  of  55  mi/h 
must  be  maintained  and  enforced  and  a further 
reduction  in  speed  to  45  to  50  mi/h  should  be 
imposed  on  urban  expressways  where  small- 
car  traffic  is  heavy. 

3.  Improvements  in  Small-Gir  Structural  Design  - 
Small-car  structures  must  be  designed  to  some 
degree  of  crashworthiness.  It  is  recommended 
that  a 30-mi/h  barrier  crash  be  the  goal,  with 
emphasis  on  passenger  safety  rather  than  dam- 
age reduction.  The  structure  must  provide 
adequate  crush  distance.  Additionally,  it  is 
recommended  that  the  structure  of  large  cars 
be  limited  in  front-end  stiffness  so  as  to 
reduce  their  aggressiveness  to  the  smaller  cars. 


Comfort.  Anthropometric  studies  suggest  that 
an  automobile  with  interior  dimensions  adequate  to 
provide  a commonly  accepted  level  of  comfort  for 
the  middle  90%  of  the  riding  public  can  be  designed 
for  automobile  weights  as  low  as  1 500  lb.  Maximum 
comfort  and  convenience,  on  the  other  hand,  cannot 
be  obtained  without  significant  changes  in  extant  com- 
partment configurations.  Further  reductions  in  vehicle 
weight,  to  1000  lb  for  instance,  can  be  expected  to 
reduce  the  percentage  of  the  riding  public  that  will 
accept  the  comfort  level;  however,  no  quantitative 
estimates  have  been  formulated. 

Marketing  and  Policy  Considerations.  Selec- 
tion of  an  optimal  policy  to  effect  a shift  to  small<ar 
ownership  depends  heavily  on  empirical  information 
about  the  price  and  income-elasticities  of  demand  for 
gasoline,  and  the  cross-elasticities  between  new  and 
used  cars  and  between  cars  and  other  modes  of  transit. 
Recent  empirical  studies  provide  a very  wide  range 
of  estimates,  none  of  which  includes  price  increases 
as  large  as  those  recently  experienced  in  the  U;S.  We 
therefore  have  very  little  hard  evidence  from  which 
to  work  in  estimating  the  effects  of  alternative  manda- 
tory and  price  (tax)  policies. 

On  the  basis  of  recent  observations,  the  price- 
elasticities  of  demand  for  gasoline  are  relatively  low, 
even  in  the  long  run,  and  tax  policies  must  incorporate 
fairly  large  price  increases  in  gasoline  and  in  cars 
having  low  fuel  economy,  together  with  careful  sub- 
sidy and  rebate  schemes,  to  effect  a significant  reduc- 
tion in  gasoline  consumption  without  adverse  effects 
on  income.  A shift  to  small  cars  can  best  be  effected 
by  tax  jxilicies  designed  to  raise  exponentially  the 
original  investment  cost  of  automobiles  as  their  fuel 
economy  declines.  Such  a policy  will  have  strong  em- 
ployment and  income  impacts  that  must  be  cushioned 
by  redistribution  of  the  tax  revenues  through  other 
programs. 
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INTRODUCTION 


APL  conducts  a vigorous  program  of  exploratory  and  advanced  develop- 
ment in  support  of  improved  missile  and  defense  systems.  Identifying  critical 
problems  that  may  limit  the  effectiveness  of  weapon  systems  and  applying  new 
science  and  technology  to  provide  solutions  in  a timely  fashion  can  both  in- 
crease the  technical  options  available  and  improve  the  expertise  of  the  APL 
staff  for  advanced  system  development.  APL  has  long  been  in  the  forefront 
of  aerodynamics,  structures,  propulsion,  controls,  and  guidance  technology  and 
design  (or  missile  systems.  Advanced  development  tasks  in  such  areas  as 
remote  sensor  technology  and  computer  applications  to  data  and  systems  man- 
agement are  also  conducted  to  support  APL’s  mission. 

The  items  presented  here  are  only  a small  sample  of  this  class  of  work, 
selected  for  technical  interest  and  possible  applicability  in  nondefense  areas. 
Nap-of-the-Earth  Radar  is  a concept  that  may  be  potentially  applicable  to  civil 
flight  safety.  Electronic  beamsteering  increases  the  capability  of  a parametric 
sonar  for  high-resolution  ocean  bottom  surveys.  The  earth’s  electrostatic  field 
may  provide  a simple,  cheap  system  for  obstacle  sensing  and  avoidance  for 
pilotless  or  small  piloted  aircraft  under  certain  adverse  visibility  conditions, 
e.g.,  at  night  or  in  fog  or  low  cloud  cover.  The  error  correction  encoder/error 
correction  decoder  technique  is  generally  applicable  to  high-density,  multitrack, 
magnetic  tape  recording  and  reproduction. 


NAP-OF-THE-EARTH  RADAR  (NOTER) 
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The  Nap-of-the-Earth  Radar  (NOTER)  described 
herein  was  developed  as  a demonstration  of  proof  of 
principle.  The  project  was  based  on  a feasibility  study 
of  meeting  requirements  for  helicopter  terrain  scan- 
ning and  wire  detection  with  a 95-GHz  radar.  Terrain 
avoidance  and  the  avoidance  of  collision  with  unseen 
wires  is  a critical  issue  with  military  helicopter 
operations. 


BACKGROUND 

NOTER  was  intended  to  be  capable  of  integration 
with  any  of  the  terrain  imaging  systems  currently  in 
advanced  state  of  development  by  the  U.S.  Army  to 
provide  a day  and  night  capability  of  nap-of-the- 
earth  flight.  At  the  time  the  development  was  xinder- 
taken,  APL  was  aware  of  the  Army’s  laser  terrain 
warning  system  (LOTAWS)  program.  The  objective 
of  the  NOTER  program  was  to  provide  an  alternate, 
near-term  technique  for  evaluation.  The  Army  subse- 
quently decided  to  devote  all  its  efforts  to  the  laser 
approach. 

A most  pressing  problem  for  military  helicopter  op- 
erations is  the  danger  of  collision  with  unseen  wires. 
Army  strike  reports  (Ref,  1)  list  collisions  even  with 
wires  whose  locations  were  known,  as  well  as  those 
intentionally  or  unintentionally  deployed  in  a tactical 
environment.  Discussions  with  operational  personnel 
at  the  time  of  die  feasibility  study  indicated  that  even 
very  small,  soft  wires  can  disable  a helicopter.  The 
development  of  the  effective  air-defense  systems  now 
deployed  worldwide  makes  very-low-altitude  flight 
mandatory,  which  in  turn  necessitates  the  employment 
of  wire  detection  and  avoidance  systems. 


DISCUSSION 

A study  indicated  that  a radar  system  in  the  95-GHa 
region  could  detect  small  wires,  even  in  adverse 
weather  conditions.  To  assemble  and  demonstrate  such 
a system  with  available  components  was  reasonably 
straightforward.  Although  components  are  currently 
somewhat  expensive  because  of  limited  demand,  the 
system  has  the  potential  of  modest  production  cost. 
NO'TER  has  no  critical  alignment  and  maintenance 
features  and,  accordingly,  is  well-adapted  to  the  heli- 
copter’s harsh  vibration  and  operational  environment. 
Millimeter  wave  radar  systems  are  basically  compact, 
lightweight,  and  reliable  and  can  be  mounted  on  a 
72  helicopter  as  small  as  the  Gibra.  The  operational  char- 


acteristics and  general  specifications  of  NOTER  are 
listed  in  Tables  1 and  2. 


TABLE  1 

OPERATIONAL  CHARACTERISTICS 

Detect  Army  field  wire  (V^-in.  twisted  pair)  at  range  of 
2 km  and  S/N  of  > 10  dB 

Detect  #22  wire  at  range  of  1 km  and  S/N  of  > 10  dB 

Covert  operation  as  a result  of  limited  range  and  pencil- 
beam  RF  propagation 

Useful  in  moderately  foul  weather  (e.g.,  fog,  medium  rain- 
fall) 

Field  of  view:  4°  elevation  by  200°  azimuth 

Automatic  pilot  alert  system  via  visual  indicator 

Display:  variable  persistence  CRT  giving  both  pictorial  rendi- 
tion of  entire  field  of  view  and  analog  display  of  range  to 
wires  or  other  objects. 

Human  interface:  no  special  interpretive  operator  training 
required 


TABLE  2 

GENERAL  SPECIFICATIONS 


Frequency:  95  GHz 

Transmitter:  Amperex  DX-287  magnetron 
Peak  power:  5kW 
PRF:  2 kHz 

Pulse  length:  >100  ns 

Average  power:  Sl.O  W 

Antenna  size:  18-in.  diameter 

Antenna  type:  Cassegrainian 

Antenna  scan;  Palmer  (via  nutating  subreflector) 

Radar  beamwidlh:  9-4  mrad 

Scan  time:  S:  1 s for  total  field  of  view 

Range  resolution:  ~50  ft 

Display:  CRT  variable  persistence  (Tektronix  605) 
Operational  environment:  airborne/helicopter 
Size:  nose  pod  ~6  ft* 

Weight:  nose  pod  ~85  lb;  display/processor  ~50  lb 
Power:  2S1.5  kW,  prime 
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The  current  system  has  an  easily  comprehended 
visual  pilot  display  with  an  audio/flashing-light  warn- 
ing system,  and  requires  little  or  no  operator  training. 
The  system  consists  of  a nose  pod  weighing  85  lb 
with  a volume  of  6 ft*  (Fig.  1).  The  nose  pod  con- 
tains a scanning  parabolic  dish  antenna,  the  95-GHz 
magnetron,  a klystron  local  oscillator,  and  the  balance 
of  the  radar  components  (Fig.  2).  'The  radar  per- 
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Fig.  2 Front  view  of  rmdax. 


forms  a Palmer  scan  through  an  angle  of  4°  by  200° 
centered  around  the  helicopter  velocity  vector.  This 
is  required  so  that  a near-broadside  aspect  of  the 
wire  is  always  presented,  irrespective  of  the  relative 
geometry.  Figure  3 shows  the  visibility  of  several 
different-sized  wires  in  angle  from  broadside  for 
various  antenna  diameters. 

Figure  4 is  a plan  view  plot  of  the  total  NOTER 
azimuth  scan  (200°)  used  in  procuring  a radar  image 
of  object  space  and  of  the  radar’s  instantaneous  field 


r 

I 9-mtad  BEAM.9SGHZ- 

L _x  . ._i_.  f. 

-30°  -20°  -10° 


v\azimuth  scan 

I 

1-  1 

10°  ,20°  30° 


Fig.  4 Wire  detectioa  range  versus  angular  incidence 
of  radar  beam. 


of  view.  The  angular  visibility  of  a wire  with  respect 
to  orthogonal  radar  beam  impingement  varies  in- 
versely with  range  and  directly  with  wire  size. 
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Fig.  5 NOTER  display. 


tion  of  the  display.  A variable  range  threshold  line 
is  available  for  easy  interpretation  of  object  returns 
and  proximity.  The  display,  associated  electronics,  and 
power  supplies  are  located  widiin  the  helicopter. 

The  radar  has  been  bench-tested  with  limited  field 
tests  at  APL,  and  operates  as  designed.  Flight  opera- 
tional field  testing  depends  upon  the  support  of  the 
military  services. 
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Figure  5 is  a rendition  of  what  the  operator  might 
see  on  the  system’s  variable  persistence  display.  The 
lower  half  of  the  display  represents  a pictorial  view 
of  object  space.  The  upper  half  of  the  screen  displays 
range  to  wires  or  objerts  shown  in  the  pictorial  posi- 
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ELECTRONIC  BEAMSTEERING  FOR  A 
PARAMETRIC  SONAR 


An  Innovativt  digital  baamstaaring  and  signal- 
ganaration  system  for  a parametric  sonar  has  bean 
dasignad,  tabricatad,  and  tasted  by  APL.  In-water  test- 
ing of  the  system  demonstrated  Its  versatility  and 
marked  the  first  time,  to  our  knowledge,  that  a para- 
metric sonar  has  bean  steered  electronically  over  such 
a wide  range  of  angles. 


BACKGROUND 

For  the  past  several  years,  the  Development  En- 
gineering Division  of  the  Naval  Oceanographic  Office 
(NAVOCEANO)  has  been  developing  a parametric 
sonar  for  high- resolution  bathymetric  surveillance.  In 
order  to  stabilize  the  sonar  beam  on  a spot  at  the 
bottom  of  the  ocean  as  the  survey  ship  pitches  and 
rolls,  it  was  desirable  to  develop  an  electronic  beam- 
steering  system  for  the  sonar  capable  of  steering  in 
two  dimensions.  Thus,  in  the  fall  of  1973  NAV- 
OCEANO resjuested  the  assistance  of  APL  in  devel- 


oping portions  of  a beamsteering  and  signal-genera- 
tion system  for  the  sonar. 

APL  subsequently  has  designed  and  fabricated  an 
experimental  digital  beamsteering  and  signal -genera- 
tion system  and  simultaneously  designed  and  procured 
a 61-element  transducer  array  (Fig.  1)  from  Gould, 
Inc.,  to  be  used  for  demonstrating  the  system.  Fabri- 
cation of  the  array  and  the  associated  beamsteering 
electronics  was  completed  in  Fiscal  Year  1975.  Suc- 
cessful in-water  tests  of  the  integrated  system  were 
conducted  in  June  1975. 


DISCUSSION 

In  a parametric  sonar,  the  nonlinear  compressibility 
of  water  or  another  medium  is  utilized  to  generate 
difference  frequencies  between  two  or  more  high- 
frequency  transmitted  sonar  beams.  The  advantage  of 
this  arrangement  is  that  the  beamwidth  at  the  differ- 
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Fig.  1 Sixty-oae  element  experimenul  array;  faceplate 
and  individual  element  are  shown. 


ence  frequencies  can  be  made  as  narrow  as  that  of 
the  transmitted  higher-frequency  beams. 

The  sound  pattern  in  the  water  consists  of  regions 
of  high-pressure  fluctuations  where  there  is  construc- 
tive interference  of  the  two  signals  and  regions  of 
little  sound  intensity  where  there  is  destructive  inter- 
ference. These  regions  propagate  along  the  beam  at 
the  velocity  of  sound.  In  the  regions  where  the  sound 
intensity  is  high,  the  average  density  of  the  water, 
averaged  over  several  cycles,  is  reduced  from  the 
ambient  because  of  the  nonlinearity  of  the  compressi- 
bility of  the  water.  Consequently,  there  is  a pattern  of 
moving  density  changes  that  correspond  to  the  density 
and  pressure  variations  of  the  desired  low-frequency 
(or  secondary)  beam.  The  pattern  generates  the  actual 
transmitted  low-frequency  beam,  and  it  is  easy  to  see 
why  the  beamwidth  is  the  same  as  that  of  the  higher- 
frequency  beam. 

Recognizing  that  the  sound  envelope  pattern  in  the 
water  determines  the  characteristics  of  the  low- 
frequency  beam,  an  innovative  technique  was  devised 
at  APL  to  create  this  envelope  pattern  in  the  water  by 
modulating  a smg/f  carrier  frequency  on  and  off.  In 
this  way,  the  beamsteering  and  the  beam  modulation 
are  separate  functions.  It  is  not  necessary  to  have  a 
linear  system  transmit  two  frequencies  through  the 
beamsteering  and  power  amplifier  circuits  to  generate 
the  sound  in  the  water.  It  is  necessary  only  to  use  a 
single  high-frequency  sonar  beam  steered  in  the  de- 
sit^  direction  and  then  to  modulate  the  transmitting 
amplifiers  to  generate  the  sound  pattern  desired  in  the 
water.  This  permits  the  use  of  efficient  amplifiers  and 
simplified  beamsteering  circuits  for  a phased  array, 
since  the  beam  direction  may  be  determined  by  phase 


control  (and  this  for  a single  frequency)  rather  than 
time-delay  control.  Only  the  pulse  shaping  at  the  de- 
sired difference  frequency  need  be  done  with  time- 
delay  modulation  control. 

Once  the  principles  of  the  APL  signal -generation 
and  beamsteering  techniques  were  established,  the  tech- 
niques were  implemented  using  digital  electronics. 
A square  wave  was  used  for  the  high-frequency  drive 
signal  so  that  it  could  be  represented  solely  by  a series 
of  I’s  and  0‘s  in  digital  logic.  For  the  same  reason, 
pulsewidth  modulation  rather  than  amplitude  modu- 
lation was  used  to  modulate  the  drive  signal.  Dur- 
ing pulsewidth  modulation,  the  power  amplifiers  are 
turned  on  and  off  with  a varying  duty  cycle  so  that  the 
distribution  of  sound  energy  in  the  seawater  closely 
approximates  that  obtained  with  amplitude  modula- 
tion. 

Implementation  of  the  beamsteering  was  achieved 
by  passing  the  drive  signal  through  a two-dimensional 
matrix  of  shift  registers  controlled  by  two  variable- 
frequency  clocks.  The  phase  of  the  drive  signal  dif- 
fered between  each  port  of  the  shift-register  matrix 
by  a known  and  controllable  amount.  The  phase- 
delayed  signals  from  each  port  of  the  shift-register 
matrix  were  then  applied  to  the  proper  elements  of 
the  array  to  steer  the  beam  in  pitch  and  roll. 

After  a preliminary  checkout  of  the  6l-element  ar- 
ray and  the  supf>orting  electronics,  a series  of  in-water 
tests  was  conducted  at  the  Millstone  Test  Facility  of 
NUSC  near  New  London,  CT.  The  APL  digital  sig- 
nal generator  was  compared  with  a conventional 
sinusoidal  signal  generator,  array  beam  patterns  were 
measured  at  various  frequencies,  and  nonlinear  mixing 
was  demonstrated.  The  APL  signal  generator  produced 
the  equivalent  of  two  high-frequency  sound  beams  in 
the  water  (Fig.  2).  Furthermore,  as  desired,  the  pat- 
terns did  not  differ  significantly  from  those  obtained 
with  a simple  sine-wave  drive.  Finally,  the  APL  beam- 
steering system  was  shown  to  function  well  by  steering 
the  primary  beams  in  fine  increments  between  ±75®. 
A typical  beam  pattern  is  shown  in  Fig.  3 for  a roll 
steering  angle  of  30*. 

Beam  patterns  were  also  recorded  at  the  difference 
frequency  during  the  tests  even  though  the  maximum 
available  range  at  Millstone  was  too  short  to  allow 
full  development  of  the  nonlinear  beam.  Neverthe- 
less, the  beam  patterns  of  the  nonlinear  beam  shown 
in  Fig.  4 display  the  narrow  beamwidth  characteristic 
of  parametric  patterns.  Electronic  beamsteering  of  the 
difference  frequency  was  also  demonstrated  during  the 
tests,  as  was  chirping  of  the  difference  frequency  by 
varying  the  modulation  frequency  of  the  carrier. 

Based  on  the  Millstone  tests,  the  unique  techniques 
devised  at  APL  for  signal  generation  and  beamsteer- 
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(a)  24.2-kHz  PATTERN 


30°  ROLL  STEERING  ANGLE 

Fig.  3 Roll-steering  bean,  faftern  obtained  with  un- 
modulated signal  generKtoi. 


ing  of  a parametric  .sonar  were  shown  to  work  well  in 
practice  as  well  as  in  principle.  The  system  provided 
good  tjuality  steerable  sonar  beams  at  both  primary 
frequencies,  thus  resulting  in  a well-developed  steer- 
able nonlinear  beam  in  the  water. 


VOC-B5.0V 

Fig.  4 Difference-frequency  beam  pattern  obtained 
during  high-power  mi^ulated  operation  of  signal  gen- 
erator. 


(b)  19.2-kHz  PATTERN 

Fig.  2 Broadside  beam  patterns  obtained  with  modu- 
lated signal  generator. 


THE  EARTH'S  ELECTROSTATIC  FIELD  AS 
AN  OBSTACLE  SENSOR 


The  perturbations  in  the  atmospheric  electric  field 
in  the  vicinity  of  orographic  protrusions,  as  computed 
by  Poisson’s  equation,  have  been  experimentally  veri- 
fied for  a limited  range  of  parameters.  These  atmos- 
pheric field  perturbations  can  be  used,  at  least  for  a 
limited  group  of  vehicles  and  protrusions,  for  terrain 
avoidance  in  "fair-weather"  conditions. 


BACKGROUND 


Using  relatively  simple  electronics  and  servomecha- 
nisms as  an  autopilot,  M.  L.  Hill  (Ref.  1)  showed 
that  it  was  possible  to  fly  a remotely  piloted  vehicle 
(RPV)  along  the  equipotential  contours  that  exist  in 
the  atmosphere,  parallel  to  the  earth's  surface.  This 
"electrostatic  stabilization"  of  flight  over  level  terrain 
has  been  shown  to  be  possible  whenever  the  so-called 
"fair-weather”  electric  atmospheric  field  exists,  typic- 
ally about  90%  of  the  time.  Since  the  conductivity  of 
the  ground  (or  any  orographic  protrusion)  is  so  much 
greater  than  that  of  the  atmosphere,  the  boundary 
condition  for  the  Poisson’s  equation  describing  the 
atmospheric  potential  is  that  the  value  of  the  poten- 
tial at  the  earth's  surface  and  over  the  surface  of  any 
protrusion  is  a constaiit  conventionally  taken  as  zero. 
Thus  in  the  vicinity  o'  all  protrusions,  natural  or  man- 
made, the  atmospheric  field  will  be  distorted  com- 
pared with  the  field  that  would  exist  for  a flat  earth. 
The  possibility  of  utilizing  the  atmospheric  electric 
field  for  obstacle  avoidance  depends  upon  the  detec- 
tion and  use  of  these  spatial  disturbances. 


DISCUSSION 


The  first-order  attempt  at  a model  for  the  atmo- 
spheric electric  neld  posits  a conducting  layer  at  about 
90  km  altitude  that  is  maintained  at  a positive  poten- 
tial relative  to  the  earth  of  about  300  000  V by 
thunderstorm  pumping.  The  spherical  capacitor  formed 
by  the  electrosphere  and  the  earth  is  discharged  by 
^air-weather  leakage  current  censitics  of  about  2 X 
lO"’*  A/m*.  Under  fair-weather  conditions,  the  cur- 
rent density,  J,  is  constant  with  altitude  and  the  con- 
ductivity, A,  increases  with  altitude.  Since  J = AE, 
where  E is  the  electric  field,  the  electric  field  has  a 
divergence,  V'E  “ p/t„,  where  p is  the  charge 
density.  Thus  the  atmospheric  potential,  is  given 
by  Poisson’s  equation 

W = -p/to- 


If  it  is  assumed  that  the  conductivity  is  a simple  ex- 
ponential function  of  the  altitude,  z,  or  the  form 
A = A„e^,  it  can  be  shown  (Ref.  2)  that  the  atmo- 
spheric potential  is  a solution  of 


Vv-P^|^  = o. 


Computer  and  analytic  solutions  of  Eq.  (1)  are 
discussed  in  Refs.  2 and  3.  Figure  1 shows  a compu- 
ter solution  of  Eq.  (1)  for  a 15-m  wall  (Ref.  3).  It 
shows  the  contours  of  constant  potential  under  the 
realistic  assumption  that  the  electric  field  strength  at 
the  surface  of  the  earth  is  100  V/m.  These  contours 
show  the  typical  behavior  of  the  electric  field  near  a 
protrusion;  as  one  approaches  the  wall,  the  horizontal 
component  of  the  atmospheric  field  increases  from 
zero.  The  vertical  component  decreases  as  one  ap- 
proaches the  wall  near  the  ground  but  is  much  higher 
than  the  unperturbed  field  directly  above  the  wall. 
Figure  2 shows  an  analytic  calculation,  also  from  Ref. 
3,  of  the  potential  contours  near  a cross  section  of 
South  Mountain,  Maryland.  The  horizontal  com- 
ponent of  the  electric  field  is  5.5  V/m  measured  at  a 
distance  five  times  the  mountain  height  at  an  altitude 
one-half  the  mountain  height. 

The  computations  all  assume  an  idealized  fair-wea- 
ther atmospheric  field  with  a vertical  scale  factor  p, 
as  contained  in  Eq.  (1),  of  2.5  X 1(H  m-*.  In  the 
study  of  atmospheric  electricity,  the  term  fair  weather 
refers  to  occasions  when  the  electric  field  vector  is 
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Fig.  1.  Potential  contours  about  a wall  15m  high. 
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Fig.  2 Equipotential  surfaces  near  South  Mountain. 


pointed  toward  the  earth  and  there  is  no  precipita- 
tion or  electrical  storm  activity  in  the  area.  However, 
during  times  of  meteorological  fair  weather,  the 
atmospheric  field  can  reverse  and  is  always  subject  to 
significant  fluctuations.  These  fluctuations  are  discussed 
in  Ref.  4;  their  possible  effects  on  obstacle  avoid- 
ance are  discussed  in  Ref.  3.  As  Ref.  3 shows,  the 
fair-weather  field  exists  at  least  90%  of  the  time  over 
much  of  the  earth,  and  fluctuations  in  the  horizontal 
component  of  the  field  are  rare  near  ground  level. 

The  measurement  of  the  atmospheric  electric  field 
or  potential  is  difficult  because  the  atmospheric  resis- 
tivity is  typically  5 X 10**  n/m.  For  the  experimental 
work  reported  here,  three  methods  of  measuring  the 
atmospheric  field  were  used  (Ref.  3).  For  ground- 
based  measurements,  a very  high  impedance  volt- 
meter, coupled  to  the  atmospheric  field  with  a radio- 
active ecjudizer,  and  a field  mill  were  used.  The  field 
mill  measures  the  electric  field  by  means  of  the  rela- 
tionship E = a/«o  where  <r  is  the  bound  charge  on  a 
conductor  that  is  periodically  exposed  to  the  atmo- 
spheric field.  Measurements  from  a Cessna  337  aircraft 
used  both  the  high-impedance  voltmeter  and  the 
low-impedance  electric  field-sensing  circuit  described 
in  Ref.  1.  All  of  these  instruments  are  difficult  to 
use.  From  an  aircraft,  there  is  the  additional  prob- 
lem caused  by  the  charging  of  the  aircraft  surface 
relative  to  the  ambient  atmosphere.  Methods  of  mini- 
mizing or  eliminating  these  difficulties  are  discussed  in 
Ref.  3. 

As  discussed  ia  Ref.  3,  ground-based  measure- 
ments of  field  perturbations  about  obstacles  are  best 
carried  out  using  a field  mill.  Figure  3 shows  the 
change  in  the  vertical  and  horizontal  components  of 
the  atmospheric  field  as  a stepladder  covered  with  a 
78  conducting  mesh  is  moved  toward  a field  mill 


Fig.  3 Vertical  and  horizontal  electric  field  versus 
time;  measurements  made  with  a field  mill  at  vaiiou) 
distances  in  front  of  a 2.4-m  laddM  covered  wi^  con- 
ducting mesh. 
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Fig.  4 Lower  wing  differencial  signal  versus  liffle; 
flights  south  and  nonh  along  South  Mountain  at  1300 


mounted  1 m above  the  earth.  The  results  conform 
with  those  predicted  from  Fig.  2. 

Figure  4 shows  the  potential  difference  measured 
wingtip  to  wingtip  for  an  instrumented  Cessna  337 
aircraft  flying  parallel  to  South  Mountain.  Measured 
horizontal  gradients  of  13  to  20  V/m  agree  well 
with  the  predictions  of  Fig.  2 for  a flight  path  just 
below  the  mountain  peak  and  200  to  300  ft  from  the 
slope.  Since  a variation  in  the  horizontal  field  of  5 
V/m  can  be  detected,  field  disturbances  for  this  type 
of  mountain  can  be  detected  at  horizontal  distances 
of  five  times  the  mountain  height  for  horizontal  flight 
paths  lying  below  the  peak  of  the  ridge. 
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ERROR  CORRECTION  ENCODER/ 
ERROR  CORRECTION  DECODER 


A high-density  digital  tape  recording/reproducing 
systam  planned  for  use  in  one  of  APL’s  programs  has 
an  inherent  bit  error  rate  (BER)  on  the  order  of  10-*. 
APL  designed  and  developed  an  Error  Correction  En- 
coder (ECE)  and  an  Error  Correction  Decoder  (ECD) 
that  reduce  the  effective  BER  to  a value  on  the  order 
of  The  lowered  BER  is  achieved  by  generating 

coded  information  (parity  bits)  from  the  input  data, 
recording  that  information  on  tracks  that  would  other- 
wise be  unused,  and  then,  on  playback,  decoding  the 
information  and  using  it  to  correct  errors  in  the  data. 


BACKGROUND 

In  the  process  of  recording  high-density  pulse 
coded  modulation  (pcm)  data  on  magnetic  tape  and 
subsequently  playing  back  the  data,  the  tape  recorder/ 
reproducer  and  the  tape  itself  may  introduce  errors. 
Such  errors  are  caused  by  the  practical  limits  of  tape 
recorder  performance  and  minute  imperfections  in  the 
magnetic  tape.  The  measure  of  these  errors,  the  ber, 
is  the  number  of  errors  that  occur  in  a given  number 


of  bits.  The  high-quality,  14-track  tape  recorder/ 
reproducer  planned  for  use  in  an  APL  program 
records  data  at  a density  of  30  kilobits  per  inch  and 
has  a BER  on  the  order  of  10~®.  In  this  program,  data 
will  sometimes  be  received  from  only  8 or  12  input 
channels,  thus  leaving  six  or  two  tape  tracks  free  for 
error  correction  information. 

An  innovative  encoding  and  decoding  technique  was 
derived  for  accomplishing  error  correction  by  using 
the  otherwise  free  tracks.  The  technique  has  been  im- 
plemented in  the  ECE  and  the  ECD.  The  ECE  and 
the  ECD  will  lower  the  overall  ber  to  10**^. 


DISCUSSION 

Figure  1 i;  a simplified  block  diagram  of  the  sys- 
tem in  which  the  ECE  and  the  ECD  are  designed  to 
operate.  For  simplicity,  control  lines  are  omitted. 
Fourteen  data  channels  (0  to  13)  can  be  accommo- 
dated although  error  correction  is  possible  only  if  8 


1 
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or  12  are  used.  The  record  digital  interface  converts 
the  data  to  a recording  code  and  controls  the  tape 
recorder  speed  synchronization.  The  playback  digital 
interface  controls  the  tape  recorder  playback  speed, 
synchronizes  and  deskews  the  data,  and  checks  the  data 
for  dropouts  by  means  of  signal  monitoring  elec- 
tronics. As  a result  of  this  last  operation,  the  play- 
back digital  interface  generates  data  quality  informa- 
tion (</„  to  dii)  on  a line-for-line,  bit-for-bit  basis. 

The  ECE  and  the  ECD  operate  at  rates  up  to  10 
MHz.  Several  different  modes  are  available,  depend- 
ing on  the  number  of  active  input  channels  (8,  12, 
or  14)  and  the  method  of  correction.  (If  all  14 
channels  of  input  data  are  present,  then  the  ECE  and 
ECD  merely  pass  the  data  straight  through. ) 

The  addition  of  an  external  multiplexer  and  de- 
multiplexer would  allow  handling  data  at  120  MHz 
from  a single  input  channel,  but  would  require  a 
42-track  tape  recorder/reproducer  for  subsequent 
recording. 

Recording  Operations.  If  eight  channels  of  in- 
put data  are  present  (14,8  mode),  six  modulo-2 
adder  circuits  in  the  ECE  generate  6 parity  bits  for 
each  8-bit  word.  In  this  case,  the  data  appear  in 
channels  V,  to  I/,,  and  the  parity  bits  in  channels 
Uo  to  U2  and  I/,,  to  U„.  Each  14-bit  word  from  the 
ECE  is  complete  in  itself  and  independent  of  adja- 
cent words. 

If  12  channels  of  input  data  are  present  (14,12 
mode),  nine  modulo-2  adder  circuits  in  the  ECE 
operate  on  a block  of  four  contiguous  words,  gen- 
erating 8 parity  bits  that  ate  packed  2 in  each  word. 
The  bits  across  four  contiguous  words  in  a single 
track  constitute  a byte.  In  this  case,  the  data  appear 
in  channels  to  f/i,  and  the  parity  bits  in  channels 
I/o  Sind  t^is-  Each  14-bit  word  from  the  ECE  is  not 
independent  in  itself  (so  far  as  parity  is  concerned), 
but  rather  the  four-word  block  is  an  independent 
unit. 

Reproducing  Operations.  On  playback,  the 
ECD  can  use  the  data  and  parity  bits  to  detect  errors 


TABLE  1 

ECE/ECD  CORKECnON  CAPABIUTIES 


Mode  and 
Technique 

Type  of  Dau  Errors  Corrected 

(14,8) 

Error  correction 

Three  adjacent  track  errors 
Two  adjacent  track  errors 
Two  track  errors  separated  by  one  good 
track 

Single  track  errors 

Double  errors,  parity  tracks  0 and  1.1  (no 
error  indication) 

(14,8) 

Erasure  correction 

Single  track  error,  where  the  track  num- 
t^r  matches  the  data-quality  line  indi- 
cating "bad  data ’’  or  matches  one  of  two 
data-quality  lines  indicating  “bad  data*’ 
Any  two  track  errors  where  the  track  num- 
l^rs  match  the  data-quality  lines  in 
dicating  “bad  data" 

(14,12) 

Error  correction 

Byte  errors  on  anv  single  data  track 
(a  byte  error  can  be  one  of  the  15  non- 
zero combinations  of  the  four  bits) 

(14,12) 

Erasure  correction 

Single  byte  errors  on  any  single  data  track, 
where  the  track  number  matches  the 
data-quality  line  indicating  "bad  data" 
or  the  track  number  matches  either  one 
of  two  data-quality  lines  indicating  "bad 
data" 

Two  byte  errors  on  any  two  tracks  ■where 
tbe  track  numbers  match  the  data- 
quality  lines  indicating  "bad  data" 

and  correct  the  data  by  an  "error  correction”  tech- 
nique that  is  essentially  the  reverse  of  the  encoding 
operation  performed  by  ECE.  Alternatively,  the  ECD 
can  use  the  data-quality  information  from  the  play- 
back digital  interface,  which  indicates  tracks  or  bit 
positions  likely  to  contain  erroneous  data,  to  cor- 
rect the  data  by  an  "error  erasure”  technique.  A sum- 
mary of  the  correction  capabilities  is  given  in  Table  1. 
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INTRODUCTION 


APL  requires  thorough,  up-to-date  competence  in  the  most  recent  advances 
in  basic  knowledge,  which  can  only  be  maintained  by  carrying  out  creative 
research  in  the  frontier  areas  of  physical  science.  Applied  research  and  devel- 
opment programs  draw  continually  on  advances  in  knowledge  and  under- 
standing of  basic  physical  phenomena  but  are  often  constrained  by  their 
product  or  application  orientation  from  contributing  broadly  to  the  advance- 
ment of  this  knowledge.  The  APL  Research  Center  provides  this  capability 
through  innovative  work  in  selected  areas  most  pertinent  to  the  APL  tech- 
nical programs.  These  comprise  principally  the  mechanisms  whereby  matter 
and  energy  interact  at  the  atomic  or  molecular  and  quantum  levels,  and  their 
application  to  both  microscopic  and  macroscopic  physical  and  biological 
processes.  In  general,  this  research  is  aimed  at  void  areas  between  the  often 
highly  idealized  research  problems  characteristic  of  universities  and  the  spe- 
cific application  or  product  - oriented  research  characteristic  of  industrial 
laboratories.  Much  of  the  work  is  multidisciplinary  in  order  to  fill  the  gaps 
between  the  traditional  scientific  disciplines. 

The  productivity,  originality,  and  quality  of  this  work  are  attested  by 
more  than  50  publications  per  year  in  leading  professional  journals  and  by 
many  invited  presentations.  Only  a small  sample  of  this  work  is  presented 
herein.  Other  research  - areas  not  covered  here  or  elsewhere  in  this  volume 
include  applied  mathematics,  lasers,  microwaves,  and  solid-state  physics. 


ORGANIC  CONDUCTORS 


Racent  dfscovaries  of  several  organic  compounds 
with  metallic  properties  at  room  temperature  are  of 
great  potential  technological  significance.  However, 
materials  previously  studied  have  been  characterized 
by  a sharp  metal-to-insulator  transition  at  tempera- 
tures above  those  where  a transition  to  a supercon- 
ducting state  might  be  expected  to  occur.  A new 
organic  conductor,  hexamethylene-tetraselenafulvalin- 
lum  tetracyanoquinodimethanide  (HMTSF-TCNQ), 
has  been  discovered  that  remains  metallic  to  tempera- 
tures below  i * K. 


BACKGROUND 

The  original  development  of  the  salt  of  tetrathio- 
fulvalene  and  7,7,8,8-tetracyanoquinodimethane 
(TTF-TCNQ)  raised  the  electrical  conductivity  of 
the  organic  metallic  state  by  two  orders  of  magni- 
tude (Ref.  1).  Since  that  development,  related  or- 
ganic compounds  of  high  conductivity  have  also  been 
synthesized  in  this  and  other  laboratories  (Table  1). 
The  possibility  of  creating  organic  compounds  with 
superconducting  properties  at  relatively  high  tem- 


peratures has  been  precluded  by  the  occurrence  of 
metal-to-insulator  transitions  in  all  materials  studied 
previously. 


DISCUSSION 

The  newly  discovered  salt  HMTSF-TCNQ  is  an 
addition  to  the  class  of  quasi-one-dimensional  con- 
ductors based  upon  the  prototype  TTF-TCNQ,  but 
its  properties  represent  a radical  departure  from 
those  of  previous  compounds  (Ref.  2).  It  is  the  first 
organic  material  to  continue  conducting  as  the  tem- 
perature, T,  approaches  0°K,  with  conductivity  on 
the  order  of  10’  n^'cm-’  at  1°K.  This  behavior  is  in 
sharp  contrast  to  the  other  materials  in  Table  1 whose 
electrical  conductivity,  a,  is  characterized  by  some 
iit  some  temperature  T„  and  then  a sharp  drop 
below  some  critical  temperature  T„.  Its  room-temper- 
ature conductivity  of  2200  n^’cm"’  is  the  highest  of 
any  known  organic  material. 


TABLE  I 


ORGANIC  SEMIMETALS;  1:1  TCNQ  SALTS  OF 


x:kx'. 


COMPOUND 

R 

X 

<7jj,  (O  ' »"•  ') 

Tm 

TTF 

-H 

S 

500 

20 

59 

DMTTF 

2-H.2-CH3 

s 

50 

25 

'V50 

TMTTF 

-CHs 

s 

1000 

1 

100 

-CHa 

s 

350 

15 

60 

HMTTF 

-CHsCHaCHi- 

s 

500 

4 

80 

TSF 

-H 

Se 

800 

12 

40 

DTDSF 

-H 

2S,2Se 

500 

7 

64 

TMTSF 

-CH.1 

Se 

1000 

6 

71 

HMTSF 

-CH^CHaCHa- 

Se 

2000 

No  transition 

The  molecular  structure  of  HMTSF-TCNQ  is 
illustrated  in  Fig.  1.  The  new  donor  molecule  is  the 
logical  extension  of  the  tetramethyl  analogs  TMTTF 
(Ref.  3)  and  TMTSF,  which  b«ome  insulating  at 
low  temperatures.  However,  the  larger  nonplanar 
molecule  HMTSF  is  incompatible  with  the  crystal 
structure  of  TMTSF-TCNQ.  This  leads  to  a dis- 
ordered structure  along  one  principal  direction  per- 
pendicular to  the  stadcing  axis.  The  higher  molecular 
symmetry,  reduced  coupling  between  conducting 
chains,  and  lack  of  three-dimensional  long-range  or- 
der suggest  natural  rationalizations  for  suppression 
of  the  insulating  state. 

From  a room-temperature  value  of  2200  n"‘cm*', 
the  resistivity  p(T)  drops  rapidly  as  the  sample  is 
cooled  (Fig.  2).  However,  below  110“K,  the  curve 
flattens  and  passes  through  a broad  minimum  be- 
tween 45®  and  75°K,  where  the  conductivity  is  about 
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Fig.  2 Normalized  DC  resiitivities  (solid  line)  of  nine 
samples  and  microwave  resistivities  (a)  of  two  addition- 
al samples  of  HMTSF-TCN(}.  (Inset:  DC  conductivity 
of  typiM  crystal  near  T - 0*K) 

3.4  times  the  room-temperature  value.  With  further 
cooling  the  resistivity  rises  again,  but  extrapolation 
from  T — 1®K  to  T — 0 (inset)  yields  values  still 
in  the  metallic  range. 

In  dose  agreement  with  the  DC  results  is  the  nor- 
malized microwave  resistivity  at  10  GHz,  presented 


in  Fig.  2 for  two  additional  crystals  (Refs.  4 and  5). 
The  microwave  dielectric  constant  remains  large  and 
negative,  as  for  a metal,  throughout  the  temperature 
range.  No  electron  spin  resonance  signal  is  ob^rvable 
in  this  compound  at  any  temperature. 

The  resistivity  appears  to  consist  of  at  least  two 
distinct  contributions.  Above  100°K,  the  normalized 
curves  are  highly  reproducible;  below  45®K,  varia- 
tions are  substantially  greater  where  conduction  is 
presumably  limited  by  impurities  and  lattice  defects. 
The  simplest  assumption  is  that  p(T)  is  the  sum 
of  an  instrinsic  part  that  decreases  and  an  "impurity” 
part  that  increases  as  the  temperature  is  lowered, 
forming  a minimum  at  their  crossing  point.  Recent 
data  for  the  quasi-one-dimensional  inorganic  polymer 
(SN),  invite  similar  analysis  (Ref.  6). 

The  intrinsic  resistivity  of  HMTSF-TCNQ  above 
110®  K is  described  by 

p(T)/p  (300®K)  = 4 ^T’’  (1) 

where  4 -=  0.264,  ^ - 8.84  X 10-»,  and  y - 2.39. 
The  other  compounds  in  this  class  obey  a similar  re- 
lation with  values  fory  of  2.33  (TTF-TCNQ),  2.34 
(TMTTF-TCNQ),  and  2.4  (TMTSF-TCNQ). 
From  the  trend  in  the  coefficient  i,  which  decreases 
with  molecular  weight,  it  appears  that  the  second 
term  in  Eq.  (1)  is  phonon  dominated.  For  a simple 
one-dimensional  metal,  the  usual  one  electron-phonon 
interaction  yields  a resistivity  proportional  to  T.  The 
large  value  for  y found  in  these  materials  is  not  a 
general  feature  of  one-dimensional  metals  but  rather 
is  a feature  of  two-band  organic  systems. 

HMTSF  also  forms  a conducting  salt  with  the  elec- 
tron acceptor  1 1 , 1 1 , 1 2, 1 2-tetracyano-2,6-naphtho- 
quinodim^ane  (TNAP) . Preliminary  Dc  and  micro- 
wave measurements  indicate  that  HMTSF-TNAP 
remains  metallic  to  1®K  and  dis{^ays  a resistance 
minimum  at  a somewhat  higher  temperature  and 
lower  conductivity  than  does  HMTSF-TCNQ. 
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PHOTOEXCITED  TRIPLET  MECHANISM  OF 
CHEMICALLY  INDUCED  ELECTRON  SPIN 
POLARIZATION 


A study  has  bean  mads  of  the  role  of 

photoexefted  triplet-state  molecules  in  the  production 
of  electron-spin-polarized  radicals  during  the  photo- 
decomposltion  of  carbonyl  compounds.  The  predicted 
dependence  of  the  electron  spin  polarization  on  the 
relative  orientations  of  the  electric  vector  of  a polar- 
ized photoexcitation  light  source  and  the  magnetic 
field  has  bean  observed,  thus  confirming  the  triplet 
mechanism. 


Wan  and  his  co-workers  at  Queens  University, 
Ontario,  Canada  (Ref.  5).  This  not  only  confirns^ 
the  TM  as  the  source  of  the  electron  spin  polarization 
in  the  systems  but  gave  considerable  informsyon 
about  the  electronically  excited  states  involved  and 
their  chemical  reactivit)',  which  must  be  considerably 
faster  than  was  believed  from  earlier  wet-chemistry 
experiments. 


1 
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BACKGROUND 

Chemically  induced  magnetic  polarization  is  the 
generic  name  for  a number  of  recently  discovered 
processes  whereby  free-radical  reactions  in  liquids 
yield  abnormal  populations  of  the  nuclear  spin  states 
of  the  diamagnetic  products  and  reactants,  and  ab- 
normal populations  of  the  electron  spin  states  of  the 
free-radical  intermediates  (Ref.  1 ) . TTie  polarizations, 
which  are  readily  observed  by  nuclear-magnetic-res- 
onance and  electron-spin-resonance  spectroscopy,  yield 
information  about  very  rapid  individual  steps  in  the 
overall  reaction  mechanism  that  are  difficult  to  study 
directly. 

The  origin  of  the  strong  electron  spin  polarizations 
produced  by  photodecomposition  of  aromatic  car- 
bonyl compounds  is  of  interest,  since  these  systems 
have  both  practical  and  theoretical  importance  as  pho- 
tochemical models.  One  possibility  was  the  combined 
effects  of  chemical  bonding  interactions  and  magnetic 
interactions  between  two  radicals  (radical  pair  mech- 
anism, or  rpm);  the  other  was  the  production  of 
polarization  in  a photoexcited  triplet-state  molecule 
that  then  reacted  to  give  a spin-polarized  radical  pair 
(triplet  mechanism,  or  tm)  (Ref.  2).  A definitive 
experimental  test  of  these  two  mechanisms  was  sug- 
gested by  detailed  analysis  of  the  TM  using  a time- 
dependent  quantum  mechanical  formulation  to  ac- 
count for  rotation  of  the  molecule  during  the  evo- 
lution of  the  electron  spin  states.  This  led  to  several 
interesting  conclusions,  the  most  notable  being  that 
the  electron  spin  polarization  depends  on  the  rela- 
tive orientation  of  a polarized  photoexciting  light 
and  the  external  magnetic  field  (Ref.  3).  Shortly  after 
preprints  of  Ref.  3 were  sent  to  scientists  working 
in  the  field,  the  predicted  light-polarization  depen- 
dence was  observed  in  the  photodecomposition  of  var- 
ious quinones  by  A.  J.  Dobbs  and  K.  A.  McLauchlan 
at  Oxford  University,  England  (Ref.  4)  and  J.  K.  S. 


DISCUSSION 

In  the  electronic  ground  states  of  most  molecules, 
all  electron  spins  are  piaired  (singlet  spin  states), 
each  pair  forming  one  chemical  bond.  If  an  electron 
pair  bond  is  broken,  each  of  the  resulting  free-radical 
fragments  has  one  unpaired  electron.  TTie  magnetic 
moment  associated  with  the  spin  of  this  unpaired 
electron  may  be  parallel  or  antiparallel  to  an  external 
magnetic  field;  thus,  a free  radical  has  a doublet  elec- 
tron spin  state.  In  an  excited  electronic  state  of  a 
molecule,  the  spin  of  the  excited  electron  may  be 
paired  with  the  spin  of  its  original  partner  in  the 
fully  paired  ground  state  (excited  singlet)  or  may 
be  antipaired  with  it.  In  the  latter  case,  the  total 
magnetic  moment  of  the  two  electrons  may  be  par- 
allel, antiparallel,  or  perpendicular  to  the  external 
magnetic  field,  and  the  resulting  state  is  a triplet  spin 
state. 

The  electron  spin  state  cannot  change  during  photo- 
excitation, and  so  the  singlet  excited  state  is  formed 
initially.  Usually,  however,  the  corresponding  excited 
triplet  state  lies  at  a lower  energy  level  than  the 
singlet,  and  there  is  a tendency  to  slide  "downhill" 
into  the  triplet  state.  The  change  in  electron  spin  state 
(intersystem  crossing,  or  ISC)  is  effected  by  the  spin- 
orbit  interaction,  i.e.,  the  interaction  between  the 
electron  spin  magnetic  moment  and  the  magnetic 
field  produced  by  its  orbital  motion  in  Ae  electro- 
static field  of  the  molecule.  The  interaction  is  nat- 
urally spin  selective;  however,  since  Ae  starting 
singlet  state  is  nonmagnetic,  Ae  initial  triplet  polar- 
ization must  also  be  field  independent  and  Aerefore 
resides  in  Ae  purely  intermolecular  interaction  be- 
tween Ae  magnetic  moments  of  Ae  two  electron  spins 
and  the  associated  zero-field  triplet  spin  states.  After 
the  "birth”  of  the  triplet,  the  magnetic  moments  of 
these  spins  begin  to  interact  wiA  the  magnetic  field. 
Provided  the  magnetic  field  is  strong  enough,  the 
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Fig.  1 Schematic  description  of  the  triplet  mechanism 
of  electron  spin  polarization. 


This  process  is  indicated  schematically  in  Fig  1 
where  <azrs  's  the  spin-spin  interaction  energy; 
Tj,  Ty,  and  are  the  zero-field  triplet  spin  states; 
<i>z  is  the  energy  of  interaction  with  the  external  mag- 
netic field;  and  T,,  T„,  and  T.,  are  the  2feeman  spin 
states  in  a strong  magnetic  field.  Since  the  zero-field 
and  Zeeman  spin  states  are  the  same  for  some  but 
not  all  orientations  of  the  molecule  with  respect  to 
the  magnetic  field,  the  amount  of  polarization  that 
can  be  transferred  from'  zero-field  to  Zeeman  spin 
states  depends  on  the  initial  orientation  of  the  triplet 
molecule.  This  orientation  can  be  partially  controlled 
by  using  polarized  light,  which  will  excite  only  those 
molecules  whose  transition  dipole  moment  has  a 
component  in  the  direction  of  the  light  polarization. 
This  leads  to  a dependence  of  the  observed  polar- 
ization on  the  relative  orientations  of  the  light- 
tiola'ization  axis  and  the  external  magnetic  field. 


purely  intermolecular  initial  polarization  is  partially 
transferred  to  a polarization  with  respect  to  the  mag- 
netic field  and  the  associated  triplet  2^eeman  spin 
states. 

The  strength  of  the  magnetic  field  is  important 
because  the  rate  of  transfer  of  the  polarization  is 
determined  by  the  rate  at  which  the  spin  magnetic 
moments  begin  to  precess  about  the  magnetic  field, 
thereby  competing  with  the  loss  of  the  intermolecular 
polarization  by  -random  rotational  motions  of  the 
molecules  that  average  to  zero  the  dipolar  interaction 
between  the  electron  magnetic  moments.  Finally,  the 
chemical  reaction  of  the  triplet  molecule  transfers  the 
polarization  to  a pair  of  doublet  radicals.  The  last 
step  must  occur  rapidly  (within  10^  s);  otherwise, 
the  modulation  of  the  dipolar  interaction  by  the  ran- 
dom motions  of  the  triplet  molecule  destroys  the 
spin  polarization  by  disorienting  the  steady  pre- 
cession of  the  spin  magnetic  moments  about  the  mag- 
netic field. 
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REACTION  RATES  OF  H AND  O WITH 
METHYL  HALIDES 


Rate  constants  for  raact/ons  of  the  typo  H + CHsX 
-»  HX  -h  CHs  and  0 -h  CHsX  -»  OH  -h  CHsX  fwhere 
X Is  Br,  Cl,  or  F)  have  been  maasured  over  tha  tarn- 
paratura  range  of  300  to  lOOO'K.  These  reactions  are 


Important  as  prototypes  In  tha  basic  understanding  of 
flame  Inhibition  by  halogen  compounds  and  may  also 
be  of  Intarast  In  tha  chamistry  of  tha  stratosphara 
perturbed  by  such  compounds. 
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BACKGROUND 

Among  the  many  substances  that  have  been  used  as 
chemical  inhibitors  of  flame  propagation,  some  of  the 
best  are  relatively  simple  compounds  that  contain 
halogen  atoms,  particularly  bromine.  The  usual  as- 
sumption is  that  such  materials  function  as  "radical 
traps,”  i.e.,  converting  or  diverting  free  atoms  and 
radicals  such  as  H,  O,  OH,  etc.  from  their  crucial  role 
as  chain  carriers  in  the  fast  reactions  upon  which  the 
combustion  process  depends.  The  more  detailed  under- 
standing of  such  a mechanism  will  depend  heavily  on 
reliable  knowledge  of  the  rates  of  the  various  radical- 
inhibitor  reactions.  In  another  area,  there  has  been 
much  recent  concern  over  the  effects  of  species  such 
as  the  chlorofluorocarbons  and  methyl  halides  on  the 
stratospheric  ozone  layer,  which  is  so  vital  to  life  on 
earth.  The  reaction  rates  of  these  species  with  simple 
atmospheric  atoms  and  radicals  (particularly  O and 
OH)  are  of  prime  importance  in  assessing  these 
effects. 

For  a number  of  years,  the  chemical  kinetics  re- 
search program  at  APL  has  been  actively  involved  in 
measuring  such  atom-molecule  reaction  rates  over  a 
wide  temperature  range  and  under  well-defined  condi- 
tions. The  basic  apparatus  (Fig.  1)  consists  of  a fast- 
flow  system  in  which  atoms  are  generated  by  a micro- 
wave  discharge  in  a small  concentration  of  parent  gas 
(Hj  or  Oj)  heavily  diluted  with  an  inert  gas  such  as 
helium.  The  gases  are  pumped  down  a quartz  reactor 
past  a movable  injector  through  which  a stable  reac+ant 
gas  is  added,  and  then  to  an  electron  spin  resonance 
(esr)  cavity  for  atom  detection  and  a quadrupole 
mass  spectrometer  inlet  for  stable  species  detection. 
The  position  of  the  movable  injector  determines  the 
reaction  time  and  allows  the  rate  constant  of  the  reac- 
tion to  be  measured.  The  reactor  can  be  uniformly 
heated  or  cooled  over  a 200  to  1000°K  range. 


MOVABLE 

INJECTOR 


DISCHARGE 

CAVITY 


TO 

PUMP 


HEATER  OR  COOLER 


ESR 

CAVITY 


TO 

MASS 

SPECTROMETER 


DISCUSSION 

For  the  present  reactions,  it  was  desirable  to 
operate  the  apparatus  with  the  concentration  of  H or 
O atoms  in  large  excess  over  that  of  the  CH,X  added 
through  the  movable  injector.  The  mass  spectrometer 
was  then  used  to  monitor  the  decay  of  CH,X  as  a 
function  of  reaction  time  with  an  essentially  constant 
atom  concentration,  thereby  considerably  simplifying 
the  data  analysis.  Results  for  the  rate  constant  of  the 
reaction  H -f  CH,  X -»  HX  -I-  CH,  are  presented  in 
Arrhenius  form  in  Fig.  2 as  the  plotted  lines  through 
the  experimental  points  (squares  and  circles).  The 
differences  in  reactivity  of  CHaBr,  CH,C1,  and  CHjF 
are  immediately  apparent  and  reflect  the  increasing 
C-X  bond  energy  in  that  series.  The  greater  effective- 
ness of  CHjBr  as  a flame  inhibitor  is  also  clearly  indi- 
cated, since  its  appreciably  larger  rate  constant  permits 
the  much  more  efficient  removal  of  H atoms  from  the 
reaction  zone.  A subsequent  fast  step  H HBr  — > 
H,  + Br  means  that  two  H atoms  are  effectively  re- 
combined to  H^  with  the  substitution  of  relatively 
unreactive  CH.,  and  Br.  The  tendency  for  the  plots  to 
curve  somewhat  upward  at  the  h/gh-temperafure  end 
may  be  an  indication  that  the  alternate  reaction  path 
H + CH,X  -p  CH4  -I-  X begins  to  contribute  slightly 
under  these  conditions. 

The  triangular  points  in  Fig.  2 are  data  obtained 
directly  from  the  addition  of  CH,Br  as  a trace  flow 
from  a "point"  source  in  an  actual  Hj-O,  flame  at 
temperatures  in  the  900  to  1100°K  range.  These  data 
were  from  an  entirely  different  experiment  performed 
at  APL  (Ref.  1).  The  flame  experiment  was  carried 
out  using  a mixture  that  was  very  H,  rich,  so  that  the 
concentrations  of  O and  OH  were  greatly  suppressed 
compared  with  H,  and  the  contributions  of  O -b 
CH,X  and  OH  + CH,X  were  minimized.  The  agree- 
ment within  a factor  of  2 or  3 between  the  two  experi- 
ments on  H -t-  CHjBr  in  the  overlapping  tempera- 
ture range  is  satisfactory  and  tends  to  corrolx)rate  some 
of  the  assumptions  made  in  the  flame  experiment. 
Even  more  striking  is  the  agreement  shown  in  Fig.  3, 
where  data  obtained  as  ratios  of  rate  constants  are  pre- 
sented from  the  two  types  of  experiment.  As  i^io 
data,  certain  inherent  errors  in  both  techniques  tend 
to  cancel  out,  and  the  agreement  is  excellent.  This  is 
one  of  the  very  rare  examples  of  a direct  check  on 
kinetic  data  obtained  in  an  actual  flame  by  an  inde- 
pendent method  at  the  same  temperature. 

Similar  experiments  were  carried  out  on  the  reaction 
O + CH,X  -*  CH,X  + OH  where  X is  Br  or  O. 
The  O atom  reactions  differ  fundamentally  from  the 
foregoing  H reactions  since  the  O reactions  involve 
the  abstraction  of  H from  CH,X  while  the  H + 
CH,X  reactions  are  halogen  (X)  abstractions.  Thus 
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Fi*.  1 Schematic  of  fast  flow  teaaor. 
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Fig.  2 Arrhenius  plot  of  rate  constants  for  H + CHsX. 
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Fig.  4 Arrhenius  plot  of  rate  constants  for  O + CHsBr. 
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radical.  Results  for  the  O + CHsBr  rate  constant  are 
plotted  as  the  solid  line  in  Fig.  4.  Also  shown  (the 
dashed  line)  are  earlier  data  from  APL  (Ref.  2)  on 
the  reaction  O + CH4  -*  CH,  + OH.  This  is  also  a 
simple  H abstraction  step,  and  the  similarity  in  slope 
(i.e.,  the  activation  energy)  in  the  two  cases  is  clear. 

Results  for  both  H + CHjX  and  O + CH,X  re- 
actions have  been  published  (Refs.  3 and  4). 
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ATTENUATION  OF  SOUND  BY  RIGID  SPHERES 


An  acoustic  cavity  technique  was  used  to  measure 
separately  the  viscous  and  thermal  components  of  the 
attenuation  of  sound  by  rigid  spheres.  Experiments  on 
spheres  of  known  size  in  various  gases  provided,  for 
the  first  time,  a rigorous  quantitative  test  of  theoreti- 
cal predictions  of  the  individual  attenuation  compo- 
nents. Good  agreement  was  obtained  between  theory 
and  experiment. 


BACKGROUND 

The  absorption  of  sound  by  particles  suspended  in 
a gas  is  an  interesting  process  that  plays  an  important 
role  in  determining  the  stability  of  solid-propellant 
rocket  motors.  Theory  (Refs.  1 to  3)  predicts  an 
acoustic  attenuation  that  can  be  resolved  into  two 
terms:  a viscous  component  caused  by  the  relative 
motion  of  the  gas  with  respect  to  the  particles,  and 
a thermal  component  caused  by  irreversible  heat 
transfer  between  the  gas  and  the  particles. 

Despite  extensive  previous  experimental  studies  on 
the  attenuation  of  sound  by  particles  of  various  sizes 
in  suspensions  and  emulsions,  a rigorous  quantitative 
verification  of  the  theoretical  predictions  for  the  indi- 
vidual components  of  attenuation  had  not  been  ob- 
tained. Under  typical  experimental  conditions  the 
particles  are  subjected  simultaneously  to  both  viscous 
and  thermal  damping  effects,  so  that  the  individual 
effects  are  not  separable.  In  the  experiments  reported 
here  (Ref.  4),  spheres  of  known  size  were  suspended 
in  an  acoustic  cavity  so  that  thq'  could  be  selectively 
subjected  to  either  viscous  or  thermal  damping, 
thereby  permitting  the  two  components  of  attenuation 
to  be  measured  independently. 


DISCUSSION 

The  principle  of  the  method  can  be  understood  by 
considering  the  standing  wave  patterns  (Fig.  1)  for 
a cylindrical  acoustic  cavity  excited  in  the  funda- 
mental and  second  harmonic  longitudinal  modes.  At 
positions  corresponding  to  velocity  antiriodes,  particles 
would  be  subjected  only  to  viscous  damping,  while 
at  velocity  nodes  (pressure  antinodes),  th^  would  be 
subjected  only  to  thermal  damping.  For  example,  if 
a particle  were  translated  along  the  axis  of  the  cavity 
excited  in  the  second  harmonic  mode,  it  would  be 
subjected  to  maximum  thermal  damping  at  X >■  0, 
L/2,  and  L;  to  maximum  viscous  damping  at 


PRESSURE  AMPLITUDE 

VELOCITY  AMPLITUDE 


(a)  FUNDAMENTAL  MODE 


Fig.  1 Acoustic  cavity  excited  in  longitudinal  modes. 


X — Z./4  and  3L/4;  and  to  readily  calculable  mix- 
tures of  viscous  and  thermal  damping  at  other  posi- 
tions. 

The  acoustic  cavity  arrangement  is  shown  sche- 
matically in  Fig.  2.  Standing  waves  are  generated  by 
a piston -type  vibration  generator.  The  small  change 
in  cavity  damping  produced  by  one  or  more  spheres 
suspended  by  the  translating  rod  is  measured.  To  do 
this,  the  acoustic  driver  (vibration  generator)  is  pro- 
grammed to  excite  a particular  mode  for  a prescribed 
time,  typically  10  ms,  allowing  the  energy  in  that 
mode  to  build  up.  The  driver  is  then  turned  off, 
causing  an  exponential  free  decay  of  the  acoustic 
pressure,  P,  in  accordance  with  the  expression 
P — Pt,e~”\  where  P„  is  the  pressure  at  time  / — 0, 
and  a is  the  decay  constant  in  nepers  per  second 
(Np/s).  In  order  to  obtain  accurate  measurements 
of  the  decay  constant,  a special  electronic  system  had 
to  be  designetl  (Ref.  4).  With  this  system,  changes 
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Fig.  2 Acoustic  cavity  arraagement  for  atteouatioo 
measurements. 


in  the  decay  constant  of  about  0.002  Np/s,  repre- 
senting about  0.01%  of  the  empty  cavity  damping, 
could  be  readily  observed. 

Experiments  were  carried  out  on  rigid  spheres 
having  diameters  of  0.794  cm  (3/l6  in.)  and  0.953 
cm  (3/8  in.)  in  argon,  nitrogen,  and  dichlorodifluo- 
romethane  (Freon-12).  These  gases  were  selected  in 
order  to  cover  a wide  range  in  the  relative  influence 
of  viscous  to  thermal  damping. 

Experimental  data  for  the  damping  produced  by  a 
single  0.953— cm-diameter  ball  bearing  in  argon  are 
shown  in  Fig.  3.  The  cavity  was  excited  in  the  second 
harmonic  mode  at  a frequency  of  1.1195  kHz.  Curve 
A shows  the  cavity  damping  produced  by  the  ball 
bearing  and  the  translating  rod,  curve  B shows  the 
damping  produced  by  the  translating  rod  alone,  and 
curve  C shows  the  difference  between  curves  A and 
B and  gives  the  attenuation  produced  by  the  ball 
bearing  alone.  The  solid  curve  is  the  attenuation  pre- 
dicted by  theory.  The  data  are  in  good  quantitative 
agreement  with  theoretical  calculations  for  both  vis- 


Fig.  3 Experimental  data  for  the  damping  of  a single 
0.953-cm-diameter  sphere  in  argon. 


Fig.  4 Comparison  of  experiment  with  theory  for  the 
difference  between  viscous  and  thermal  atteouatioo  for 
spheres  in  various  gases. 
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cous  and  thennal  attenuations.  The  sinusoidal  varia- 
tion of  attenuation  with  distance  is  as  expected,  with 
maxima  at  velocity  antinodes  and  minima  at  velocity 
nodes. 

The  difference  between  the  viscous  and  thermal 
components  of  attenuation  in  be  measured  with 
good  precision,  without  requiring  the  cavity  to  be 
opened.  Results  for  spheres  in  nitrogen,  argon,  and 
Freon-12  are  given  in  Fig.  4.  Figure  4a  shows  the 
variation  in  attenuation  with  displacement  for  three 
0.953^-cm  ball  bearings  in  nitrogen.  Fig.  4b  for  a 
single  0.953-cm  ball  bearing  in  argon,  and  Fig.  4c 
for  three  0.953-cm  ball  bearings  in  Freon-12.  The 
theoretically  calculated  values  in  each  case  are  shown 
as  solid  curves.  The  experiments  cover  a wide  range 
in  the  relative  importance  of  thermal  to  viscous  damp- 
ing, ranging  from  0.544  for  argon,  to  0.322  for  nitro- 
gen, to  0.0975  for  Freon-12.  The  agreement  between 


theory  and  experiment  is  very  good  and  provides 
quantitative  verification  of  the  theoretical  analysis. 


REFERENCES 

1.  P.  S.  Epstein  and  R.  R.  Carhart,  "The  Absorption  of  Sound  in 
Suspensions  and  Emulsions.  I.  water  Fog  in  Air,'*  /.  Acoust. 
Soc.  Am..  25,  195),  pp.  553-565. 

2.  ).  C.  F.  Chow,  '‘Attenuation  of  Acoustic  Waves  in  Dilute  Emul- 
sions and  Suspensions,*'  /.  Acoust.  Soc.  Am.,  36,  pp. 

2395-2401. 

3.  S.  Temkin  and  R.  A.  Dobbins.  "Attenuation  and  DisMrsion  of 
Sound  by  Particulate  Relaxation  Processes,"  /.  Acoust.  Soc.  Am., 
40.  1966.  pp.  317-324. 

4.  S.  N.  Foner  and  B.  H.  Nall,  "Attenuation  of  Sound  by  Rigid 
Spheres:  Measurentent  of  the  Viscous  and  Thermal  Components 
of  Attenuation  and  Comparison  with  Theory,"  /.  Acoust.  Soc. 
Am.,  57,  1975,  pp.  59-66. 


Authors:  S.  N.  Foner  and  B.  H.  Nall 
Support:  NA  VSEA  SYSCOM 


MOLECULAR  ENERGIES  AND  STRUCTURE 


The  diagrammatic  many-body  perturbation  theory 
has  been  formulated  through  third  order  and  appiied 
to  the  Calculation  of  electron  pair  correiation  energies 
of  several  diatomic  molecules.  The  high  accuracy 
achieved  indicates  that  this  technique  shouid  be  useful 
for  the  determination  of  other  correlation-dependent 
properties  such  as  electronic  spectra,  photoionization 
potentials,  polarizabilities,  and  susceptibilities. 


BACKGROUND 

Chemical  reactions  play  an  important,  fundamental, 
and  pervading  role  in  nature  and  in  man’s  use  of 
nature  in  such  areas  as  fuels,  synthetic  materials, 
chemical  lasers,  and  interactions  of  chemicals  in  the 
atmcephere  A deeper  undentanding  of  the  intimate 
letaih  of  the  underlying  chemical  reactions  would 
maMe  muih  pmgreia  to  be  nude  in  terms  of  energy 
unproved  materials,  and  redesign  of 
optiiniae  certain  desirable  prod- 
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-d  m anmaMr  hemiiai  wiiennediale  is  re- 
thnwgh  rapenmental 
H.  levliup  theoretical 
TV  ihsnrr 


the  methods  can  be  extended  systematically  to  more 
complex  and  difficult  problems.  The  present  work 
indicates  that  the  diagrammatic  many-body  perturba- 
tion theory  offers  a very  promising  approach  to  these 
chemical  problems. 


DISCUSSION 

The  many-body  perturbation  theory  had  been  de- 
veloped originally  for  nuclear  physics.  Later  it  was 
used  successfully  to  describe  the  electronic  structure  of 
atoms.  Current  work  (Refs.  1 through  4)  brings 
these  methods  to  bear  on  molecular  systems.  The 
essential  difficulty  with  molecules  is  that  the  electrons 
are  distributed  among  several  nuclear  centers.  This 
problem  has  been  overcome  in  the  present  work  by 
the  use  of  basis  sets  of  Slater-type  atomic  orbitals 
located  on  all  nuclear  centers. 

Preliminary  applications  have  been  made  to  the 
calculation  of  pair-correlation  energies  in  the  dia- 
tomic hydrides  LiH,  BH,  and  HF  (Ref.  1)  and 
NaH  (Ref.  2),  using  small  tusis  sets.  More  ex- 
tensive calculations  have  been  carried  out  (Refs.  3 
and  4)  using  quite  large  basis  sets  consisting  of  46 
irhftali  for  the  closed  shell  first-row  diatomic  hy- 


TVw 


drides.  Even  though  the  first-order  perturbation  ex- 
pansion yields  99.62%  of  the  energy  of  the  HF 
molecule,  this  energy  is  so  large  that  the  remaining 
0.38%,  called  the  correlation  energy,  represents  an 
error  of  over  225  kcal/mole,  which  is  orders  of  mag- 
nitude greater  than  the  tolerable  level  requited  for 
"chemical  accuracy"  (about  1 to  10  kcal/mole).  To 
obtain  this  degree  of  accuracy,  higher  orders  in 
the  perturbation  expansion  need  to  be  considered 
(Ref.  5). 

In  the  many-body  approach,  the  correlation  energy 
is  developed  in  terms  of  diagrams  representing  vari- 
ous types  of  interelectronic  interaction.  All  the  rele- 
vant diagrams  are  presented  in  Fig.  1.  The  second- 
order  correlations  are  derived  from  diagrams  A and 
B in'  the  figure.  Diagrams  C and  D are  of  the  third- 
order  "particle-ladder”  type,  and  diagrams  E and  F 
are  the  corresponding  third-order  "hole-ladder" 
terms.  The  third-order  "ring"  diagrams  G through 
N represent  "hole-particle”  interactions.  Finally  the 
third-order  diagrams  Q through  T represent  the  inter- 
action between  excited  state  orbitals  and  an  arbitrary 
potential  (Ref.  6). 
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I Fig.  I Ail  second-order  (A  and  B)  and  third-order 

i (C  through  T)  energy  diagrams  contributing  to  the  cor- 

f relation  energy  of  molecular  systems  in  the  many-body 


perturbation  theory. 

The  perturbation  expansions  make  use  of  all  pos- 
sible N-electfon  determinantal  wave  functions  that 
can  be  constructed  from  a given  basis  set.  Through 
third  order,  this  collection  of  excited-state  determi- 
nants consists  of  double  excitations  from  the  reference 


state,  i.e.,  determinants  formed  from  the  reference 
determinant  by  deleting  orbitals  i and  / and  replac- 
ing them  with  orbitals  a and  b,  for  all  unique  com- 
binations of  the  indexes  /,  a,  and  b.  The  total  pos- 
sible number  of  these  double  excitations  arising  from 
the  46  orbital  basis  sets  (Refs.  3 and  4)  is  22  968 
for  LiH,  54  825  for  BH,  and  149  445  for  HF.  How- 
ever, because  of  the  spin  and  spatial  symmetry  ortho- 
gonalities of  the  molecular  orbitals,  the  number  of 
determinants  entering  the  perturbation  expansions 
with  nonzero  interactions  is  reduced  to  2396  for  LiH, 
5355  for  BH,  and  14  377  for  HF.  One  of  the  power- 
ful attributes  of  the  diagrammatic  perturbation  method 
(Ref.  5)  is  that  all  of  these  nonzero  interactions  can 
be  efficiently  taken  into  account.  This  is  in  striking 
contrast  to  other  methods  for  handling  the  electron 
correlation  problem,  where  it  is  feasible  to  include  the 
effects  of  only  a small  subset  of  these  interactions. 

All  of  the  nonzero  determinantal  states  hive  been 
used  to  calculate  the  various  correlation  energy  cor- 
rections indicated  in  Fig.  1.  For  each  of  the  closed- 
shell  first-row  diatomic  hydrides,  this  work  recovers 
more  than  94%  of  the  experimental  correlation  energy 
when  the  diagrams  of  Fig.  1 are  calculated  with  pair 
restrictions  and  include  selections  of  higher-order  dia- 
grams via  denominator  shifts.  For  the  HF  molecule 
in  particular,  97%  of  the  correlation  energy  is  recov- 
ered, representing  99.99%  of  the  total  energy  with  a 
residual  error  of  only  about  6 kcal/mole.  For  the  LiH 
and  BH  molecules,  the  residual  errors  are  about  3 and 
5 kcal/mole,  respectively. 

The  application  of  the  many-body  perturbation 
theory  has  been  able  to  produce  energy  results  that  are 
within  chemical  accuracy  for  these  systems.  This  woik 
provides  the  most  accurate  numerical  determination  of 
correlation  energies  ever  obtained  for  these  systems. 
Moreover,  the  results  are  amenable  to  a pair-by-pair 
analysis,  which  permits  an  understanding  of  the 
various  interaction  processes. 
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MASTER’S  DEGREES  AWARDED  AT 
APL  EVENING  COLLEGE  CENTER 
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As  early  as  1958,  a few  technical  courses,  but  no 
degree  programs,  were  offered  after  regular  working 
hours  at  the  Applied  Physics  Laboratory.  These  courses 
were  taught  by  APL  staff  members  under  sponsorship 
of  the  Evening  0>llege  of  The  Johns  Hopkins  Univer- 
sity. In  1963,  as  a result  of  increasing  staff  interest  in 
expanding  the  advanced  education  opportunities  at 
APL  and  the  rapidly  increasing  number  of  persons  in 
the  adjacent  area  desiring  advanced  technical  courses, 
a Center  of  the  Evening  College  of  the  University  was 
authorized  at  APL. 

In  the  fall  of  1964,  the  APL  Evening  College 
Center  formally  began  operation.  A program  leading 
to  the  Master  of  Science  in  Electrical  Engineering  was 
instituted,  and  all  courses  were  opened  to  the  public. 
One  hundred  and  sixteen  individuals  complete  149 
registrations  for  the  six  courses  offered.  (Some  in- 
dividuals registered  for  more  than  one  course.)  In 
1965,  the  courses  offered  increased  to  11,  the  number 
of  individuals  enrolled  to  193,  and  class  registration 
to  241. 

Since  1965,  ad  hoc  committees  of  APL  scientists 
and  engineers  have  from  time  to  time  proposed  and 
helped  to  implement  four  additional  master's  degree 
programs  to  fill  perceived  needs  of  the  APL  staff  and 
die  neighboring  technical  community.  These  programs 
lead  to  Master  of  Science  degrees  with  majors  in  nu- 
merical science  (first  offered  in  1966),  applied  physics 
(1967),  space  technology  (1967),  and  computer  sci- 
ence (1971). 

Class  enrollment  grew  steadily  to  about  600,  involv- 
ing more  than  400  individuals  by  1967.  Both  totals 
remained  nearly  constant  through  1974,  but  1975 


brought  a new  record  high  class  enrollment  of  696, 
with  494  individual  students. 

In  June  1968,  The  Johns  Hopkins  University 
granted  the  first  degrees  awarded  for  the  completion 
of  programs  offered  by  the  Evening  College  Center. 
Since  then,  the  University  has  granted  a total  of  577 
master’s  degrees  to  students  of  the  Center.  These  de- 
grees constitute  71%  of  all  master’s  degrees  in  tech- 
nical areas  granted  during  that  period  by  the  Evening 
College  of  The  Johns  Hopkins  University  and  cover 
five  academic  areas  of  specialization; 

Numerical  science  249, 

Electrical  engineering  224, 

Computer  science  60, 

Applied  physics  23,  and 
Space  technology  21. 

The  APL  Evening  College  Center  program  serves 
not  only  APL  staff  members  desiring  advanced 
scientific  education  but  also  the  entire  technical  com- 
munity of  the  Washington-Baltimore  area,  including 
nearby  Virginia.  The  ratio  of  outside  participants  to 
APL  staff  members  enrolled  has  risen  steadily  and 
86%  of  the  students  now  come  from  outside  APL. 

The  Center  is  especially  fortunate  in  its  faculty. 
Most  of  the  instructors  are  highly  qualified  members 
of  the  APL  Senior  Staff;  thus  the  students  benefit  by 
having  research-minded  teachers  who  are  involved 
with  the  subjects  they  teach  and  understand  the  needs 
of  part-time  students  who  work  in  research  and  de- 
velopment organizations. 


96 
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